Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 848,362 A→G 6.8% intergenic (‑358/+46) glnH ← / ← dps glutamine transporter subunit/Fe‑binding and storage protein; stress‑inducible DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913848,3620AG 6.8% 301.3 / 9.5 103intergenic (‑358/+46)glnH/dpsglutamine transporter subunit/Fe‑binding and storage protein; stress‑inducible DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (55/41);  new base G (5/2);  total (60/43)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01

Only 100 of 104 total aligned reads displayed.

GTAAAACCTGATGAAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTTTCGCTTCGCC  >  NC_000913/848213‑848503
                                                                                                                                                     |                                                                                                                                             
gTAAAACCTGATGAAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGaa                                                                                                                                               <  1:1251331/150‑1 (MQ=255)
 tAAAACCTGATGAAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAg                                                                                                                                              >  2:1788823/1‑150 (MQ=255)
   aaaCCTGATGAAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGcc                                                                                                                                            <  1:690446/150‑1 (MQ=255)
     aCCTGATGAAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGc                                                                                                                                          >  2:577357/1‑150 (MQ=255)
     aCCTGATGAAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGc                                                                                                                                          >  2:1424897/1‑150 (MQ=255)
      ccTGATGAAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCt                                                                                                                                         <  2:502414/150‑1 (MQ=255)
        tgatgaAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCttt                                                                                                                                       <  2:1824728/150‑1 (MQ=255)
          attaaaaaaCTATCTTTATGTAACGATTGTGTGATGATGTGTATACAAAAAATATACTTAATCAGCTGCAGCTAGCAGATGGATGCACTAAATACGTGCGTTGAGGTGGGCTGAAAAACAAAACGGCTTCCTGACAGGAGGCCGttttgt                                                                                                                                     <  2:103095/147‑3 (MQ=255)
               aaaaCAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCggg                                                                                                                                <  1:226355/150‑1 (MQ=255)
               aaaaCAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCggg                                                                                                                                <  2:554129/150‑1 (MQ=255)
                aaaCAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGt                                                                                                                               >  2:1811751/1‑150 (MQ=255)
                 aaCAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCATCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTa                                                                                                                              <  2:288699/150‑1 (MQ=255)
                     aTCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTaa                                                                                                                          >  2:40146/1‑150 (MQ=255)
                      tCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTaaa                                                                                                                         <  2:46480/150‑1 (MQ=255)
                       cTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAg                                                                                                                        >  2:1613755/1‑150 (MQ=255)
                       cTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAg                                                                                                                        <  2:1088842/150‑1 (MQ=255)
                           aTGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCt                                                                                                                    >  2:1572009/1‑150 (MQ=255)
                            tGTAACGATTGTGTGATGATGTGGATACTAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTg                                                                                                                   >  1:233889/1‑150 (MQ=255)
                             gTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGc                                                                                                                  <  2:655001/150‑1 (MQ=255)
                             gTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGc                                                                                                                  <  1:1333053/150‑1 (MQ=255)
                              tAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGTGTACTAAAGTTCTGCa                                                                                                                 >  1:1645708/1‑150 (MQ=255)
                              tAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCa                                                                                                                 >  1:924502/1‑150 (MQ=255)
                              tAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGAGGGCTGCAAAACAAAACGGCCTCCTGACAGGAAGCCGCTTTTATCGGGTACTAAACTTCTACa                                                                                                                 >  1:337505/1‑150 (MQ=255)
                                aCGATTTTGTGATGATGTGGATAAAAAAAATTTACTTAATCAGCTTGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCAcc                                                                                                               <  2:1197459/150‑1 (MQ=255)
                                  gATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCAt                                                                                                             >  1:17255/1‑150 (MQ=255)
                                     tgtgtgATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATAGGGTAATAAAGATCTGCACCATCac                                                                                                          >  2:104683/1‑149 (MQ=255)
                                      gtgtgATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGATGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGc                                                                                                         <  1:219174/150‑1 (MQ=255)
                                      gtgtgATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGc                                                                                                         >  2:1576394/1‑150 (MQ=255)
                                       tgtgATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCg                                                                                                        <  2:817094/150‑1 (MQ=255)
                                       tgtgATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCg                                                                                                        <  1:1330218/150‑1 (MQ=255)
                                        gtgATGATGTGGATACAAAAAATTTACTTAATCAGCTGTAGATAGCAGATGGATGCACTAATTAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGATCTGAACAATCAGCGa                                                                                                       >  1:136082/1‑150 (MQ=255)
                                         tgatgaTGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGAt                                                                                                      >  1:213647/1‑150 (MQ=255)
                                         tgatgaTGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGAt                                                                                                      <  1:697434/150‑1 (MQ=255)
                                         tgatgaTGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCAGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGAt                                                                                                      <  2:106524/150‑1 (MQ=255)
                                          gatgatGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATg                                                                                                     <  1:1656801/150‑1 (MQ=255)
                                          gatgatGTGGATAAAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATg                                                                                                     <  2:362899/150‑1 (MQ=255)
                                              atgTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAACACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGAtt                                                                                                 >  2:710687/1‑150 (MQ=255)
                                              atgTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGAtt                                                                                                 <  2:105233/150‑1 (MQ=255)
                                               tgtgGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGAttt                                                                                                >  2:530027/1‑150 (MQ=255)
                                                   gATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTc                                                                                            <  1:1789186/150‑1 (MQ=255)
                                                   gATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTc                                                                                            <  1:732009/150‑1 (MQ=255)
                                                     tACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGTGCTGCAAAACAAACCGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGa                                                                                          <  2:288493/150‑1 (MQ=255)
                                                     tACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGa                                                                                          <  2:509127/150‑1 (MQ=255)
                                                     tACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGa                                                                                          >  1:1140044/1‑150 (MQ=255)
                                                      aCAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGAt                                                                                         <  1:1238408/150‑1 (MQ=255)
                                                       cAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGCCCTCCTGTCAGGAAGCCGCTTATATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGagg                                                                                        >  2:1675614/1‑148 (MQ=255)
                                                       cAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGAAGCACTAAATAAGTGCGTTGAGGAGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATg                                                                                        >  1:1471512/1‑150 (MQ=255)
                                                        aaaaaaTTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGt                                                                                       >  2:1300790/1‑150 (MQ=255)
                                                          aaaaTTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTa                                                                                     <  2:1258118/150‑1 (MQ=255)
                                                                tACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCg                                                                               >  2:1316761/1‑150 (MQ=255)
                                                                 aCTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGa                                                                              >  2:461479/1‑150 (MQ=255)
                                                                  cTTAATCAGCTGGAGCTAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGAt                                                                             <  2:965867/150‑1 (MQ=255)
                                                                    tAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATaa                                                                           >  1:1188662/1‑150 (MQ=255)
                                                                    tAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATaa                                                                           >  2:387153/1‑150 (MQ=255)
                                                                    tAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATaa                                                                           >  1:1040007/1‑150 (MQ=255)
                                                                    tAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATaa                                                                           >  1:1413834/1‑150 (MQ=255)
                                                                     aaTCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATaaa                                                                          >  2:1097062/1‑150 (MQ=255)
                                                                       tCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAAcc                                                                        >  1:389285/1‑150 (MQ=255)
                                                                           cTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTTGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAg                                                                    >  2:1765110/1‑150 (MQ=255)
                                                                               agaTAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAAt                                                                >  1:1415307/1‑150 (MQ=255)
                                                                                   aGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTAt                                                            <  2:616153/150‑1 (MQ=255)
                                                                                   aGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCAACAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTaa                                                            >  2:640312/1‑149 (MQ=255)
                                                                                      aGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCa                                                         <  1:805631/150‑1 (MQ=255)
                                                                                           gatgCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTcg                                                    >  2:1716633/1‑150 (MQ=255)
                                                                                              gCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCga                                                 >  2:475038/1‑150 (MQ=255)
                                                                                               cACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCgag                                                >  2:273956/1‑150 (MQ=255)
                                                                                               cACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCgag                                                <  1:583444/150‑1 (MQ=255)
                                                                                                   aaaTAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGAcgc                                            >  1:706730/1‑150 (MQ=255)
                                                                                                             gTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGTTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAgg                                  <  2:1206198/150‑1 (MQ=255)
                                                                                                             gTTGAGGGGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTTGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGGCGCGAGACGCGGCGGTCAgg                                  <  1:1804513/150‑1 (MQ=255)
                                                                                                                  ggTCGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGTGTACTAAAGCTCTGCACCATCAGCGAAGGATTTATTCGATGCTAGACTCGATAAACCACAGGAATTTATCCAGGCCGCGAGACGCGGCGGTCAGGATATc                             <  1:73134/150‑1 (MQ=255)
                                                                                                                    tGGGCTGCAACACAATACGGCCTCCTGTCAGGAAGCCGCTTTTATCTGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCTCAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTg                           <  1:104683/150‑1 (MQ=255)
                                                                                                                    tGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATAACTg                           >  2:666069/1‑150 (MQ=255)
                                                                                                                     gggCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGc                          >  1:1642551/1‑150 (MQ=255)
                                                                                                                       gCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCgg                        >  1:761705/1‑150 (MQ=255)
                                                                                                                          gCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGgtgt                     >  1:690701/1‑150 (MQ=255)
                                                                                                                          gCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGgtgt                     >  2:1267365/1‑150 (MQ=255)
                                                                                                                          gCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGgtgt                     >  2:1077430/1‑150 (MQ=255)
                                                                                                                            cgataaaagcGGCTTCCTGACAGGAGGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTgtcg                   >  2:1115326/6‑150 (MQ=255)
                                                                                                                             aaacaaaacGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTgtcgt                  >  1:582388/1‑150 (MQ=255)
                                                                                                                              aacaaaacGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTgtcgtc                 >  1:1276889/1‑150 (MQ=255)
                                                                                                                              aacaaaacGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTgtcgtc                 <  1:1576394/150‑1 (MQ=255)
                                                                                                                              aacaaaacGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTgtcgtc                 >  2:1734616/1‑150 (MQ=255)
                                                                                                                              aacaaaacGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTAGCGAGACGCGGCGGTCAGGATATCTGCGGTgtcgtc                 >  1:1798908/1‑150 (MQ=255)
                                                                                                                              aacaaaacGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATAAGCGATGGATTTATTCGATGTTAGACTCGATAATCCACAGGAATTTATACAGGTCGCGAGACGCGGCGGTAAGGAGATCTGCGgtggcctc                 >  2:135977/1‑145 (MQ=255)
                                                                                                                               acaaaacGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCa                <  1:1097062/150‑1 (MQ=255)
                                                                                                                                caaaacGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCAt               >  2:376928/1‑150 (MQ=255)
                                                                                                                                 aaaacGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATc              >  2:378172/1‑150 (MQ=255)
                                                                                                                                 aaaacGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCAACATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATc              >  2:387278/1‑150 (MQ=255)
                                                                                                                                  aagcGGCTTCCTGACAGGAGGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCt             >  1:1039385/4‑150 (MQ=255)
                                                                                                                                  aagcGGCTTCCTGACAGGAGGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCt             >  1:183441/4‑150 (MQ=255)
                                                                                                                                  aagcGGCTTCCTGACAGGAGGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCt             >  1:988599/4‑150 (MQ=255)
                                                                                                                                   agcGGCTTCCTGACAGGAGGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCtt            >  2:485048/3‑150 (MQ=255)
                                                                                                                                     cGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTttc          >  1:164707/1‑150 (MQ=255)
                                                                                                                                     cGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTttc          <  2:1185510/150‑1 (MQ=255)
                                                                                                                                       gCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTttcgc        <  1:667482/150‑1 (MQ=255)
                                                                                                                                          tcctGTCAGGAAGCCGCTTTTATCGGGTACTAATGTTCTGCACCTTCATCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCTTCTttcgcttc     <  2:1413088/150‑1 (MQ=255)
                                                                                                                                           cctGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGATTTTATCCAGGTCGCGAGCCGCGGCGGTCGGGGGATCTGCGGGGGCGGCATCTttcgcgtgg    >  2:1674935/1‑146 (MQ=255)
                                                                                                                                            ctGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTttcgcttcgc   <  1:1300790/150‑1 (MQ=255)
                                                                                                                                             tGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTTTCGCTTCGcc  >  2:1458667/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                             
GTAAAACCTGATGAAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTTTCGCTTCGCC  >  NC_000913/848213‑848503

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: