Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,677,182 (T)6→7 coding (1169/2061 nt) yhjG ← putative inner membrane‑anchored periplasmic AsmA family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,677,1761.T98.6% 271.1 / ‑3.0 72coding (1175/2061 nt)yhjGputative inner membrane‑anchored periplasmic AsmA family protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base T (40/31);  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 = 7.68e-01

ATGATATGAGTAGAAAGATCGCTAATCGGCAGGCTACTGCCATGCTCAATGCGCCGCCCTTTGAAGCGAACATCGGCATCCATAACGTCCCATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCGGAATCAACGCCAATCAACGGTCCCAGGTCCGCCAGCCGCAATTGCCGCGACTCGACATCACCTTCCAGTTTTGGTCGTGGCTTTCCGGT  >  NC_000913/3677032‑3677318
                                                                                                                                                 |                                                                                                                                              
atgatATGAGTAGAAAGATCGCTAATCGGCAGGCTACTGCCATGCTCAATGCGCCGCCCTTTGAAGCGAACATCGGCATCCATAACGTCCCATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCtttttt                                                                                                                                          >  1:965689/1‑151 (MQ=255)
   atatGAGTAGAAAGATCGCTAATCGGCAGGCTACTGCCATGCTCAATGCGCCGCCCTTTGAAGCGAACATCGGCATCCATAACGTCCCATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCtttttttcg                                                                                                                                       >  2:311096/1‑148 (MQ=255)
    tatGAGTAGAAAGATCGCTAATCGGCAGGCTACTGCCATGCTCAATGCGCCGCCCTTTGAAGCGAACATCGGCATCCATAACGTCCCATTTATCGGTTTCGAAGCGGTCATAAGACAGCACTTTGCCCGCAGGCTGAACGCtttttttcg                                                                                                                                       >  1:450557/1‑147 (MQ=255)
       gAGTAGAAAGATCGCTAATCGGCAGGCTACTGCCATGCTCAATGCGCCGCCCTTTGAAGCGAACATCGGCATCCATAACGTCCCATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGccc                                                                                                                                     <  2:652979/150‑1 (MQ=255)
        aGTAGAAAGATCGCTAATCGGCAGGCTACTGCCATGCTCAATGCGCCGCCCTTTGAAGCGAACATCGGCATCCATAACGTCCCATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCctt                                                                                                                                   >  1:562759/1‑151 (MQ=255)
          tAGAAAGATCGCTAATCCGCAGGCTACTGCCATGCTCACTGCGCCGCCCTTTGAAGCGAACATCGGCATCCATAACGTCCCATTTATCGGTTTCGAAGCGGGCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCcttc                                                                                                                                  >  1:302628/1‑150 (MQ=255)
             aaaGATCGCTAATCGGCAGGCTACTGCCATGCTCAATGCGCCGCCCTTTGAAGCGAACATCGGCGTCCATAACGTCCCATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTtcttct                                                                                                                              <  2:1107655/151‑1 (MQ=255)
              aaGATCGCTAATCGGCAGGCTACTGCCATGCTCAATGCGCCGCCCTTTGAAGCGAACATCGGCATCCATAACGTCCCATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTg                                                                                                                             <  2:813764/151‑1 (MQ=255)
              aaGATCGCTAATCGGCAGGCTACTGCCATGCTCAATGCGCCGCCCTTTGAAGCGAACATCGGCATCCATAACGTCCCATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTg                                                                                                                             >  2:771851/1‑151 (MQ=255)
               aGATCGCTAATCGGCAGGCTACTGCCATGCTCAATGCGCCGCCCTTTGAAGCGAACATCGGCATCCATAACGTCCCATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGt                                                                                                                            <  1:911185/151‑1 (MQ=255)
                gATCGCTAATCGGCAGGCTACTGCCATGCTCAATGCGCCGCCCTTTGAAGCGAACATCGGCATCCATAACGTCCCATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGtt                                                                                                                           <  2:29083/151‑1 (MQ=255)
                 aTCGCTAATCGGCAGGCTACTGCCATGCTCAATGCGCCGCCCTTTGAAGCGAACATCGGCATCCATAACGTCCCATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTACCCGCAGGCTGCACGCTTTTTTTCGCCCTTCTTCTGTTc                                                                                                                          >  1:330123/1‑151 (MQ=255)
                  tCGCTAATCGGCAGGCTACTGCCATGCTCAATGCGCCGCCCTTTGAAGCGAACATCGGCATCCATAACGTCCCATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGt                                                                                                                            >  2:1120646/1‑148 (MQ=255)
                  tCGCTAATCGGCAGGCTACTGCCATGCTCAATGCGCCGCCCTTTGAAGCGAACATCGGCATCCATAACGTCCCATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTc                                                                                                                          >  1:1052804/1‑150 (MQ=255)
                  tCGCTAATCGGCAGGCTACTGCCATGCTCAATGCGCCGCCCTTTGAAGCGAACATCGGCATCCATAACGTCCCATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTc                                                                                                                          >  1:176289/1‑150 (MQ=255)
                      tAATCGGCAGGCTACTGCCATGCTCAATGCGCCGCCCTTTGAAGCGAACATCGGCATCCATAACGTCCCATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGAcc                                                                                                                     >  2:1119051/1‑151 (MQ=255)
                      tAATCGGCAGGCTACTGCCATGCTCAATGCGCCGCCCTTTGAAGCGAACATCGGCATCCATAACGTCCCATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGAc                                                                                                                      >  1:332713/1‑150 (MQ=255)
                       aaTCGGCAGGCTGCTGCCATGCTCAATGCGCCGCCCTTTGAAGCGAACATCGGCATCCATAACGTCCCATTTATCGGTTTCGAAGCGGTCATAGGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGAc                                                                                                                      <  1:617936/149‑1 (MQ=255)
                          cGGCAGGCTACTGCCATGCTCAATGCGCCGCCCTTTGAAGCGAACATCGGCATCCATAACGTCCCATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGttt                                                                                                                 >  2:689520/1‑151 (MQ=255)
                          cGGCAGGCTACTGCCATGCTCAATGCGCCGCCCTTTGAAGCGAACATCGGCATCCATAACGTCCCATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGtt                                                                                                                  <  1:569813/150‑1 (MQ=255)
                           ggCAGGCTACTGCCATGCTCAATGCGCCGCCCTTTGAAGCGAACATCGGCATCCATAACGTCCCATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTc                                                                                                                >  1:109815/1‑151 (MQ=255)
                              aGGCTACTGCCATGCTCAATGCGCCGCCCTTTGAAGCGAACATCGGCATCCATAACGTCCCATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGAc                                                                                                             >  1:869373/1‑151 (MQ=255)
                                 cTACTGCCATGCTCAATGCGCCGCCCTTTGAAGCGAACATCGGCATCCATAACGTCCCATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACttt                                                                                                          >  1:543448/1‑151 (MQ=255)
                                    cTGCCATGCTCAATGCGCCGCCCTTTGAAGCGAACATCGGCATCCATAACGTCCCATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCt                                                                                                       >  1:639540/1‑151 (MQ=255)
                                      gCCATGCTCAATGCGCCGCCCTTTGAAGCGAACATCGGCATCCATAACGTCCCATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGc                                                                                                     >  1:615804/1‑151 (MQ=255)
                                        cATGCTCAATGCGCCACCCTTTGAAGCGAACATCGGCATCCATAACGTCCCATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGccc                                                                                                   >  2:829347/1‑151 (MQ=255)
                                         aTGCTCAATGCGCCGCCCTTTGAAGCGAACATCGGCATCCATAACGTCCCATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGcccc                                                                                                  <  1:89460/151‑1 (MQ=255)
                                           gCTCAATGCGCCGCCCTTTGAAGCGAACATCGGCATCCATAACGTCCCATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCt                                                                                                 <  1:615090/150‑1 (MQ=255)
                                           gCTCAATGCGCCGCCCTTTGAAGCGAACATCGGCATCCATAACGTCCCATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCt                                                                                                 <  2:1109427/150‑1 (MQ=255)
                                              cAATGCGCCGCCCTTTGAAGCGAACATCGGCATCCATAACGTCCCATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTcc                                                                                             <  1:41169/151‑1 (MQ=255)
                                              cAATGCGCCGCCCTTTGAAGCGAACATCGGCATCCATAACGTCCCATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTc                                                                                              <  2:178006/150‑1 (MQ=255)
                                                aTGCGCCGCCCTTTGAAGCGAACATCGGCATCCATAACGTCCCATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCg                                                                                           >  2:113832/1‑151 (MQ=255)
                                                aTGCGCCGCCCTTTGAAGCGAACATCGGCATCCATAACGTCCCATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTGCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCg                                                                                           <  1:637283/151‑1 (MQ=255)
                                                 tGCGCCGCCCTTTGAAGCGAACATCGGCATCCATAACGTCCCATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCgg                                                                                          >  2:755002/1‑151 (MQ=255)
                                                   cgccgcCCTTTGAAGCGAACATCGGCATCCATAACGTCCCATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCGGa                                                                                         >  2:565900/1‑150 (MQ=255)
                                                     ccgccCTTTGAAGCGAACATCGGCATCCATAACGTCCCATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCGGAATc                                                                                      <  2:1149576/151‑1 (MQ=255)
                                                         ccTTTGAAGCGAACATCGGCATCCATAACGTCCCATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCTACTTTTCTGCCCCTTTCCCGGAATCAACg                                                                                  <  1:922380/151‑1 (MQ=255)
                                                               aaGCGAACATCGGCATCCATAACGTCCCATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCGGAATCAACGCCAATc                                                                            <  1:771790/151‑1 (MQ=255)
                                                                aGCGAACATCGGCATCCATAACGTCCCATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCGGAATCAACGCCAATCa                                                                           <  2:777415/151‑1 (MQ=255)
                                                                  cGAACATCGGCATCCATAACGTCCCATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCGGAATCAACGCCAATCAAc                                                                         <  2:925743/151‑1 (MQ=255)
                                                                     aCATCGGCATCCATAACGTCCCATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCGGAATCAACGCCAATCAACGGt                                                                      <  1:189636/151‑1 (MQ=255)
                                                                          ggCATCCATAACGTCCCATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCGGAATCAACGCCAATCAACGGTCCCAg                                                                 <  2:942715/151‑1 (MQ=255)
                                                                           gCATCCATAACGTCCCATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCGGAATCAACGCCAATCAACGGTCCCAgg                                                                <  2:406845/151‑1 (MQ=255)
                                                                                 aTAACGTCCCATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCGGAATCAACGCCAATCAACGGTCCCAGGTCCGcc                                                          >  2:526458/1‑151 (MQ=255)
                                                                                 aTAACGTCCCATTTATCGGGTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCGGAATCAACGCCAATCAACGGTCCCAGGTcc                                                             >  1:837189/1‑148 (MQ=255)
                                                                                       tCCCATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCGGAATCAACGCCAATCAACGGTCCCAGGTCCGCCAGCCGc                                                    >  2:75427/1‑150 (MQ=255)
                                                                                          cATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCGGAATCAACGCCAATCAACGGTCCCAGGTCCGCCAGCCGCAAt                                                 >  2:939841/1‑151 (MQ=255)
                                                                                             ttATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCGGAATCAACGCCAATAAACGGTCCCAGGTCCGCCAGCCGCAATTGc                                              >  2:667761/1‑151 (MQ=255)
                                                                                              tATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCGGAATCAACGCCAATCAACGGTCCCAGGTCCGCCAGCCGCAATTGc                                              <  1:606959/150‑1 (MQ=255)
                                                                                                       cGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCGGAATCAACGCCAATCAACGGTCCCAGGTCCGCCAGCCGCAATTGCCGCGACTCg                                     <  2:696200/150‑1 (MQ=255)
                                                                                                       cGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCGGAATCAACGCCAATCAACGGTCCCAGGTCCGCCAGCCGCAATTGCCGCGACTCg                                     >  1:875112/1‑150 (MQ=255)
                                                                                                        gAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCGGAATCAACGCCAATCAACGGTCCCAGGTCCGCCAGCCGCAATTGCCGCGACTCGAc                                   >  1:648734/1‑151 (MQ=255)
                                                                                                         aaGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCGGAATCAACGCCAATCAACGGTCCCAGGTCCGCCAGCCGCAATTGCCGCGACTCGAc                                   >  1:694757/1‑150 (MQ=255)
                                                                                                          aGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCGGAATCAACGCCAATCAACGGTCCCAGGTCCGCCAGCCGCAATTGCCGCGACTCGACAt                                 <  2:944828/151‑1 (MQ=255)
                                                                                                               tCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCGGAATCAACGCCAATCAACGGTCCCAGGTCCGCCAGCCGCAATTGCCGCGACTCGACATCACCt                            >  1:753026/1‑151 (MQ=255)
                                                                                                                 aTAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCGGAATCAACGCCAATCAACGGTCCCAGGTCCGCCAGCCGCAATTGCCGCGACTCGACATCACCTTc                          >  2:460367/1‑151 (MQ=255)
                                                                                                                 aTAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCGGAATCAACGCCAATCAACGGTCCCAGGTCCGCCAGCCGCAATTGCCGCGACTCGACATCACCTTc                          >  2:188744/1‑151 (MQ=255)
                                                                                                                   aaGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCGGAATCAACGCCAATCAACGGTCCCAGGTCCGCCAGCCGCAATTGCCGCGACTCGACATCACCTTcc                         >  1:139357/1‑150 (MQ=255)
                                                                                                                   aaGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCGGAATCAACGCCAATCAACGGTCCCAGGTCCGCCAGCCGCAATTGCCGCGACTCGACATCACCTTcc                         <  2:798790/150‑1 (MQ=255)
                                                                                                                    aGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCGGAATCAACGCCAATCAACGGTCCCAGGTCCGCCAGCCGCAATTGCCGCGACTCGACATCACCTTcc                         >  1:798658/1‑149 (MQ=255)
                                                                                                                       cagcaCTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCGGAATCAACGCCAATCAACGGTCCCAGGTCCGCCAGCCGCAATTGCCGCGACTCGACATCACCTTCCAGttt                    >  2:1201158/1‑151 (MQ=255)
                                                                                                                         gcaCTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCGGAATCAACGCCAATCAACGGTCCCAGGTCCGCCAGCCGCAATTGCCGCGACTCGACATCACCTTCCAGTTTTg                  <  1:844540/151‑1 (MQ=255)
                                                                                                                           aCTTTGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCGGAATCAACGCCAATCAACGGTCCCAGGTCCGCCAGCCGCAATTGCCGCGACTCGACATCACCTTCCAGTTTTGGt                >  2:791453/1‑151 (MQ=255)
                                                                                                                              ttGCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCGGAATCAACGCCAATCAACGGTCCCAGGTCCGCCAGCCGCAATTGCCGCGACTCGACATCACCTTCCAGTTTTGGTCGt             >  1:1171266/1‑151 (MQ=255)
                                                                                                                                gCCCGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCGGAATCAACGCCAATCAACGGTCCCAGGTCCGCCAGCCGCAATTGCCGCGACTCGACATCACCTTCCAGTTTTGGTCGTg            <  1:699076/150‑1 (MQ=255)
                                                                                                                                  ccGCTGGCTGAGCGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCGGAATCAACGCCAATCAACGGTCCCAGGTCCGCCAGCCGCAATTGCCGCGACTCGACATCACCTTCCATTTTTGGTCGTGGCt         <  1:1029843/151‑1 (MQ=255)
                                                                                                                                   cGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCGGAATCAACGCCAATCAACGGTCCCAGGTCCGCCAGCCGCAATTGCCGCGACTCGACATCACCTTCCAGTTTTGGTCGTGGCtt        >  2:121364/1‑151 (MQ=255)
                                                                                                                                   cGCAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCGGAATCAACGCCAATCAACGGTCCCAGGTCCGCCAGCCGCAATTGCCGCGACTCGACATCACCTTCCAGTTTTGGTCGTGGCtt        >  1:917817/1‑151 (MQ=255)
                                                                                                                                     cAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCGGAATCAACGCCAATCAACGGTCCCAGGTCCGCCAGCCGCAATTGCCGCGACTCGCCATCACCTTCCAGTTTTGGTCGTGGCTTTc      >  2:671847/1‑151 (MQ=255)
                                                                                                                                     cAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCGGAATCAACGCCAATCAACGGTCCCAGGTCCGCCAGCCGCAATTGCCGCGACTCGACATCACCTTCCAGTTTTGGTCGTGGCttt       <  2:196421/150‑1 (MQ=255)
                                                                                                                                     cAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCGGAATCAACGCCAATCAACGGTCCCAGGTCCGCCAGCCGCAATTGCCGCGACTCGACATCACCTTCCAGTTTTGGTCGTGGCTTTc      <  2:504172/151‑1 (MQ=255)
                                                                                                                                     cAGGCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCGGAATCAACGCCAATCAACGGTCCCAGGTCCGCCAGCCGCAATTGCCGCGACTCGACATCACCTTCCAGTTTTGGTCGTGGCTTTc      >  2:482039/1‑151 (MQ=255)
                                                                                                                                       ggCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCGGAATCAACGCCAATCAACGGTCCCAGGTCCGCCAGCCGCAATTGCCGCGACTCGACATCACCTTCCAGTTTTGGTCGTGGCTTTCCg    <  2:639558/151‑1 (MQ=255)
                                                                                                                                        gCTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCGGAATCAACGCCAATCAACGGTCCCAGGTCCGCCAGCCGCAATTGCCGCGACTCGACATCACCTTCCAGTTTTGGTCGTGGCTTTCCgg   >  1:310049/1‑151 (MQ=255)
                                                                                                                                         cTGAACGCTTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCGGAATCAACGCCAATCAACGGTCCCAGGTCCGCCAGCCGCAATTGCCGCGACTCGACATCACCTTCCAGTTTTGGTCGTGGCTTTCCggt  >  2:178599/1‑151 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
ATGATATGAGTAGAAAGATCGCTAATCGGCAGGCTACTGCCATGCTCAATGCGCCGCCCTTTGAAGCGAACATCGGCATCCATAACGTCCCATTTATCGGTTTCGAAGCGGTCATAAGGCAGCACTTTGCCCGCAGGCTGAACGCTTTTTTCGCCCTTCTTCTGTTCAGACCGTTTCGACTTTTCTGCCCCTTTCCCGGAATCAACGCCAATCAACGGTCCCAGGTCCGCCAGCCGCAATTGCCGCGACTCGACATCACCTTCCAGTTTTGGTCGTGGCTTTCCGGT  >  NC_000913/3677032‑3677318

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