Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 848,356 T→A 7.1% intergenic (‑352/+52) 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,3560TA 7.1% 525.1 / 24.9 155intergenic (‑352/+52)glnH/dpsglutamine transporter subunit/Fe‑binding and storage protein; stress‑inducible DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (79/65);  new base A (7/4);  total (86/69)
Fisher's exact test for biased strand distribution p-value = 7.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.96e-02

Only 100 of 156 total aligned reads displayed.

TAAAACCTGATGAAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTTTCGCTTCGCCA  >  NC_000913/848214‑848504
                                                                                                                                              |                                                                                                                                                    
tAAAACCTGATGAAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGaa                                                                                                                                                >  1:589233/1‑149 (MQ=255)
tAAAACCTGATGAAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGaa                                                                                                                                                <  3:970909/149‑1 (MQ=255)
     ccTGATGAAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGc                                                                                                                                           <  4:346095/149‑1 (MQ=255)
      cTGATGAAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCt                                                                                                                                          <  1:369549/149‑1 (MQ=255)
      cTGATGAAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCt                                                                                                                                          >  2:1025025/1‑149 (MQ=255)
       tgatgaAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCTTCCTGACAGGAGGCCGttt                                                                                                                                         <  2:975371/149‑1 (MQ=255)
       tgatgaAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCtt                                                                                                                                         >  2:206557/1‑149 (MQ=255)
              aaaaCAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGc                                                                                                                                           >  1:353355/1‑140 (MQ=255)
              aaaaCAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGc                                                                                                                                           <  2:353355/140‑1 (MQ=255)
              aaaaCAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCgg                                                                                                                                  >  3:224365/1‑149 (MQ=255)
                aaCAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGt                                                                                                                                >  1:100312/1‑149 (MQ=255)
                 aCAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTa                                                                                                                               >  3:456868/1‑149 (MQ=255)
                   aaTCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACt                                                                                                                             >  2:832370/1‑149 (MQ=255)
                      cTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTaaa                                                                                                                          <  2:15920/149‑1 (MQ=255)
                      cTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTaaa                                                                                                                          <  4:122766/149‑1 (MQ=255)
                           tGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCg                                                                                                                                   >  1:231903/1‑135 (MQ=255)
                           tGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCt                                                                                                                     <  4:559688/149‑1 (MQ=255)
                            gTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTg                                                                                                                    >  2:120773/1‑149 (MQ=255)
                            gTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTg                                                                                                                    <  2:831220/149‑1 (MQ=255)
                             tAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGAGGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGc                                                                                                                   >  2:775209/1‑149 (MQ=255)
                              aaCGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCa                                                                                                                  <  2:23375/149‑1 (MQ=255)
                                        tgatgaTGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGa                                                                                                        >  2:563535/1‑149 (MQ=255)
                                         gatgatGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGAt                                                                                                       <  2:96315/149‑1 (MQ=255)
                                         gatgatGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGAt                                                                                                       <  3:627939/149‑1 (MQ=255)
                                                 ggATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTAt                                                                                               <  2:28081/149‑1 (MQ=255)
                                                     aCAAAAAATTTACTTAATCAGCTGGAGATAGTAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTTAAGTTCTGCACCATCAGCGATGGATTTCTTCGa                                                                                           <  1:315177/149‑1 (MQ=255)
                                                     aCAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGa                                                                                           <  3:197219/149‑1 (MQ=255)
                                                      cAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGAt                                                                                          <  1:717271/149‑1 (MQ=255)
                                                       aaaaaaTTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATg                                                                                         >  1:830583/1‑149 (MQ=255)
                                                       aaaaaaTTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATg                                                                                         <  4:456868/149‑1 (MQ=255)
                                                       aaaaaaTTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATg                                                                                         >  3:948742/1‑149 (MQ=255)
                                                       aaaaaaTTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATg                                                                                         >  4:301457/1‑149 (MQ=255)
                                                        aaaaaTTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGt                                                                                        >  2:1020388/1‑149 (MQ=255)
                                                        aaaaaTTTACTTAACCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGt                                                                                        >  3:83085/1‑149 (MQ=255)
                                                                    aaTCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGAAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGTTGGATTTATTCGATGTTAGACTCGATaa                                                                            >  4:213443/1‑149 (MQ=255)
                                                                     aTCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATaaa                                                                           >  2:719949/1‑149 (MQ=255)
                                                                     aTCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATaaa                                                                           <  1:1020388/149‑1 (MQ=255)
                                                                      tCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAAc                                                                          <  2:910331/149‑1 (MQ=255)
                                                                      tCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAAc                                                                          <  2:348919/149‑1 (MQ=255)
                                                                          cTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGGCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACcaca                                                                      <  2:600787/149‑1 (MQ=255)
                                                                          cTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACcaca                                                                      <  3:230505/149‑1 (MQ=255)
                                                                           tGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAg                                                                     >  4:232885/1‑149 (MQ=255)
                                                                              agaTAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGaa                                                                  <  2:968858/149‑1 (MQ=255)
                                                                               gaTAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAAt                                                                 <  4:450269/149‑1 (MQ=255)
                                                                               gaTAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAAt                                                                 >  1:85992/1‑149 (MQ=255)
                                                                               gaTAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAAt                                                                 <  1:1083401/149‑1 (MQ=255)
                                                                                  aGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATg                                                                                         <  2:456895/122‑1 (MQ=255)
                                                                                  aGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATg                                                                                         >  1:456895/1‑122 (MQ=255)
                                                                                     aGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATcc                                                           <  3:341130/149‑1 (MQ=255)
                                                                                      gatggatgCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATcc                                                           <  1:49202/148‑1 (MQ=255)
                                                                                          gatgCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTc                                                      <  3:690002/149‑1 (MQ=255)
                                                                                           atgCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTcg                                                     >  4:1045307/1‑149 (MQ=255)
                                                                                               aCTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCgag                                                 >  1:168016/1‑149 (MQ=255)
                                                                                               aCTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCgag                                                 >  3:466350/1‑149 (MQ=255)
                                                                                               aCTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCgag                                                 >  2:1098015/1‑149 (MQ=255)
                                                                                               aCTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCgag                                                 >  1:123043/1‑149 (MQ=255)
                                                                                                cTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCgaga                                                >  1:75519/1‑149 (MQ=255)
                                                                                                 tAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATTAACCACAGGAATTTATCCa                                                          <  3:635351/138‑1 (MQ=255)
                                                                                                 tAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGAc                                               >  4:72460/1‑149 (MQ=255)
                                                                                                 tAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGAc                                               >  2:274613/1‑149 (MQ=255)
                                                                                                 tAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCCCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGaa                                                                  >  1:1085598/1‑130 (MQ=255)
                                                                                                   aaTAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGAcgc                                             <  3:821437/149‑1 (MQ=255)
                                                                                                   aaTAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGAcgc                                             >  1:562278/1‑149 (MQ=255)
                                                                                                   aaTAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTCTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGAcgc                                             <  2:546531/149‑1 (MQ=255)
                                                                                                          gCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGAt                                                                                                   >  1:455180/1‑88 (MQ=255)
                                                                                                            gTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAg                                    >  2:225019/1‑149 (MQ=255)
                                                                                                            gTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCAAGGTCGCGAGACGCGGCGGTCAg                                    <  4:67175/149‑1 (MQ=255)
                                                                                                               gAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGat                                 >  3:314988/1‑149 (MQ=255)
                                                                                                                 ggTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCa                                                          >  4:570551/1‑122 (MQ=255)
                                                                                                                   tGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCt                             <  1:179386/149‑1 (MQ=255)
                                                                                                                     ggCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGc                           >  3:111220/1‑149 (MQ=255)
                                                                                                                       cTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGGATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCgg                         >  1:564102/1‑149 (MQ=255)
                                                                                                                        tGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGgt                        >  3:106213/1‑149 (MQ=255)
                                                                                                                         gCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGgtg                       >  4:329375/1‑149 (MQ=255)
                                                                                                                            aaacaaaacGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTgtcg                    <  1:372948/149‑1 (MQ=255)
                                                                                                                            aaacaaaacGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTgtcg                    >  1:471734/1‑149 (MQ=255)
                                                                                                                             aacaaaacGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTgtcgt                   >  4:291625/1‑149 (MQ=255)
                                                                                                                             aacaaaacGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTgtcgt                   >  3:164167/1‑149 (MQ=255)
                                                                                                                              acaaaacGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTgtcgtc                  <  3:889723/149‑1 (MQ=255)
                                                                                                                              acaaaacGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTgtcgtc                  >  1:973081/1‑149 (MQ=255)
                                                                                                                              acaaaacGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTgtcgtc                  >  3:886577/1‑149 (MQ=255)
                                                                                                                               caaaacGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCa                 <  4:497503/149‑1 (MQ=255)
                                                                                                                               caaaacGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCa                 >  1:887209/1‑149 (MQ=255)
                                                                                                                               caaaacGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCa                 >  3:1010007/1‑149 (MQ=255)
                                                                                                                                 aagcGGCTTCCTGACAGGAGGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGa                                                                                           <  1:976492/70‑1 (MQ=37)
                                                                                                                                 aagcGGCTTCCTGACAGGAGGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAg                                    >  2:529598/4‑128 (MQ=255)
                                                                                                                                 aagcGGCTTCCTGACAGGAGGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATc               >  2:746584/4‑149 (MQ=255)
                                                                                                                                 aaacGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTCCTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATATGCGGTGTCGTCATc               >  4:958888/1‑149 (MQ=255)
                                                                                                                                  aacGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCt              <  4:948742/149‑1 (MQ=255)
                                                                                                                                    cGGCTTCCTGACAGGAGGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCttt            >  3:348241/1‑149 (MQ=255)
                                                                                                                                        ttcctGACAGGAGGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTttcgct        >  1:608777/2‑149 (MQ=255)
                                                                                                                                        ctcctGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAgg                                                        >  4:305425/1‑101 (MQ=255)
                                                                                                                                        ctcctGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTttcgct        >  3:297747/1‑149 (MQ=255)
                                                                                                                                        ctcctGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTttcgct        >  2:249177/1‑149 (MQ=255)
                                                                                                                                        ctcctGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTttcgct        >  1:11992/1‑149 (MQ=255)
                                                                                                                                         tcctGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCCCCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTc                                                      <  2:432436/102‑1 (MQ=255)
                                                                                                                                         tcctGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAggcc                                                      >  1:432436/1‑100 (MQ=255)
                                                                                                                                         tcctGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTttcgctt       >  4:66140/1‑149 (MQ=255)
                                                                                                                                          cctGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTttcgcttc      <  3:291625/149‑1 (MQ=255)
                                                                                                                                          cctGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTttcgcttc      <  2:471734/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
TAAAACCTGATGAAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTTTCGCTTCGCCA  >  NC_000913/848214‑848504

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: