Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 848,350 C→T 8.8% 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 8.8% 276.7 / 13.9 91intergenic (‑346/+58)glnH/dpsglutamine transporter subunit/Fe‑binding and storage protein; stress‑inducible DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (41/42);  new base T (5/3);  total (46/45)
Fisher's exact test for biased strand distribution p-value = 7.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

ATTATTTAGCGTAAAACCTGATGAAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTTT  >  NC_000913/848203‑848494
                                                                                                                                                   |                                                                                                                                                
attattTAGCGTAAAACCTGATGAAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGcct                                                                                                                                                 >  1:627196/1‑149 (MQ=255)
    tttAGCGTAAAACCTGATGAAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTg                                                                                                                                             >  4:633967/1‑149 (MQ=255)
     ttAGCGTAAAACCTGATGAAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGt                                                                                                                                            <  1:678252/149‑1 (MQ=255)
                 cTGATGAAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCt                                                                                                                                <  2:434125/149‑1 (MQ=255)
                   gatgaAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCttt                                                                                                                              <  3:633967/149‑1 (MQ=255)
                    atgaAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCtttt                                                                                                                             >  4:175940/1‑149 (MQ=255)
                      gaAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTAt                                                                                                                           >  2:666332/1‑149 (MQ=255)
                       aaaaaaCAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATc                                                                                                                          >  3:391314/1‑149 (MQ=255)
                           aaCAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGt                                                                                                                      >  1:186190/1‑149 (MQ=255)
                                      tGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCt                                                                                                           >  3:210684/1‑149 (MQ=255)
                                        tAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGc                                                                                                         >  3:687923/1‑149 (MQ=255)
                                         aaCGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTa                                                                                                                     >  1:211206/1‑136 (MQ=255)
                                         aaCGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTa                                                                                                                     <  2:211206/136‑1 (MQ=255)
                                            gATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCa                                                                                                     >  3:37272/1‑149 (MQ=255)
                                                gtgtgATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAg                                                                                                 <  4:97767/149‑1 (MQ=255)
                                                   tgatgaTGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGa                                                                                              >  1:355338/1‑149 (MQ=255)
                                                   tgatgaTGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGa                                                                                              <  4:123017/149‑1 (MQ=255)
                                                    gatgatGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCtt                                                                                                                               >  2:481047/1‑115 (MQ=255)
                                                    gatgatGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCtt                                                                                                                               <  1:481047/115‑1 (MQ=255)
                                                    gatgatGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGAt                                                                                             <  2:611390/149‑1 (MQ=255)
                                                     atgatgTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATg                                                                                            <  2:627196/149‑1 (MQ=255)
                                                       gatgTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTc                                                                                                            <  3:208895/131‑1 (MQ=255)
                                                       gatgTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTc                                                                                                            >  4:208895/1‑131 (MQ=255)
                                                         tgtgGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGAtt                                                                                        <  4:566343/149‑1 (MQ=255)
                                                            ggATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTAt                                                                                     >  2:341663/1‑149 (MQ=255)
                                                             gATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTAtt                                                                                    >  4:508389/1‑149 (MQ=255)
                                                             gATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTAtt                                                                                    >  2:464794/1‑149 (MQ=255)
                                                              aTACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTc                                                                                   >  3:255683/1‑149 (MQ=255)
                                                               tACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCg                                                                                  <  4:210684/149‑1 (MQ=255)
                                                                 cAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATCAGTGCGTTGAGGAGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGAt                                                                                >  4:507906/1‑149 (MQ=255)
                                                                 cAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGAt                                                                                <  3:543311/149‑1 (MQ=255)
                                                                  aaaaaaTTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATg                                                                               >  4:423217/1‑149 (MQ=255)
                                                                   aaaaaTTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGt                                                                              >  1:565973/1‑149 (MQ=255)
                                                                   aaaaaTTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGt                                                                              >  1:457006/1‑149 (MQ=255)
                                                                     aaaTTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCTTCCTGACAGGAGGCCGttt                                                                                                                               <  4:524209/98‑1 (MQ=39)
                                                                     aaaTTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCTTCCTGACAGGAGGCCGttt                                                                                                                               >  3:524209/1‑98 (MQ=38)
                                                                     aaaTTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTa                                                                            >  4:260416/1‑149 (MQ=255)
                                                                       aTTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGa                                                                          <  4:37272/149‑1 (MQ=255)
                                                                       aTTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGa                                                                          <  3:175940/149‑1 (MQ=255)
                                                                          tACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTc                                                                       <  4:564731/149‑1 (MQ=255)
                                                                            cTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGa                                                                     <  2:520339/149‑1 (MQ=255)
                                                                             ttAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGAt                                                                    >  3:171015/1‑149 (MQ=255)
                                                                                aTCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATaaa                                                                 >  1:499033/1‑149 (MQ=255)
                                                                                 tCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAAc                                                                <  1:158461/149‑1 (MQ=255)
                                                                                 tCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCAGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAAc                                                                >  3:645091/1‑149 (MQ=255)
                                                                                     cTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACcaca                                                            <  4:255683/149‑1 (MQ=255)
                                                                                         agaTAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGaa                                                        <  1:359375/149‑1 (MQ=255)
                                                                                          gaTAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAAt                                                       <  3:508389/149‑1 (MQ=255)
                                                                                          gaTAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAAt                                                       <  1:464794/149‑1 (MQ=255)
                                                                                              gCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGaa                                                        <  4:187781/144‑1 (MQ=255)
                                                                                              gCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGaa                                                        >  3:187781/1‑144 (MQ=255)
                                                                                               cAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATc                                                  <  2:355338/149‑1 (MQ=255)
                                                                                                aGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGAt                                                                                             >  4:64922/1‑105 (MQ=255)
                                                                                                aGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGAt                                                                                             <  3:64922/105‑1 (MQ=255)
                                                                                                aGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATcc                                                 <  4:534623/149‑1 (MQ=255)
                                                                                                    ggatgCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGt                                             >  3:600011/1‑149 (MQ=255)
                                                                                                     gatgCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTc                                            <  4:391314/149‑1 (MQ=255)
                                                                                                       tgCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTcgc                                          <  1:341663/149‑1 (MQ=255)
                                                                                                         cACTAAATAAGTGCGTTGAGGTTGGCTGCAAAACAAAACAGCCACCTGTCAGGAAGCCGCTTTTATCGGGTACTAAACTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAAATTATCCAGGTCGCga                                        <  4:645091/149‑1 (MQ=255)
                                                                                                          aCTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTAt                                                                                     <  1:348514/103‑1 (MQ=255)
                                                                                                          aCTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTAt                                                                                     >  2:348514/1‑103 (MQ=255)
                                                                                                          aCTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACcaca                                                            >  4:350968/1‑128 (MQ=255)
                                                                                                          aCTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACcaca                                                            <  3:350968/128‑1 (MQ=255)
                                                                                                          aCTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTAt                                                   <  1:682042/137‑1 (MQ=255)
                                                                                                          aCTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTAt                                                   >  2:682042/1‑137 (MQ=255)
                                                                                                          aCTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCgag                                       >  1:503419/1‑149 (MQ=255)
                                                                                                              aaTAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGAcgc                                   >  4:226617/1‑149 (MQ=255)
                                                                                                                  aGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACgcggcg                               <  1:666332/149‑1 (MQ=255)
                                                                                                                   gTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGcggcgg                              >  1:139881/1‑149 (MQ=255)
                                                                                                                    tGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGa                                                                                          >  1:57883/1‑88 (MQ=255)
                                                                                                                    tGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGa                                                                                          <  2:57883/88‑1 (MQ=255)
                                                                                                                       gTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATcc                                                 >  4:271292/1‑126 (MQ=255)
                                                                                                                       gTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATcc                                                 <  3:271292/126‑1 (MQ=255)
                                                                                                                       gTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAg                          >  4:542039/1‑149 (MQ=255)
                                                                                                                       gTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAACGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAg                          <  2:457006/149‑1 (MQ=255)
                                                                                                                                   tGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAAc                                                                >  1:252111/1‑99 (MQ=255)
                                                                                                                                   tGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAAc                                                                <  2:252111/99‑1 (MQ=255)
                                                                                                                                   tGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGgt              >  1:9514/1‑149 (MQ=255)
                                                                                                                                      aaaacaaaacGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACt                                                                        >  3:398556/1‑88 (MQ=255)
                                                                                                                                      aaaacaaaacGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTgtc           <  1:34239/149‑1 (MQ=255)
                                                                                                                                      aaaacaaaacGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTgtc           <  3:260416/149‑1 (MQ=255)
                                                                                                                                      aaaacaaaacGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTgtc           <  3:584160/149‑1 (MQ=255)
                                                                                                                                      aaaacaaaacGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTgtc           <  2:186190/149‑1 (MQ=255)
                                                                                                                                      aaaacaaaacGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTGCTGCACCATCAGCGATGGATTTATTCGATGTTAGACt                                                                        <  4:398556/88‑1 (MQ=255)
                                                                                                                                         acaaaacGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTgtcgtc        >  3:661287/1‑149 (MQ=255)
                                                                                                                                            aagcGGCTTCCTGACAGGAGGCCGTTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATc     >  3:630983/4‑149 (MQ=255)
                                                                                                                                              gcGGCTTCCTGACAGGAGGCCGTTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCtt   >  2:36756/2‑149 (MQ=255)
                                                                                                                                              gcGGCTTCCTGACAGGAGGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACcaca                                                            <  2:590130/91‑1 (MQ=255)
                                                                                                                                              gcGGCTTCCTGACAGGAGGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACcaca                                                            >  1:590130/2‑92 (MQ=255)
                                                                                                                                               cGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCttt  >  1:185970/1‑149 (MQ=255)
                                                                                                                                                ggCTTCCTGACAGGAGGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATc                                                  <  4:50752/100‑1 (MQ=255)
                                                                                                                                                ggCTTCCTGACAGGAGGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATc                                                  >  3:50752/1‑100 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
ATTATTTAGCGTAAAACCTGATGAAAAAACAATCTTTATGTAACGATTGTGTGATGATGTGGATACAAAAAATTTACTTAATCAGCTGGAGATAGCAGATGGATGCACTAAATAAGTGCGTTGAGGTGGGCTGCAAAACAAAACGGCCTCCTGTCAGGAAGCCGCTTTTATCGGGTACTAAAGTTCTGCACCATCAGCGATGGATTTATTCGATGTTAGACTCGATAAACCACAGGAATTTATCCAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTTT  >  NC_000913/848203‑848494

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