Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,444,192 C→G 5.5% G27G (GGC→GGG) ‡ 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,1920CG 5.5% 254.4 / 4.1 72G27G (GGC→GGG) ‡tamBtranslocation and assembly module for autotransporter export, inner membrane subunit
Reads supporting (aligned to +/- strand):  ref base C (39/29);  new base G (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 = 8.21e-01

GATAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAA  >  NC_000913/4444062‑4444338
                                                                                                                                  |                                                                                                                                                  
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                                                                                                                                          <  2:1348278/111‑3 (MQ=255)
                              aTCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGTGCCCACCAgaa                                                                                                                                          >  1:1348278/1‑109 (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                                                                                                                                               <  2:636940/83‑1 (MQ=255)
                                                     aGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGcaccac                                                                                                                                               >  1:636940/1‑83 (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                                                                                                                            <  3:214547/91‑1 (MQ=255)
                                                                aaaaaTCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTa                                                                                                                            >  4:214547/1‑91 (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                                                                                                                                      <  3:634401/58‑1 (MQ=255)
                                                                                       tcgttatcTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCt                                                                                                                                      >  4:634401/1‑58 (MQ=255)
                                                                                         gttatcTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGt                                                <  2:565269/142‑1 (MQ=255)
                                                                                         gttatcTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGt                                                >  1:565269/1‑142 (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                                                       >  4:359447/3‑104 (MQ=255)
                                                                                                                        cgcTGGTGGTGCCCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGAtctc                                                       <  3:359447/102‑1 (MQ=255)
                                                                                                                           tggtggGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGAt                                                          >  3:993093/1‑98 (MQ=255)
                                                                                                                           tggtggGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGAt                                                          <  4:993093/98‑1 (MQ=255)
                                                                                                                            ggtggGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTAAAAGCgg      <  1:1429000/149‑1 (MQ=255)
                                                                                                                             gtggGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGcc                                                                                                  >  3:90640/1‑56 (MQ=255)
                                                                                                                             gtggGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGcc                                                                                                  <  4:90640/56‑1 (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                                           <  2:63623/107‑1 (MQ=255)
                                                                                                                                 gCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGAc                                           >  1:63623/1‑107 (MQ=255)
                                                                                                                                  |                                                                                                                                                  
GATAAAGACGAACACGGGTTACAGTTTTACATCGGTCTGGGGCCAGAATTATGAGTTTATGGAAAAAAATCAGCCTCGGCGTGGTTATCGTTATCTTACTGTTGCTGGGATCGGTGGCGTTTCTGGTGGGCACCACCAGCGGCCTGCATCTGGTATTTAAAGCGGCGGATCGCTGGGTGCCAGGACTGGATATTGGCAAGGTCACCGGCGGCTGGCGCGATCTCACCTTGTCTGACGTTCGTTATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAA  >  NC_000913/4444062‑4444338

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: