Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,052,261 G→A 9.4% intergenic (+21/+28) gnsA → / ← yccM putative phosphatidylethanolamine synthesis regulator/putative 4Fe‑4S membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,052,2610GA 9.4% 337.0 / 16.8 116intergenic (+21/+28)gnsA/yccMputative phosphatidylethanolamine synthesis regulator/putative 4Fe‑4S membrane protein
Reads supporting (aligned to +/- strand):  ref base G (59/46);  new base A (6/5);  total (65/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.89e-01

Only 100 of 117 total aligned reads displayed.

GACTTTATCGCGCAAAAAATCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAAC  >  NC_000913/1052112‑1052408
                                                                                                                                                     |                                                                                                                                                   
gACTTTATCGCGCAAAAAATCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTg                                                                                                                                                     <  1:967097/150‑1 (MQ=255)
  cTTTATCGCGCAAAAAATCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCg                                                                                                                                                   >  1:424330/1‑150 (MQ=255)
     tATCGCGCAAAAAATCGCCGAGCTGAACAAGAATACTGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAgg                                                                                                                                                >  2:1625207/1‑150 (MQ=255)
     tATCGCGCAAAAAATCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAgg                                                                                                                                                <  1:1679870/150‑1 (MQ=255)
     tATCGCGCAAAAAATCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCGTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAgg                                                                                                                                                <  1:1211300/150‑1 (MQ=255)
      aTCGCGCAAAAAATCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGt                                                                                                                                               <  1:216093/150‑1 (MQ=255)
       tCGCTCAAAAAATCGCCGAGCTGAACAAGAATACAGTGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTTATTTTGTTAACATCACCCCGTGCGAGGTg                                                                                                                                              <  2:1713388/150‑1 (MQ=255)
       tCGCGCAAAAAATCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTg                                                                                                                                              <  1:564763/150‑1 (MQ=255)
          cgcAAAAAATCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCTCGCACGGGGTGATg                                                                                                                                           <  2:1460286/150‑1 (MQ=255)
          cgcAAAAAATCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATg                                                                                                                                           <  1:1055028/150‑1 (MQ=255)
           gcAAAAAATCGGCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTGGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCTCGCACGGGGTGATGt                                                                                                                                          <  2:1358286/150‑1 (MQ=255)
            cAAAAAATCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCTCGCACGGGGTGATGtt                                                                                                                                         <  2:678670/150‑1 (MQ=255)
            cAAAAAATCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGAGGtt                                                                                                                                         >  1:823249/1‑150 (MQ=255)
             aaaaaaTCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCTCGCACGGGGTGATGTTc                                                                                                                                        <  1:755415/150‑1 (MQ=255)
              aaaacTCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCTCGCACGGGGTGATGTTca                                                                                                                                       <  1:1112248/150‑2 (MQ=255)
                aaaTCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGc                                                                                                                                     <  1:307940/150‑1 (MQ=255)
                  aTCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCtt                                                                                                                                   >  2:188749/1‑150 (MQ=255)
                     gCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCATCCCGTGCGAGGTGATGTTCCGCttgtt                                                                                                                                >  2:546943/1‑150 (MQ=255)
                     gCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCttgat                                                                                                                                >  2:273949/1‑148 (MQ=255)
                      ccGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCttgttg                                                                                                                               <  1:908677/150‑1 (MQ=255)
                            tGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAAtt                                                                                                                         >  1:603271/1‑150 (MQ=255)
                             gAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTTTGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTAGTTGCTAAttt                                                                                                                        >  2:1089051/1‑150 (MQ=255)
                             gAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAAttt                                                                                                                        <  2:1328193/150‑1 (MQ=255)
                             gAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACACCGTGCGAGGTGATGTTCCGCTTGTTGCTAAttt                                                                                                                        >  2:73052/1‑150 (MQ=255)
                                cAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGCCCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCTCTTGTTGCTAATTTAGt                                                                                                                     >  2:887832/1‑150 (MQ=255)
                                cAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATTATGTCAAAATCAAAATAATGTCATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGt                                                                                                                     >  2:1187137/1‑150 (MQ=255)
                                cAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGt                                                                                                                     >  1:554156/1‑150 (MQ=255)
                                 aaGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTg                                                                                                                    <  1:754067/150‑1 (MQ=255)
                                 aaGAATACAGGGAAAGAAGTCGCTGAAATTCGCTTCACCGCACGAGAAAAAATGCCCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTg                                                                                                                    <  1:827692/150‑1 (MQ=255)
                                  aGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGa                                                                                                                   >  2:1104029/1‑150 (MQ=255)
                                   gAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGAc                                                                                                                  >  1:817549/1‑150 (MQ=255)
                                   gAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGAc                                                                                                                  <  2:1392847/150‑1 (MQ=255)
                                    aaTACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGAcc                                                                                                                 <  1:1104029/150‑1 (MQ=255)
                                        cAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATc                                                                                                             >  1:559606/1‑150 (MQ=255)
                                        cAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATc                                                                                                             <  1:1220732/150‑1 (MQ=255)
                                        aaGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACTATAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATc                                                                                                             <  1:340559/149‑1 (MQ=255)
                                           ggAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCAtt                                                                                                          >  1:1437993/1‑150 (MQ=255)
                                                gaagTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGgcgc                                                                                                     >  1:41055/1‑150 (MQ=255)
                                                gaagTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGgcgc                                                                                                     >  2:117653/1‑150 (MQ=255)
                                                 aagGCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGCCCGGGCTGGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCt                                                                                                    <  1:1404794/150‑1 (MQ=255)
                                                  agTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATTATGTGATTTTGTGAACATCACCACGTGCGAGGTGATCTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCtt                                                                                                   >  2:603181/1‑150 (MQ=255)
                                                  agTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCtt                                                                                                   >  2:848169/1‑150 (MQ=255)
                                                  agTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCGTGTTGCTAATTTAGTGACCAATCATTGGCGCtt                                                                                                   >  2:1270570/1‑150 (MQ=255)
                                                     ctctGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTtgtg                                                                                                >  1:220682/1‑150 (MQ=255)
                                                      tctGATATTCGCTTCTACGCACGAGAAAATATGTCCGGGCTTGAAAGTTATTATGTTAATATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTgg                                                                                               <  1:784404/150‑1 (MQ=255)
                                                       ctGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGa                                                                                              >  1:1277318/1‑150 (MQ=255)
                                                       ctGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGa                                                                                              >  1:88369/1‑150 (MQ=255)
                                                         gAAATTCGCTTCACCGCACGAGAAAATATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAAt                                                                                            <  1:1586894/150‑1 (MQ=255)
                                                              tCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGc                                                                                       <  1:1411997/150‑1 (MQ=255)
                                                               cGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCg                                                                                      <  2:424330/150‑1 (MQ=255)
                                                                 cTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTc                                                                                    <  2:1516311/150‑1 (MQ=255)
                                                                        gCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGGTAATTTAGTGACCAATCATTGGCGCTTGAGGAATTAAGCGTCGGTACaa                                                                             >  2:753362/1‑150 (MQ=255)
                                                                         cACGAGAAAAAATGACCGGGCTTGAAAGTAATGATGTCAAAATCAAAATAATGAGAATTTCTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCAATGCAGCTAGTGGAATTAAGAGTCGGTACaaa                                                                            >  2:521204/1‑149 (MQ=255)
                                                                             agaAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCAGCTTGTTGCTAATTTAGTGACCAATCATTGGCACTTGAGGAATTAAGCGTCGGTACAATtcct                                                                        >  1:1686495/1‑150 (MQ=255)
                                                                              gaAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATtcctc                                                                       <  2:744230/150‑1 (MQ=255)
                                                                               aaaaaaTGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATtcctcc                                                                      >  1:917558/1‑150 (MQ=255)
                                                                               aaaaaaTGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATtcctcc                                                                      <  1:361484/150‑1 (MQ=255)
                                                                                  aaaTGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGc                                                                   <  1:1625207/150‑1 (MQ=255)
                                                                                       aCCGGGCGTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGCCCAATCATTGGCGCTTGTGGAATTAAGCTGCGGTACAATTCCTCCGGCACCgg                                                              >  1:163859/1‑150 (MQ=255)
                                                                                            gCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCttt                                                         <  1:1345994/150‑1 (MQ=255)
                                                                                             ctggAAAGTTATGATGTCAAAATGATAATAATGTGATTTTGTTAACATCACCCCGTGCGTGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAAGCATTGGCGCTTGTGGAATTCAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTg                                                        <  1:209601/147‑1 (MQ=255)
                                                                                             cTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTACTCCGGCACCGGGCTTTg                                                        >  2:1671863/1‑150 (MQ=255)
                                                                                               tGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGcc                                                      >  1:738896/1‑150 (MQ=255)
                                                                                               tGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGcc                                                      <  2:1109129/150‑1 (MQ=255)
                                                                                                gAAAGTTATGATGACAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAAGTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCa                                                     <  1:1260548/150‑1 (MQ=255)
                                                                                                 aaaGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCAt                                                    >  2:1239198/1‑150 (MQ=255)
                                                                                                      tATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTcc                                               <  2:63946/150‑1 (MQ=255)
                                                                                                      tATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACACCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTcc                                               >  1:651904/1‑150 (MQ=255)
                                                                                                         gatgTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGc                                            <  1:73052/150‑1 (MQ=255)
                                                                                                          atgTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCa                                           >  2:314443/1‑150 (MQ=255)
                                                                                                           tgTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATCCTCCCGCAt                                          >  2:87108/1‑150 (MQ=255)
                                                                                                            gTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATa                                         >  1:714348/1‑150 (MQ=255)
                                                                                                                 aatcaaaatAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTg                                    >  1:1399982/1‑150 (MQ=255)
                                                                                                                 aatcaaaatAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTg                                    >  1:1686604/1‑150 (MQ=255)
                                                                                                                 aatcaaaatAATGTGATTTTGTGAACATCACCCCGTGCGAGGGGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTg                                    <  1:1313989/150‑1 (MQ=255)
                                                                                                                    caaaatAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGaa                                                                                             <  1:1494013/90‑1 (MQ=255)
                                                                                                                    caaaatAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGaa                                                                                             >  2:1494013/1‑90 (MQ=255)
                                                                                                                    aaaaatAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGt                                 <  1:1706507/149‑1 (MQ=255)
                                                                                                                          aatGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCaa                           <  2:1726258/150‑1 (MQ=255)
                                                                                                                            tgtgATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATc                         <  1:520835/150‑1 (MQ=255)
                                                                                                                            tgtgATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTACTGACCAATCATTGGAGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTACCCAATc                         >  2:467587/1‑150 (MQ=255)
                                                                                                                                     gtgAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTTACCAATCACCGCAAAc                >  2:184365/1‑150 (MQ=255)
                                                                                                                                     gtgAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAAc                >  1:807261/1‑150 (MQ=255)
                                                                                                                                      tgAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACa               >  1:186857/1‑150 (MQ=255)
                                                                                                                                      tgAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACa               >  1:859512/1‑150 (MQ=255)
                                                                                                                                      ggAACATCACCTCGCACGGGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACa               >  1:965431/2‑150 (MQ=255)
                                                                                                                                      ggAACATCACCTCGCACGGGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGAGATTGTGGAATTAAGCGTCGGTACAAATCCTCCGGCACCGCGCTTTGACATACTCCCGCATACATTGCGTAACCAATCACCGGAAACa               >  2:1513383/2‑150 (MQ=255)
                                                                                                                                        aaCATCACCTCGCACGGGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATTCATTGCGTAACCAATCACCGCAAACATa             >  2:659534/1‑150 (MQ=255)
                                                                                                                                         aCATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATaa            <  1:1187137/150‑1 (MQ=255)
                                                                                                                                         aCATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACAAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATaa            >  1:1176560/1‑150 (MQ=255)
                                                                                                                                           aTCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATc          <  1:188749/150‑1 (MQ=255)
                                                                                                                                           aTCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATc          >  1:1281522/1‑150 (MQ=255)
                                                                                                                                            tCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATcc         <  1:62049/150‑1 (MQ=255)
                                                                                                                                            tCACCCCGTGCAAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTACCCAATCACCGCAAACATAATcc         >  2:191093/1‑150 (MQ=255)
                                                                                                                                             cACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGTAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTATCCAATCACCGCAAACATAATccc        >  2:1140387/1‑150 (MQ=255)
                                                                                                                                              aCCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTACCCAATCACGGCACGCATAATcccc       >  1:116551/1‑150 (MQ=255)
                                                                                                                                               ccccGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCa      <  2:603271/150‑1 (MQ=255)
                                                                                                                                                 ccGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGa    >  2:1551073/1‑150 (MQ=255)
                                                                                                                                                   gTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAAc  <  2:559606/150‑1 (MQ=255)
                                                                                                                                                   gTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAAc  >  1:828359/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GACTTTATCGCGCAAAAAATCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAAC  >  NC_000913/1052112‑1052408

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