Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 494,359 T→A 5.6% V95E (GTA→GAA)  ybaB → DNA‑binding protein, putative nucleoid‑associated protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913494,3590TA 5.6% 229.9 / 3.4 71V95E (GTA→GAA) ybaBDNA‑binding protein, putative nucleoid‑associated protein
Reads supporting (aligned to +/- strand):  ref base T (39/28);  new base A (2/2);  total (41/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.19e-01

TCAACGGTGCACACAACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGC  >  NC_000913/494212‑494501
                                                                                                                                                   |                                                                                                                                              
tCAACGGTGCACACAACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTa                                                                                                                                               <  1:135130/149‑1 (MQ=255)
 cAACGGTGCACACAACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTa                                                                                                                                               >  6:45713/1‑148 (MQ=255)
 cAACGGTGCACACAACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTa                                                                                                                                               <  5:45713/148‑1 (MQ=255)
       tGCACACAACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCg                                                                                                                                        >  6:11911/1‑149 (MQ=255)
           cacaACTGCCGTCGCGTCGAGATCGACCCGAGCCTGCTGGAAGACGACAACGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAAt                                                                                                                                    <  2:211289/149‑1 (MQ=255)
              aaCTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCa                                                                                                                                 <  6:215540/149‑1 (MQ=255)
                cTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGc                                                                                                                               >  7:259284/1‑149 (MQ=255)
                cTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGAATCCTCCGGAATGCAGc                                                                                                                               >  1:163679/1‑149 (MQ=255)
                cTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGAATCCTCCGGAATGCAGc                                                                                                                               >  1:169022/1‑149 (MQ=255)
                 tGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACTTCGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCCGCt                                                                                                                              >  5:211978/1‑149 (MQ=255)
                                aTCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCttt                                                                                                                 >  3:118455/1‑147 (MQ=255)
                                        gAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGaa                                                                                                                                     <  3:311090/119‑1 (MQ=255)
                                        gAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGaa                                                                                                                                     >  4:311090/1‑119 (MQ=255)
                                             tgctgGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGa                                                                                                  >  2:182937/1‑149 (MQ=255)
                                              gctgGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGAt                                                                                                 >  3:170025/1‑149 (MQ=255)
                                               ctgGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGAt                                                                                                 <  6:3617/148‑1 (MQ=255)
                                               ctgGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGAATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGAt                                                                                                 <  2:169022/148‑1 (MQ=255)
                                               ctgGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGAATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGAt                                                                                                 <  2:163679/148‑1 (MQ=255)
                                                tgGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATg                                                                                                >  8:353930/1‑148 (MQ=255)
                                                    aGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAAc                                                                                           <  4:118455/149‑1 (MQ=255)
                                                    aGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGGCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAAc                                                                                           <  3:186640/149‑1 (MQ=255)
                                                     gacgacAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAAc                                                                                           >  1:47678/1‑148 (MQ=255)
                                                        gacAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGCTGCAAACCAGc                                                                                       >  4:39795/1‑149 (MQ=255)
                                                          cAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATg                                                                                                           <  5:18631/127‑1 (MQ=255)
                                                          cAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATg                                                                                                           >  6:18631/1‑127 (MQ=255)
                                                           aaaGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTaa                                                                                                               >  4:142665/1‑122 (MQ=255)
                                                           aaaGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTaa                                                                                                               <  3:142665/122‑1 (MQ=255)
                                                           aaaGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCg                                                                                    >  1:160513/1‑149 (MQ=255)
                                                                     tGGAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAaca                                                                           >  3:76076/1‑148 (MQ=255)
                                                                     tGGAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAaca                                                                           >  1:307374/1‑148 (MQ=255)
                                                                     tGGAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAaca                                                                           >  4:324835/1‑148 (MQ=255)
                                                                       gAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAg                                                                        >  2:268217/1‑149 (MQ=255)
                                                                       gAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAg                                                                        >  2:262732/1‑149 (MQ=255)
                                                                            ccTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTAt                                                                   <  2:160513/149‑1 (MQ=255)
                                                                             cTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCt                                                                      >  8:79450/1‑145 (MQ=255)
                                                                             cTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCt                                                                      <  7:79450/145‑1 (MQ=255)
                                                                             cTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATg                                                                  >  6:117763/1‑149 (MQ=255)
                                                                             cTGGTGGCTCCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCAGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTACCACAGCTTATg                                                                  >  2:340736/1‑149 (MQ=255)
                                                                                gtggCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGaa                                                               <  2:47678/149‑1 (MQ=255)
                                                                                gaggCTGCAGCATGCACCGGCGCAGCACGTCGTATTGAAGAAACCCAGAAAGAACAACTCGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGa                                                                <  6:211978/146‑1 (MQ=255)
                                                                                       cagcaTTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACt                                                          <  8:259284/147‑1 (MQ=255)
                                                                                            ttCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGt                                                   >  2:71936/1‑149 (MQ=255)
                                                                                               aaCGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTACCACAGCTTATGGAAGCACTGCGCTGTCTg                                                >  8:83493/1‑149 (MQ=255)
                                                                                               aaCGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTg                                                >  3:132594/1‑149 (MQ=255)
                                                                                               aaCGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTg                                                >  7:72663/1‑149 (MQ=255)
                                                                                               aaCGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTg                                                >  4:327187/1‑149 (MQ=255)
                                                                                                acgacgCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGc                                               >  6:216483/1‑149 (MQ=255)
                                                                                                acgacgCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGc                                               <  4:76076/149‑1 (MQ=255)
                                                                                                   acgCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCgg                                            <  5:216483/149‑1 (MQ=255)
                                                                                                       agcaCGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGt                                        >  7:231939/1‑149 (MQ=255)
                                                                                                         caCGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTg                                      <  1:182937/149‑1 (MQ=255)
                                                                                                         caCGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTg                                      >  5:358580/1‑149 (MQ=255)
                                                                                                             tcgtATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGccc                                  >  8:233774/1‑149 (MQ=255)
                                                                                                                tATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGaa                               <  1:71936/149‑1 (MQ=255)
                                                                                                                    gaagaaACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCg                           >  3:94993/1‑149 (MQ=255)
                                                                                                                    gaagaaACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCg                           >  8:216902/1‑149 (MQ=255)
                                                                                                                     aagaaACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCg                           <  3:324835/148‑1 (MQ=255)
                                                                                                                       gaaACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGgcg                        >  2:221438/1‑149 (MQ=255)
                                                                                                                       gaaACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGgcg                        >  1:299207/1‑149 (MQ=255)
                                                                                                                         aaCGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCa                      <  7:353930/149‑1 (MQ=255)
                                                                                                                          aCGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAg                     <  1:262732/149‑1 (MQ=255)
                                                                                                                          aCGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAg                     <  1:268217/149‑1 (MQ=255)
                                                                                                                            gCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCg                   >  8:130508/1‑149 (MQ=255)
                                                                                                                              agaaagaaAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTa                 >  1:324179/1‑149 (MQ=255)
                                                                                                                                   gaaaAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCgg                          >  7:260754/1‑135 (MQ=255)
                                                                                                                                   gaaaAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCgg                          <  8:260754/135‑1 (MQ=255)
                                                                                                                                     aaaaaTGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGtt          <  7:130508/149‑1 (MQ=255)
                                                                                                                                            gCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACgctg   <  4:132594/149‑1 (MQ=255)
                                                                                                                                             ccTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACgctgc  <  5:11911/149‑1 (MQ=255)
                                                                                                                                             ccTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACgctgc  >  5:1178/1‑149 (MQ=255)
                                                                                                                                             ccTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACgctgc  <  2:307374/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
TCAACGGTGCACACAACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGC  >  NC_000913/494212‑494501

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