Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,543,754 T→C 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,7540TC96.6% 387.5 / 0.9 119T103A (ACG→GCG) narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base C (63/51);  total (65/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

Only 100 of 119 total aligned reads displayed.

CGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAAGAGTTG  >  NC_000913/1543658‑1543850
                                                                                                |                                                                                                
cgataacgataaATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCa                                                                                               >  1:2395026/1‑100 (MQ=255)
cgataacgataaATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCa                                                                                               <  2:2778409/100‑1 (MQ=255)
  ataacgataaATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATc                                                                                             >  1:1543150/1‑100 (MQ=255)
      cgataaATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCg                                                                                         >  2:1922523/1‑100 (MQ=255)
      cgataaATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCg                                                                                         <  2:2415312/100‑1 (MQ=255)
      cgataaATATCCCAAAAGGAGTATGCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCg                                                                                         <  2:1664012/100‑1 (MQ=255)
       gataaATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGt                                                                                        >  1:6980/1‑100 (MQ=255)
        ataaATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTc                                                                                       <  1:2707572/100‑1 (MQ=255)
         taaATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTcc                                                                                      >  1:562572/1‑100 (MQ=255)
             tatCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCg                                                                                  <  1:232774/100‑1 (MQ=255)
              atCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGa                                                                                 >  1:1539331/1‑100 (MQ=255)
               tCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGaa                                                                                <  1:2167282/100‑1 (MQ=255)
                cccAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAat                                                                               <  1:1982820/100‑1 (MQ=255)
                cccAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAat                                                                               <  2:2850004/100‑1 (MQ=255)
                cccAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAat                                                                               <  2:2536247/100‑1 (MQ=255)
                cccAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAat                                                                               <  1:2736120/100‑1 (MQ=255)
                 ccAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAata                                                                              >  2:2192707/1‑100 (MQ=255)
                 ccAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAata                                                                              >  2:666403/1‑100 (MQ=255)
                  cAAAAGGAGTATTCGGATTTTGCACGACAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGAAGAACTAACAACCGCCCATCGGCGTCCACCGAat                                                                               >  2:905321/1‑98 (MQ=38)
                   aaaaGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAAtata                                                                            >  2:212371/1‑100 (MQ=255)
                   aaaaGGAGTATCCGGATTTTGCACGGCTATTCTGAGCCAGACGCAAGGAATAAGCAGGATAAAAAAAAAAAAAACCGCCCATCGGCGTCCACCGAAtata                                                                            <  2:1286467/100‑1 (MQ=255)
                     aaGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAgg                                                                          <  2:1355219/100‑1 (MQ=255)
                      aGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGc                                                                         <  2:1919703/100‑1 (MQ=255)
                       ggAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCa                                                                        >  1:989210/1‑100 (MQ=255)
                        gAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCAc                                                                       >  1:638114/1‑100 (MQ=255)
                         aGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCAcc                                                                      <  1:2192720/100‑1 (MQ=255)
                         aGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCAcc                                                                      <  2:1543141/100‑1 (MQ=255)
                           tATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCAt                                                                    >  1:1101884/1‑100 (MQ=255)
                            aTTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATa                                                                   >  1:1830319/1‑100 (MQ=255)
                              tCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATaaa                                                                 <  2:203837/100‑1 (MQ=255)
                                 gATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGa                                                              >  2:927817/1‑100 (MQ=255)
                                 gATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGAAAAAAAAAAAAAAAACCGCCCATCGGCGTCCACCGAATATATGCACCATAAAgg                                                               <  2:2670008/99‑1 (MQ=25)
                                  aTTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAg                                                             <  1:2340878/100‑1 (MQ=255)
                                   ttttGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGt                                                            >  1:1745853/1‑100 (MQ=255)
                                   ttttGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGt                                                            <  2:324086/100‑1 (MQ=255)
                                    tttGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTa                                                           >  1:1172131/1‑100 (MQ=255)
                                     ttGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAg                                                          >  1:1161130/1‑100 (MQ=255)
                                      tGCACGGCAATTCCGAGCCAGACGCAAGGAATAATTAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAgg                                                         <  1:611124/100‑1 (MQ=255)
                                       gCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAggg                                                        <  2:1831782/100‑1 (MQ=255)
                                       gCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAggg                                                        >  1:2256003/1‑100 (MQ=255)
                                       gCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAggg                                                        >  2:2875592/1‑100 (MQ=255)
                                          cGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAac                                                     >  2:1270422/1‑100 (MQ=255)
                                          cGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAac                                                     <  2:699125/100‑1 (MQ=255)
                                          cGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAac                                                     >  1:2537760/1‑100 (MQ=255)
                                               aTTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCa                                                <  2:756129/100‑1 (MQ=255)
                                               aTTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCa                                                <  1:1456456/100‑1 (MQ=255)
                                                  ccGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCaata                                             >  1:1940609/1‑100 (MQ=255)
                                                   cGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCaataa                                            >  1:2486001/1‑100 (MQ=255)
                                                    gAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCaataat                                           >  1:117457/1‑100 (MQ=255)
                                                      gCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATgc                                         >  2:2775724/1‑100 (MQ=255)
                                                      gCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATgc                                         >  1:2409098/1‑100 (MQ=255)
                                                        cAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATgcgc                                       >  1:1582678/1‑100 (MQ=255)
                                                        cAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATgcgc                                       >  2:1541247/1‑100 (MQ=255)
                                                         aGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGcc                                      >  1:2259750/1‑100 (MQ=255)
                                                             gCAAGGAATAATCACGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAg                                    <  2:1790499/98‑1 (MQ=255)
                                                               aaGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAAc                                >  2:2482314/1‑100 (MQ=255)
                                                               aaGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAAc                                >  1:28722/1‑100 (MQ=255)
                                                               aaGGAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAAc                                >  1:2594176/1‑100 (MQ=255)
                                                                 ggAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACgg                              >  1:1375720/1‑100 (MQ=255)
                                                                 ggAATAATCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACCCAATAATGCGCCAGaa                                  >  2:733666/1‑96 (MQ=255)
                                                                  gattaatCAGGATTGCAGTACTAAAAACCACCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAAAACGCAATAATGCGCCAGattcgta                             <  2:937736/97‑7 (MQ=37)
                                                                   aataatCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAg                            <  2:486185/100‑1 (MQ=255)
                                                                   aataatCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAg                            <  2:523714/100‑1 (MQ=255)
                                                                     taatCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAggg                          <  2:412180/100‑1 (MQ=255)
                                                                      aatCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGt                         <  2:2693340/100‑1 (MQ=255)
                                                                      aatCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGt                         >  1:2595539/1‑100 (MQ=255)
                                                                      aatCAGGATTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGt                         <  2:2595520/100‑1 (MQ=255)
                                                                             aTTGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCgg                  <  1:2368296/100‑1 (MQ=255)
                                                                              ttGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGt                 >  2:869108/1‑100 (MQ=255)
                                                                              ttGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGt                 <  1:150257/100‑1 (MQ=255)
                                                                               tGCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGtt                <  2:1599981/100‑1 (MQ=255)
                                                                                gCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTta               <  2:6980/100‑1 (MQ=255)
                                                                                gCAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTta               >  1:1466442/1‑100 (MQ=255)
                                                                                 cAGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTtaa              >  2:2013280/1‑100 (MQ=255)
                                                                                  aGTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTtaat             >  1:1238443/1‑100 (MQ=255)
                                                                                   gTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTtaata            >  1:1566011/1‑100 (MQ=255)
                                                                                   gTACTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTtaata            >  1:1317336/1‑100 (MQ=255)
                                                                                     aCTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATaaa          >  2:593784/1‑100 (MQ=255)
                                                                                     aCTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATaaa          <  1:1130996/100‑1 (MQ=255)
                                                                                     aCTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATaaa          >  2:1047438/1‑100 (MQ=255)
                                                                                     aCTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATaaa          <  1:2482330/100‑1 (MQ=255)
                                                                                      cTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATaaaa         >  1:922470/1‑100 (MQ=255)
                                                                                      cTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATaaaa         >  2:1706548/1‑100 (MQ=255)
                                                                                      cTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATaaaa         <  2:1278578/100‑1 (MQ=255)
                                                                                      cTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATaaaa         >  1:184137/1‑100 (MQ=255)
                                                                                      cTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATaaaa         >  1:185831/1‑100 (MQ=255)
                                                                                      cTAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATaaaa         >  1:2170100/1‑100 (MQ=255)
                                                                                       tAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATaaaaa        <  1:798336/100‑1 (MQ=255)
                                                                                       tAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATaaaaa        >  1:167375/1‑100 (MQ=255)
                                                                                       tAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATaaaaa        <  2:185830/100‑1 (MQ=255)
                                                                                       tAAAAACCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATaaaaa        <  2:184136/100‑1 (MQ=255)
                                                                                       tAAAAACCGCCCATCGGCGTACACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGAGGTTAATaaaaa        <  2:328902/100‑1 (MQ=255)
                                                                                        naaaaCCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATACAAag       >  2:1257968/2‑100 (MQ=255)
                                                                                        naaaaCCGCCCATCGGCGTCCACCGAATACAGGCACCACAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGAGGTTAATAAAAag       >  2:655291/2‑100 (MQ=255)
                                                                                        aaaaaCCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAag       >  1:379168/1‑100 (MQ=255)
                                                                                        aaaaaCCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAag       >  1:331449/1‑100 (MQ=255)
                                                                                        aaaaaCCGACCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAAAGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAag       <  2:429921/100‑1 (MQ=255)
                                                                                         aaaaCCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAaga      <  1:964765/100‑1 (MQ=255)
                                                                                          aaaCCGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGCAGGGTAATGCTGCTAATAAAAagag     <  2:1069613/100‑1 (MQ=37)
                                                                                             ccGCCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAAGAGTTg  >  1:369511/1‑100 (MQ=255)
                                                                                                |                                                                                                
CGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAAGAGTTG  >  NC_000913/1543658‑1543850

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