Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 848,350 C→T 6.9% intergenic (‑346/+58) 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,3500CT 6.9% 462.5 / 17.5 144intergenic (‑346/+58)glnH/dpsglutamine transporter subunit/Fe‑binding and storage protein; stress‑inducible DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (73/61);  new base T (6/4);  total (79/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.33e-01

Only 100 of 152 total aligned reads displayed.

GATTATTTAGCGTAAAACCTGATGAAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTTTCGCT  >  NC_000913/848202‑848498
                                                                                                                                                    |                                                                                                                                                    
gATTATTTAGCGTAAAACCTGATGAAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTAGTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTTCGTTGAGGTGGGCTGCAAAACAAAACGGcc                                                                                                                                                      >  3:47084/1‑149 (MQ=255)
     tttAGCGTAAAACCTGATGAAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTg                                                                                                                                                 >  1:553797/1‑149 (MQ=255)
            tAAAACCTGATGAAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGaa                                                                                                                                          >  1:475206/1‑149 (MQ=255)
            tAAAACCTGATGAAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGaa                                                                                                                                          <  3:970909/149‑1 (MQ=255)
                 ccTGATGAAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGc                                                                                                                                     <  4:346095/149‑1 (MQ=255)
                  cTGATGAAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGcctcct                                                                                                                                                  >  1:15786/1‑135 (MQ=255)
                  cTGATGAAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGcctcct                                                                                                                                                  <  2:15786/135‑1 (MQ=255)
                  cTGATGAAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCt                                                                                                                                    >  2:1025025/1‑149 (MQ=255)
                  cTGATGAAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCt                                                                                                                                    <  1:369549/149‑1 (MQ=255)
                   tgatgaAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCTTCCTGACAGGAGGCCGttt                                                                                                                                   <  2:975371/149‑1 (MQ=255)
                   tgatgaAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCtt                                                                                                                                   >  2:206557/1‑149 (MQ=255)
                   tgatgaAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCtt                                                                                                                                   <  3:426274/149‑1 (MQ=255)
                       gaAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTAt                                                                                                                               >  1:679705/1‑149 (MQ=255)
                          aaaaCAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGc                                                                                                                                     >  1:353355/1‑140 (MQ=255)
                          aaaaCAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGc                                                                                                                                     <  2:353355/140‑1 (MQ=255)
                            aaCAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGt                                                                                                                          >  1:100312/1‑149 (MQ=255)
                               aaTCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACt                                                                                                                       >  2:832370/1‑149 (MQ=255)
                                aTCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTa                                                                                                                      >  4:548720/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:831220/149‑1 (MQ=255)
                                        gTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTg                                                                                                              >  2:120773/1‑149 (MQ=255)
                                         tAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGAGGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGc                                                                                                             >  2:775209/1‑149 (MQ=255)
                                          aaCGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCa                                                                                                            <  2:23375/149‑1 (MQ=255)
                                                    tgatgaTGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGa                                                                                                  >  2:563535/1‑149 (MQ=255)
                                                     gatgatGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGTGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGAt                                                                                                 <  4:307924/149‑1 (MQ=255)
                                                     gatgatGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGAt                                                                                                 >  3:217422/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)
                                                                  cAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGAt                                                                                    <  1:206557/149‑1 (MQ=255)
                                                                   aaaaaaTTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATg                                                                                   >  3:948742/1‑149 (MQ=255)
                                                                   aaaaaaTTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATg                                                                                   >  1:830583/1‑149 (MQ=255)
                                                                   aaaaaaTTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATg                                                                                   <  4:456868/149‑1 (MQ=255)
                                                                    aaaaaTTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGt                                                                                  >  2:1020388/1‑149 (MQ=255)
                                                                    aaaaaTTTACTTAACCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGt                                                                                  >  3:83085/1‑149 (MQ=255)
                                                                             cTTAATCACCTCGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAGAACGGACTCCTGTCAGGAAGCCGCTATTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTCTTCGATGTTAGACTCGa                                                                         <  4:436729/149‑1 (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:348919/149‑1 (MQ=255)
                                                                                  tCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAAc                                                                    <  2:910331/149‑1 (MQ=255)
                                                                                      cTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGGCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACcaca                                                                <  2:600787/149‑1 (MQ=255)
                                                                                      cTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACcaca                                                                <  3:230505/149‑1 (MQ=255)
                                                                                      cTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACcaca                                                                <  4:217422/149‑1 (MQ=255)
                                                                                      cTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACcaca                                                                <  4:297199/149‑1 (MQ=255)
                                                                                          agaTAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGaa                                                            <  2:968858/149‑1 (MQ=255)
                                                                                           gaTAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAAt                                                           >  1:85992/1‑149 (MQ=255)
                                                                                           gaTAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAAt                                                           <  1:1083401/149‑1 (MQ=255)
                                                                                           gaTAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAAt                                                           <  3:1016429/149‑1 (MQ=255)
                                                                                              aGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATg                                                                                   <  2:456895/122‑1 (MQ=255)
                                                                                                      gatgCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTc                                                <  3:690002/149‑1 (MQ=255)
                                                                                                       atgCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTcg                                               >  4:1045307/1‑149 (MQ=255)
                                                                                                           aCTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCgag                                           >  2:1098015/1‑149 (MQ=255)
                                                                                                           aCTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCgag                                           >  1:123043/1‑149 (MQ=255)
                                                                                                           aCTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCgag                                           >  1:168016/1‑149 (MQ=255)
                                                                                                           aCTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCgag                                           >  3:466350/1‑149 (MQ=255)
                                                                                                           aCTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCgag                                           >  1:456418/1‑149 (MQ=255)
                                                                                                            cTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCgaga                                          <  4:466350/149‑1 (MQ=255)
                                                                                                             tAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATTAACCACAGGAATTTATCCa                                                    <  3:635351/138‑1 (MQ=255)
                                                                                                             tAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCa                                                    >  4:635351/1‑138 (MQ=255)
                                                                                                             tAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGAc                                         >  2:274613/1‑149 (MQ=255)
                                                                                                             tAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGAc                                         >  4:72460/1‑149 (MQ=255)
                                                                                                             tAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCCCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGaa                                                            >  1:1085598/1‑130 (MQ=255)
                                                                                                              aaaTAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGAcg                                        >  3:239404/1‑149 (MQ=255)
                                                                                                               aaTAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGAcgc                                       <  3:821437/149‑1 (MQ=255)
                                                                                                               aaTAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTCTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGAcgc                                       <  2:546531/149‑1 (MQ=255)
                                                                                                                      gCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGAt                                                                                             >  1:455180/1‑88 (MQ=255)
                                                                                                                      gCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGAt                                                                                             <  2:455180/88‑1 (MQ=255)
                                                                                                                        gTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAg                              <  1:586531/149‑1 (MQ=255)
                                                                                                                        gTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAg                              >  2:225019/1‑149 (MQ=255)
                                                                                                                        gTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCAAGGTCGCGAGACGCGGCGGTCAg                              <  4:67175/149‑1 (MQ=255)
                                                                                                                           gAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGat                           >  3:314988/1‑149 (MQ=255)
                                                                                                                                 ggCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGc                     >  3:111220/1‑149 (MQ=255)
                                                                                                                                   cTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGGATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCgg                   >  1:564102/1‑149 (MQ=255)
                                                                                                                                        aaacaaaacGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTgtcg              >  1:471734/1‑149 (MQ=255)
                                                                                                                                        aaacaaaacGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTgtcg              <  1:372948/149‑1 (MQ=255)
                                                                                                                                         aacaaaacGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTgtcgt             >  4:291625/1‑149 (MQ=255)
                                                                                                                                          acaaaacGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTgtcgtc            >  3:886577/1‑149 (MQ=255)
                                                                                                                                          acaaaacGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTgtcgtc            <  3:889723/149‑1 (MQ=255)
                                                                                                                                          acaaaacGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTgtcgtc            >  1:973081/1‑149 (MQ=255)
                                                                                                                                           caaaacGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCa           >  1:887209/1‑149 (MQ=255)
                                                                                                                                           caaaacGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCa           >  3:1010007/1‑149 (MQ=255)
                                                                                                                                             aagcGGCTTCCTGACAGGAGGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGa                                                                                     <  1:976492/70‑1 (MQ=37)
                                                                                                                                             aagcGGCTTCCTGACAGGAGGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAg                              <  1:529598/125‑1 (MQ=255)
                                                                                                                                             aagcGGCTTCCTGACAGGAGGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAg                              >  2:529598/4‑128 (MQ=255)
                                                                                                                                             aagcGGCTTCCTGACAGGAGGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATc         >  2:746584/4‑149 (MQ=255)
                                                                                                                                             aaacGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTCCTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATATGCGGTGTCGTCATc         >  4:958888/1‑149 (MQ=255)
                                                                                                                                              agcGGCTTCCTGACAGGAGGCCGTTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGAt                                                                                    >  3:138504/3‑73 (MQ=25)
                                                                                                                                              aacGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCt        <  4:948742/149‑1 (MQ=255)
                                                                                                                                                cGGCTTCCTGACAGGAGGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCttt      >  3:348241/1‑149 (MQ=255)
                                                                                                                                                    ctcctGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAgg                                                  >  4:305425/1‑101 (MQ=255)
                                                                                                                                                    ctcctGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAgg                                                  <  3:305425/101‑1 (MQ=255)
                                                                                                                                                    ctcctGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTttcgct  >  2:249177/1‑149 (MQ=255)
                                                                                                                                                    ctcctGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTttcgct  >  1:11992/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GATTATTTAGCGTAAAACCTGATGAAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTTTCGCT  >  NC_000913/848202‑848498

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