Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 2,801,061 G→C 5.2% R353P (CGA→CCA)  nrdE → ribonucleoside‑diphosphate reductase 2 subunit alpha

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077792,801,0610G→C 5.2% 30.9 ‑4.5 ‑0.0266R353P (CGA→CCA) nrdEribonucleoside‑diphosphate reductase 2 subunit alpha
Reads supporting (aligned to +/- strand):  new base (2/12):  ref base (127/125):  total (129/137)
Fisher's exact test strand distribution p-value = 1.13e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.98e-01

Only 100 of 281 total aligned reads displayed.

AAGCCGTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGTAAT  >  NC_007779/2800961‑2801161
                                                                                                    |                                                                                                    
aaGCCGTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCg                                                                                                      >  1:13553537/1‑101 (MQ=255)
aaGCCGTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCg                                                                                                      >  1:7777965/1‑101 (MQ=255)
aaGCCGTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCg                                                                                                      <  1:1458835/101‑1 (MQ=255)
 aGCCGTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGa                                                                                                     >  1:4709463/1‑101 (MQ=255)
  gCCGTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGAc                                                                                                    <  1:10858747/101‑1 (MQ=255)
  gCCGTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGAc                                                                                                    >  1:11079140/1‑101 (MQ=255)
   ccGTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATTAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACt                                                                                                   >  1:11826748/1‑101 (MQ=255)
    cGTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTg                                                                                                  >  1:621289/1‑101 (MQ=255)
    cGTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTg                                                                                                  <  1:7411567/101‑1 (MQ=255)
     gTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTgg                                                                                                 >  1:5993441/1‑101 (MQ=255)
     gTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTgg                                                                                                 >  1:8810269/1‑101 (MQ=255)
     gTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTgg                                                                                                 >  1:106167/1‑101 (MQ=255)
      tttGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGc                                                                                                <  1:12385723/101‑1 (MQ=255)
          ccGATGTGGCCATCAGCCAACACTATGACGAACTGTTTGCCGATGAACGCATTCGCAAAAAATCCCTCAACGCCCGTGATTTCTTCAGCGACTGGCAGaaa                                                                                           <  1:207621/101‑1 (MQ=255)
            gATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCACCGCCCGTGATTTCTTCCAGCGACTGGCAGAAAt                                                                                          >  1:3879489/1‑101 (MQ=255)
            gATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAAt                                                                                          >  1:8369237/1‑101 (MQ=255)
             aTGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATc                                                                                         >  1:1634298/1‑101 (MQ=255)
             aTGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATc                                                                                         <  1:13815055/101‑1 (MQ=255)
               gtgGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCa                                                                                       <  1:11150361/101‑1 (MQ=255)
                  gCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGtt                                                                                    >  1:4421575/1‑101 (MQ=255)
                   ccATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGttt                                                                                   <  1:2338985/101‑1 (MQ=255)
                   ccATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGttt                                                                                   <  1:1843029/101‑1 (MQ=255)
                   ccATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGttt                                                                                   <  1:5124975/101‑1 (MQ=255)
                   ccATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGttt                                                                                   >  1:6775772/1‑101 (MQ=255)
                    cATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTg                                                                                  >  1:7084108/1‑101 (MQ=255)
                     aTCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGa                                                                                 >  1:5645293/1‑101 (MQ=255)
                     aTCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGa                                                                                 <  1:2735081/101‑1 (MQ=255)
                      tCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAg                                                                                <  1:12168974/101‑1 (MQ=255)
                      tCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAg                                                                                >  1:4751739/1‑101 (MQ=255)
                           cAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCgg                                                                           <  1:6159645/101‑1 (MQ=255)
                            aaCACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGc                                                                          <  1:11193854/101‑1 (MQ=255)
                             acacTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCt                                                                         <  1:6730903/101‑1 (MQ=255)
                                cTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATc                                                                      <  1:2815374/101‑1 (MQ=255)
                                 tATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATcc                                                                     >  1:4311386/1‑101 (MQ=255)
                                 tATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATcc                                                                     >  1:5484123/1‑101 (MQ=255)
                                       gAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTAcat                                                               >  1:13123089/1‑101 (MQ=255)
                                        aaCTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTAcatc                                                              >  1:7508495/1‑101 (MQ=255)
                                        aaCTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGAAGAAATCCAGTTTGAGTCCGGCTATCCCTAcatc                                                              >  1:13903722/1‑101 (MQ=255)
                                         aCTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTAcatca                                                             >  1:10688330/1‑101 (MQ=255)
                                          cTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTAcatcat                                                            >  1:3751395/1‑101 (MQ=255)
                                          cTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTAcatcat                                                            <  1:6077380/101‑1 (MQ=255)
                                             gTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTTTTTCCGCGACTGGCAGAAAACCAGGTTTAGTCCCGGCATCCCTACATCATGTa                                                         >  1:10969050/1‑101 (MQ=255)
                                              ttGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTCAGTCCGGCTATCCCTACATCATGTAt                                                        <  1:8206832/101‑1 (MQ=255)
                                                 ccGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGaa                                                     <  1:2365246/101‑1 (MQ=255)
                                                  cGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAg                                                    <  1:256769/101‑1 (MQ=255)
                                                  cGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCAATCCCTACATCATGTATGAAg                                                    >  1:9830697/1‑101 (MQ=255)
                                                    aTGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGac                                                  >  1:10174660/1‑101 (MQ=255)
                                                     tGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGaca                                                 >  1:4013839/1‑101 (MQ=255)
                                                      gAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGacac                                                <  1:11463499/101‑1 (MQ=255)
                                                      gAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGacac                                                <  1:9261035/101‑1 (MQ=255)
                                                      gAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGacac                                                <  1:4678379/101‑1 (MQ=255)
                                                       aaCGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACg                                               >  1:13380005/1‑101 (MQ=255)
                                                       aaCGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACg                                               <  1:4170579/101‑1 (MQ=255)
                                                        aCGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACgg                                              >  1:13536094/1‑101 (MQ=255)
                                                            aTTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTaaa                                          <  1:7607085/101‑1 (MQ=255)
                                                            aTTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTaaa                                          <  1:13184624/101‑1 (MQ=255)
                                                               cGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCg                                       >  1:14218751/1‑101 (MQ=255)
                                                               cGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCg                                       >  1:1814786/1‑101 (MQ=255)
                                                                gCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATTAAGACACGGTAAACCGt                                      >  1:10297948/1‑101 (MQ=255)
                                                                gCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGt                                      >  1:4107360/1‑101 (MQ=255)
                                                                   aaaaaTACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCt                                   <  1:9035758/101‑1 (MQ=255)
                                                                   aaaaaTACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCt                                   >  1:328336/1‑101 (MQ=255)
                                                                    gcagaagatCAACGCGCGTCAGTTCTTCCAGACCATCGCTGAGATTCAGTTTGAGTCCGGCTACCCGTACATCATGTTCGAGGACACGGTTAACCGCGcga                                  <  1:4948694/93‑3 (MQ=255)
                                                                    aaaaTACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTa                                  >  1:9308373/1‑101 (MQ=255)
                                                                      agaagatCAACGCGCGTCAGTTCTTCCAGACCATCGCTGAGATTCAGTTTGAGTCCGGCTACCCGTACATCATGTTCGAGGACACGGTTAACCGCGCGAAc                                <  1:1299465/95‑1 (MQ=255)
                                                                       aTACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAAcc                               <  1:8804585/101‑1 (MQ=255)
                                                                        tACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTTAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAAccc                              >  1:6737628/1‑101 (MQ=255)
                                                                        tACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAAccc                              >  1:9514102/1‑101 (MQ=255)
                                                                        aagatCAACGCGCGTCAGTTCTTCCAGACCATCGCTGAGATTCAGTTTGAGTCCGGCTACCCGTACATCATGTTCGAGGACACGGTTAACCGCGCGAAccc                              <  1:6465489/97‑1 (MQ=255)
                                                                           atCAACGCGCGTCAGTTCTTCCAGACCATCGCTGAGATTCAGTTTGAGTCCGGCTACCCGTACATCATGTTCGAGGACACGGTTAACCGCGCGAACCCGAt                           <  1:5512905/100‑1 (MQ=255)
                                                                            tCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATc                          <  1:12864616/101‑1 (MQ=255)
                                                                            tCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATc                          >  1:2487119/1‑101 (MQ=255)
                                                                            tCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATc                          >  1:3069252/1‑101 (MQ=255)
                                                                             cAACGCGCGTCAGTTCTTCCAGACCATCGCTGAGATTCAGTTTGAGTCCGGCTACCCGTACATCATGTTCGAGGACACGGTTAACCGCGCGAACCCGAttg                         >  1:5796433/1‑99 (MQ=255)
                                                                               aCGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGcc                       >  1:12706335/1‑101 (MQ=255)
                                                                               aCGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGcc                       >  1:4357039/1‑101 (MQ=255)
                                                                                cGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTACACCGTGCTAACCCTATCGCCg                      >  1:1908588/1‑101 (MQ=255)
                                                                                 gCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCgg                     >  1:4537508/1‑101 (MQ=255)
                                                                                 gCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCgg                     >  1:6958450/1‑101 (MQ=255)
                                                                                   ccGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGgc                   <  1:11594863/101‑1 (MQ=255)
                                                                                   ccGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGgc                   <  1:12681469/101‑1 (MQ=255)
                                                                                   ccGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGgc                   <  1:12035684/101‑1 (MQ=255)
                                                                                     gTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGgcgc                 <  1:10474198/101‑1 (MQ=255)
                                                                                     gTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGgcgc                 >  1:53996/1‑101 (MQ=255)
                                                                                      tGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCa                >  1:12220525/1‑101 (MQ=255)
                                                                                      tGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCa                >  1:1991949/1‑101 (MQ=255)
                                                                                          ttcttcCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATaaa            <  1:7945697/101‑1 (MQ=255)
                                                                                          ttcttcCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATaaa            <  1:8693577/101‑1 (MQ=255)
                                                                                          ttcttcCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATaaa            <  1:9118365/101‑1 (MQ=255)
                                                                                          ttcttcCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATaaa            <  1:9210849/101‑1 (MQ=255)
                                                                                           tcttcCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATaaaa           <  1:3480332/101‑2 (MQ=255)
                                                                                           tcttcCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAat           >  1:7289024/1‑101 (MQ=255)
                                                                                            cttcCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAata          >  1:10017738/1‑101 (MQ=255)
                                                                                             ttcCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAatat         <  1:49870/101‑1 (MQ=255)
                                                                                               ccAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGa       <  1:8175816/101‑1 (MQ=255)
                                                                                               ccAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGa       >  1:669488/1‑101 (MQ=255)
                                                                                               ccAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCACATAAATATGa       >  1:8550339/1‑101 (MQ=255)
                                                                                                  gCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGCCGCATAAATATGAGTa    <  1:9580412/101‑1 (MQ=255)
                                                                                                    gACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGTAAt  >  1:12280993/1‑101 (MQ=255)
                                                                                                    gACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGTAAt  >  1:8947068/1‑101 (MQ=255)
                                                                                                    |                                                                                                    
AAGCCGTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGTAAT  >  NC_007779/2800961‑2801161

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.