Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,444,193 A→C 5.6% T28P (ACC→CCC)  tamB → translocation and assembly module for autotransporter export, inner membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,444,1930AC 5.6% 255.9 / 5.2 72T28P (ACC→CCC) tamBtranslocation and assembly module for autotransporter export, inner membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (39/29);  new base C (2/2);  total (41/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

GATAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCT  >  NC_000913/4444062‑4444341
                                                                                                                                   |                                                                                                                                                    
gATAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCAt                                                                                                                                     >  1:1240351/1‑149 (MQ=255)
     aGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGt                                                                                                                                >  1:1235432/1‑149 (MQ=255)
       aCGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTAt                                                                                                                              >  1:815411/1‑149 (MQ=255)
            cacGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTaaa                                                                                                                         >  2:922631/1‑149 (MQ=255)
                           tACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTgg                                                                                                          >  1:1068550/1‑149 (MQ=255)
                              aTCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGTGCCCACCAgaa                                                                                                                                             >  1:1348278/1‑109 (MQ=255)
                              aTCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGTGCCCACCAgaa                                                                                                                                             <  2:1348278/111‑3 (MQ=255)
                              aTCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAgcgg                                                                                                                     >  4:193446/1‑135 (MQ=255)
                              aTCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAgcgg                                                                                                                     <  3:193446/135‑1 (MQ=255)
                              aTCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGc                                                                                                             <  2:729855/143‑1 (MQ=255)
                              aTCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGc                                                                                                             >  1:729855/1‑143 (MQ=255)
                              aTCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTg                                                                                                       >  2:140219/1‑149 (MQ=255)
                                        ggCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTgg                                                                                             <  2:1240351/149‑1 (MQ=255)
                                           cAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGAc                                                                                                >  2:655304/1‑143 (MQ=255)
                                           cAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGAc                                                                                                <  1:655304/143‑1 (MQ=255)
                                            aGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGatat                                                                                         >  1:911567/1‑149 (MQ=255)
                                             gAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTaaa                                                                                                                         >  4:1412844/1‑116 (MQ=255)
                                             gAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTaaa                                                                                                                         <  3:1412844/116‑1 (MQ=255)
                                              aaTTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGGTATTg                                                                                       <  1:922631/149‑1 (MQ=255)
                                              aaTTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTg                                                                                       <  2:911567/149‑1 (MQ=255)
                                                   tGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAg                                                                                  <  2:815411/149‑1 (MQ=255)
                                                     aGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGcaccac                                                                                                                                                  >  1:636940/1‑83 (MQ=255)
                                                     aGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGcaccac                                                                                                                                                  <  2:636940/83‑1 (MQ=255)
                                                       tttATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTg                                                                                       <  2:229087/140‑1 (MQ=255)
                                                       tttATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTg                                                                                       >  1:229087/1‑140 (MQ=255)
                                                       tttATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCaa                                                                                   >  1:1380539/1‑144 (MQ=255)
                                                       tttATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCaa                                                                                   <  2:1380539/144‑1 (MQ=255)
                                                       tttATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCa                                                                              >  1:33152/1‑149 (MQ=255)
                                                             gAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACcggcg                                                                        >  2:1328206/1‑149 (MQ=255)
                                                               aaaaaaTCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCggcggc                                                                      <  2:1235432/149‑1 (MQ=255)
                                                                aaaaaTCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTa                                                                                                                               >  4:214547/1‑91 (MQ=255)
                                                                aaaaaTCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTa                                                                                                                               <  3:214547/91‑1 (MQ=255)
                                                                 aaaaTCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTg                                                                    <  1:1328206/149‑1 (MQ=255)
                                                                          cTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGAtct                                                           >  1:507326/1‑149 (MQ=255)
                                                                          cTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGAtct                                                           >  3:514863/1‑149 (MQ=255)
                                                                              gCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCAcc                                                       <  2:507326/149‑1 (MQ=255)
                                                                                   gttatcgttatcTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGa                                                                                                                 >  2:1336545/1‑86 (MQ=255)
                                                                                   gttatcgttatcTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGa                                                                                                                 <  1:1336545/86‑1 (MQ=255)
                                                                                       tcgttatcTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCt                                                                                                                                         >  4:634401/1‑58 (MQ=255)
                                                                                       tcgttatcTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCt                                                                                                                                         <  3:634401/58‑1 (MQ=255)
                                                                                         gttatcTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGt                                                   >  1:565269/1‑142 (MQ=255)
                                                                                         gttatcTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGt                                                   <  2:565269/142‑1 (MQ=255)
                                                                                         gttatcTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGAcgt                                            >  1:390177/1‑149 (MQ=255)
                                                                                                tACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGcacca                                                                                                                                                   >  2:1195479/1‑39 (MQ=255)
                                                                                                          gggATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGagca                                <  2:589986/144‑1 (MQ=255)
                                                                                                          gggATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGagca                                >  1:589986/1‑144 (MQ=255)
                                                                                                            gATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGc                         <  1:140219/149‑1 (MQ=255)
                                                                                                             aTCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCg                        >  3:520467/1‑149 (MQ=255)
                                                                                                                 gTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGc                    >  3:1172211/1‑149 (MQ=255)
                                                                                                                  tGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCg                   >  4:533175/1‑149 (MQ=255)
                                                                                                                  tGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCg                   >  3:17710/1‑149 (MQ=255)
                                                                                                                   ggCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGagcagc                              <  2:308674/137‑1 (MQ=255)
                                                                                                                   ggCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGagcagc                              >  1:308674/1‑137 (MQ=255)
                                                                                                                   ggCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCgg                  >  2:1429000/1‑149 (MQ=255)
                                                                                                                    gCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGt                 >  2:1429958/1‑149 (MQ=255)
                                                                                                                      gTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGt                                                                                                        >  1:1080486/1‑60 (MQ=255)
                                                                                                                      gTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGt                                                                                                        <  2:1080486/60‑1 (MQ=255)
                                                                                                                      gTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTaa               >  3:1368054/1‑149 (MQ=255)
                                                                                                                        cgcTGGTGGTGCCCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGAtctc                                                          <  3:359447/102‑1 (MQ=255)
                                                                                                                        cgcTGGTGGTGCCCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGAtctc                                                          >  4:359447/3‑104 (MQ=255)
                                                                                                                           tggtggGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGAt                                                             <  4:993093/98‑1 (MQ=255)
                                                                                                                           tggtggGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGAt                                                             >  3:993093/1‑98 (MQ=255)
                                                                                                                            ggtggGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTAAAAGCgg         <  1:1429000/149‑1 (MQ=255)
                                                                                                                             gtggGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGcc                                                                                                     <  4:90640/56‑1 (MQ=255)
                                                                                                                             gtggGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGcc                                                                                                     >  3:90640/1‑56 (MQ=255)
                                                                                                                             gtggGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTAAAAGCtgg        >  3:968572/1‑149 (MQ=255)
                                                                                                                              tggGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGagca                                >  2:143538/1‑124 (MQ=255)
                                                                                                                              tggGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGagca                                <  1:143538/124‑1 (MQ=255)
                                                                                                                               gggCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGAcgttcgtt                                       >  2:468480/1‑116 (MQ=255)
                                                                                                                               gggCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGAcgttcgtt                                       <  1:468480/116‑1 (MQ=255)
                                                                                                                                ggCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTgg                                                                                                          <  3:1246283/48‑1 (MQ=255)
                                                                                                                                ggCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTgg                                                                                                          >  4:1246283/1‑48 (MQ=255)
                                                                                                                                ggCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCaa     <  4:17710/149‑1 (MQ=255)
                                                                                                                                 gCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGAc                                              >  1:63623/1‑107 (MQ=255)
                                                                                                                                 gCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGAc                                              <  2:63623/107‑1 (MQ=255)
                                                                                                                                   accaccAGTGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACcgg                                                                          >  4:122884/1‑77 (MQ=255)
                                                                                                                                   accaccAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCt  <  4:1172211/149‑1 (MQ=255)
                                                                                                                                   accaccAGCGGCCTGCATCTGGTAATTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACcgg                                                                          <  3:122884/77‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                                    
GATAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCT  >  NC_000913/4444062‑4444341

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