Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,543,754 T→C 29.1% T103A (ACG→GCG)  narU ← nitrate/nitrite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,543,7540TC29.1% 166.0 / 98.2 134T103A (ACG→GCG) narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base T (47/48);  new base C (20/19);  total (67/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.24e-01

Only 100 of 134 total aligned reads displayed.

GATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAAGAGTTG  >  NC_000913/1543659‑1543850
                                                                                               |                                                                                                
gataacgataaATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCa                                                                                               >  1:726369/1‑99 (MQ=255)
   aacgataaATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCgg                                                                                           <  1:621831/100‑1 (MQ=255)
   aacgataaATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCgg                                                                                           <  1:1261500/100‑1 (MQ=255)
    acgataaATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGc                                                                                          >  2:413861/1‑100 (MQ=255)
     cgataaATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCg                                                                                         <  2:2623861/100‑1 (MQ=255)
     cgataaATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCg                                                                                         >  2:3974219/1‑100 (MQ=255)
     cgataaATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCCAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCgg                                                                                           <  1:3506223/98‑1 (MQ=255)
       ataaATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTc                                                                                       >  1:199802/1‑100 (MQ=255)
        taaATATCCTAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTcc                                                                                      <  2:2216832/100‑1 (MQ=255)
         aaaTATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCa                                                                                     <  2:1056428/100‑1 (MQ=255)
           atatCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCAcc                                                                                   <  1:905721/100‑1 (MQ=255)
            tatCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCg                                                                                  <  2:4312942/100‑1 (MQ=255)
            tatCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCg                                                                                  <  2:3654849/100‑1 (MQ=255)
              tCCCAAAAGGAGTATTCGGATTTTTCACGGCAATTCGGAGTGAGACGCAAGCAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCAcc                                                                                   <  1:1067895/97‑1 (MQ=255)
               ccgaaaaGGAGTATTTGGATTTTGCACGGCAATTCCGAGCCATACGCAAGGAATAATCAGCATTCCAGTACTAAAAACCGTCCATCGCTGTCCACCGAat                                                                               >  1:1916470/4‑100 (MQ=37)
                ccAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCCAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAata                                                                              >  2:1795402/1‑100 (MQ=255)
                ccAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTTCACCGCata                                                                              >  2:3084319/1‑100 (MQ=255)
                  aaaaGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAAtata                                                                            <  1:3975645/100‑1 (MQ=255)
                      ggAGTATTCGGATTTTGCACGGCAATTCCGAGCGAGACGCAAGGAGTAATCAGGATTGCAGTACTAAAAGCCGTCCATCGGCGTCCACCGAATATAGGCa                                                                        <  1:3755530/100‑1 (MQ=255)
                      ggAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGAGTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCa                                                                        <  1:3306726/100‑1 (MQ=255)
                       gAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCAc                                                                       >  2:2314625/1‑100 (MQ=255)
                        aGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGGCCACCGAATATAGGCAc                                                                       <  1:447244/99‑1 (MQ=255)
                         gTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCa                                                                     >  2:1352033/1‑100 (MQ=255)
                         gTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCa                                                                     >  1:2333524/1‑100 (MQ=255)
                          tATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCAt                                                                    <  1:3831582/100‑1 (MQ=255)
                          tATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCAt                                                                    >  2:3524512/1‑100 (MQ=255)
                          tATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCAt                                                                    >  2:2431587/1‑100 (MQ=255)
                          tATTCGGATTTTGCACGGCAAGTTCGCGCCAGACGCAATGAATAATAACGATTGAAATACTAATAACCGTATAACGGAGTCCACAGaag                                                                               <  1:612748/89‑2 (MQ=255)
                           aTTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATa                                                                   >  1:2835208/1‑100 (MQ=255)
                           aTTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATa                                                                   >  1:3136056/1‑100 (MQ=255)
                            ttCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATaa                                                                  >  2:498297/1‑100 (MQ=255)
                            ttCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATaa                                                                  >  2:3912012/1‑100 (MQ=255)
                            ttCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATaa                                                                  <  2:4332834/100‑1 (MQ=255)
                            ttCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATaa                                                                  >  2:3106219/1‑100 (MQ=255)
                             tCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATaaa                                                                 >  1:2501801/1‑100 (MQ=255)
                             tCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATaaa                                                                 <  2:3707591/100‑1 (MQ=255)
                             tCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATaaa                                                                 <  1:3058490/100‑1 (MQ=255)
                               ggATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAgg                                                               <  1:156627/100‑1 (MQ=255)
                               ggATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAgg                                                               >  2:3907296/1‑100 (MQ=255)
                               ggATTGTGCACGGCAATTCAGAGCGAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAgg                                                               <  1:1225239/100‑1 (MQ=255)
                                 aTTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGCATGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGa                                                              >  1:1952285/1‑99 (MQ=255)
                                 aTTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAg                                                             <  2:3145033/100‑1 (MQ=255)
                                    ttGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAg                                                          <  2:1251015/100‑1 (MQ=255)
                                      gCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAggg                                                        <  2:2215642/100‑1 (MQ=255)
                                         cGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAAGCGTCCATAGGCGTCCACAAAATATAGGCCACATAAAGGAGTAGGGAac                                                     >  1:1053284/1‑100 (MQ=38)
                                         cGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGaa                                                      <  1:646398/99‑1 (MQ=255)
                                         cGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGaa                                                      >  1:888788/1‑99 (MQ=255)
                                          ggCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAaca                                                    <  1:3486741/100‑1 (MQ=255)
                                          ggCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAaca                                                    >  2:3985315/1‑100 (MQ=255)
                                            cAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCAACGAATATTGGCACGATAACGGAGTAGGGAACACg                                                  >  1:2463855/1‑100 (MQ=37)
                                             aaTTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGc                                                 <  2:872897/100‑1 (MQ=255)
                                              aTTCCGAGCCAGACGCACGGAATAATCAGGATTGCAGTACTAAAAACCGTCCAACGGCGTGCACCGAATATACGCACCATAAAGGAGGAGTGAgcac                                                   >  1:957260/1‑97 (MQ=255)
                                              aTTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGc                                                 <  1:4166466/99‑1 (MQ=255)
                                                 ccGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCaata                                             >  1:1265301/1‑100 (MQ=255)
                                                 ccGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCaata                                             >  2:2452749/1‑100 (MQ=255)
                                                 ccGAACCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATAGGCGTCAACCGAATATAGGCACCATAAAGGAGTAGGGAAAACGCaata                                             >  1:1991310/1‑100 (MQ=37)
                                                   gAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCaataat                                           >  2:4372252/1‑100 (MQ=255)
                                                        aGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGcc                                      >  1:3157981/1‑100 (MQ=255)
                                                          aCGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAg                                    <  2:635040/100‑1 (MQ=255)
                                                           cGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCACATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGa                                   >  1:2357409/1‑100 (MQ=255)
                                                             cAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGaaa                                 >  2:1671559/1‑100 (MQ=255)
                                                             cAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGaaa                                 >  2:2282700/1‑100 (MQ=255)
                                                               aGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACg                               <  1:2100452/100‑1 (MQ=255)
                                                                ggAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACgg                              <  1:393626/100‑1 (MQ=255)
                                                                ggAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACgg                              <  2:2501814/100‑1 (MQ=255)
                                                                 gAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGa                             >  2:1370360/1‑100 (MQ=255)
                                                                  aataatCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAg                            >  2:2941774/1‑100 (MQ=255)
                                                                   ataatCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAgg                           <  2:1953360/100‑1 (MQ=255)
                                                                       tCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAggg                          >  1:3761272/1‑97 (MQ=255)
                                                                            aTTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCgg                  <  2:709482/100‑1 (MQ=255)
                                                                            aTTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCgg                  <  2:2936342/100‑1 (MQ=255)
                                                                            aTTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCgg                  <  1:2611498/100‑1 (MQ=255)
                                                                              tGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGtt                <  2:3395075/100‑1 (MQ=255)
                                                                              tGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGt                 <  1:1863101/99‑1 (MQ=255)
                                                                              tGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGtt                <  2:2833291/100‑1 (MQ=255)
                                                                               gCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTta               <  1:4079863/100‑1 (MQ=255)
                                                                               gCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTta               <  2:374603/100‑1 (MQ=255)
                                                                               gCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTta               >  2:4234201/1‑100 (MQ=255)
                                                                                 aGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTtaat             <  2:451519/100‑1 (MQ=255)
                                                                                 aGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTtaat             <  2:3395938/100‑1 (MQ=255)
                                                                                  gTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTtaata            >  2:3274821/1‑100 (MQ=255)
                                                                                  gTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTtaata            >  1:2428040/1‑100 (MQ=255)
                                                                                     cTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATaaaa         <  1:1048529/100‑1 (MQ=255)
                                                                                     cTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATaaaa         >  2:729154/1‑100 (MQ=255)
                                                                                     cTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATaaaa         >  2:14370/1‑100 (MQ=255)
                                                                                      tAAAAACCGTCGATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACTCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATaaaaa        <  1:3287434/100‑1 (MQ=255)
                                                                                      tAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATaaaaa        >  2:778734/1‑100 (MQ=255)
                                                                                       aaaaaCCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAGTGCGCCAGAAACTGCGGGTAATGCGGTTAATCAAAag       >  2:1885284/1‑100 (MQ=255)
                                                                                       aaaaaCCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAag       <  2:2563475/100‑1 (MQ=255)
                                                                                       aaaaaCCGCCGATCGGCGTCCACCAAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAag       <  1:3473508/100‑1 (MQ=255)
                                                                                       aaaaaCCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAag       >  1:772907/1‑100 (MQ=255)
                                                                                       aaaaaCCATCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAag       >  1:3396001/1‑100 (MQ=255)
                                                                                        aaaaCCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAaga      >  1:467944/1‑100 (MQ=255)
                                                                                        aaaaCCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGGAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAaga      >  2:1882516/1‑100 (MQ=255)
                                                                                        aaaaCCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAaga      >  1:3830211/1‑100 (MQ=255)
                                                                                         aaaCCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAagag     <  1:2080795/100‑1 (MQ=255)
                                                                                         aaaCCGCCCATCGTCGTTCAGCGAATATAGGCAGCACAAAGGAGTAGGGAACCCGCAATAATGCGCCAGAAACGGAGGGTAATGCGGtt                >  1:2982918/1‑89 (MQ=255)
                                                                                          aaCCGCCCATCGGCGTCCACCGAATATAGGCACCATCAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAAGAGt    <  1:2431576/100‑1 (MQ=255)
                                                                                            ccGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAAGAGTTg  <  1:3191231/100‑1 (MQ=255)
                                                                                            ccGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGCACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAAGAGTTg  >  1:3786393/1‑100 (MQ=255)
                                                                                               |                                                                                                
GATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAAGAGTTG  >  NC_000913/1543659‑1543850

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: