Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,619,262 C→T A24V (GCC→GTC)  rhsB → Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,619,2620CT97.5% 356.8 / ‑5.2 121A24V (GCC→GTC) rhsBRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (70/47);  minor base G (3/0);  total (74/47)
Fisher's exact test for biased strand distribution p-value = 2.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

Only 100 of 121 total aligned reads displayed.

ATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCAT  >  NC_000913/3619114‑3619409
                                                                                                                                                    |                                                                                                                                                   
aTTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCg                                                                                                                                                   <  1:452861/151‑1 (MQ=255)
 ttCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCgg                                                                                                                                                  >  2:301945/1‑151 (MQ=255)
   cATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCgggg                                                                                                                                                <  2:310662/151‑1 (MQ=255)
   cATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCggg                                                                                                                                                 <  1:1863484/150‑1 (MQ=255)
   cATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCggg                                                                                                                                                 <  2:1148581/150‑1 (MQ=255)
      tgtgACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTAc                                                                                                                                             >  2:1410290/1‑151 (MQ=255)
       gtgACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACg                                                                                                                                            <  1:559687/151‑1 (MQ=255)
        tgACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACg                                                                                                                                            >  2:680780/1‑150 (MQ=255)
        tgACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGt                                                                                                                                           >  1:311131/1‑151 (MQ=255)
         gACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTa                                                                                                                                          <  2:319912/151‑1 (MQ=255)
          aCCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTAt                                                                                                                                         <  2:1112664/151‑1 (MQ=255)
            cacaAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTg                                                                                                                                       <  1:1024081/151‑1 (MQ=255)
             acaAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGGTCAGTCGGGGGACGTATTgg                                                                                                                                      >  2:1198893/1‑151 (MQ=255)
              caAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGt                                                                                                                                     <  2:1864499/151‑1 (MQ=255)
               aaGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTg                                                                                                                                    <  2:866233/151‑1 (MQ=255)
               aaGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGAATTGTGCAGGGTTCAGTCGGGGTACGTATTGGt                                                                                                                                     >  2:167641/1‑150 (MQ=255)
                 gTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTTTTTACGGTTCAGTCGGGGGACGCATTGGTGcc                                                                                                                                  >  1:934441/1‑151 (MQ=25)
                 gTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGGACGTATTGGTGcc                                                                                                                                  >  2:1731368/1‑151 (MQ=255)
                  tttttCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGccc                                                                                                                                 >  2:1300397/1‑151 (MQ=255)
                  tttttCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGccc                                                                                                                                 >  1:1866580/1‑151 (MQ=255)
                  tttttCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGGTCAGTCGGGGGACGGATTGGTGccc                                                                                                                                 >  1:607535/1‑151 (MQ=255)
                   ttttCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGGACGTATTGGTGcccc                                                                                                                                >  2:833035/1‑151 (MQ=255)
                    tttCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGcccc                                                                                                                                <  2:273925/150‑1 (MQ=255)
                      tcttcGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCa                                                                                                                              <  2:2066455/150‑1 (MQ=255)
                       cttcGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCAcc                                                                                                                            >  1:147940/1‑151 (MQ=255)
                         tcGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGGTCAGTCGGGGGACGGATTGGTGCCCCCACCgg                                                                                                                          >  2:2113791/1‑151 (MQ=255)
                            cTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCg                                                                                                                        >  1:118872/1‑150 (MQ=255)
                             tttttCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTg                                                                                                                      >  2:935475/1‑151 (MQ=255)
                              ttttCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTg                                                                                                                      <  2:1492691/150‑1 (MQ=255)
                              ttttCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAACATTGTTCAGGGGTCAGTCGGGGTACGGATTTGTGCCCCCACCGGCgggg                                                                                                                     >  2:1479076/1‑151 (MQ=21)
                                  cGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTg                                                                                                                 <  1:1117883/151‑1 (MQ=255)
                                  cGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGGACGTATTGGTGCCCCCACCGGCGTGGCCt                                                                                                                  >  2:768946/1‑150 (MQ=255)
                                     aTGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTc                                                                                                              <  2:519167/151‑1 (MQ=255)
                                      tGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAACATTGTTCAGGGTTCAGTCGGGGGACGGATTGGTGCCCCCACCGGCGTGGACCGTTCg                                                                                                             >  2:775748/1‑151 (MQ=21)
                                           aTACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGtgtg                                                                                                        >  1:626481/1‑151 (MQ=255)
                                            tACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTTTTCAGGGTTCAGTCGGGGGACGTATTGGTGCCCCCACCGGCGGGGGCTGTTCGGggtg                                                                                                        >  1:116579/1‑150 (MQ=18)
                                                gTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGcccc                                                                                                    <  2:1219515/150‑1 (MQ=255)
                                                gTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGcccc                                                                                                    <  1:1044781/150‑1 (MQ=255)
                                                gTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCcg                                                                                                   <  2:628436/151‑1 (MQ=255)
                                                     tAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCCGTTCGGTGTGCCCccgcgg                                                                                               >  2:2149687/1‑150 (MQ=255)
                                                      aaaaTAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAACATTGTTCAGGGGTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCcggcgg                                                                                               >  2:517298/1‑149 (MQ=255)
                                                      aaaaTAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGAAGCAATGTTCAGGGGTCAGTCGGGGGACGTATTGGTGCCCCCACCGGCGGGGCCTGTTCGGGGTGCCCccgcgg                                                                                               >  2:141264/1‑149 (MQ=18)
                                                        aataatAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGGTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGGTCGGGGGGGCCcggggggg                                                                                             >  2:403902/1‑147 (MQ=17)
                                                          taataGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGGTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAggg                                                                                           >  1:684331/1‑147 (MQ=255)
                                                           aataGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCCGTCGGGGTACGTATTTGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCCGCGGGGTGAc                                                                                         >  1:1928612/1‑150 (MQ=255)
                                                            ataGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGGCGGGGTACGGATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGGCCCCGCGGGGTg                                                                                           >  2:120007/1‑147 (MQ=255)
                                                            ataGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCAATGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCcggcgggggg                                                                                           >  2:1903842/1‑143 (MQ=255)
                                                             taGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGGCGGGGTACGGATTGGTGCCCCCACCGGCGTGGCCTGTTCGGGGTGCCCCGGCGGAGTGACGTc                                                                                      >  2:1704052/1‑151 (MQ=255)
                                                              aGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTcc                                                                                     <  1:499833/151‑1 (MQ=255)
                                                              aGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGGTCAGTCGGGGGACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGGGACGTcc                                                                                     >  2:1245470/1‑151 (MQ=255)
                                                               gAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGGTCAGTCGGGGTACGTATTGGGGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGGcc                                                                                     >  1:657180/1‑150 (MQ=255)
                                                                aaaaGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCg                                                                                    <  1:2029678/150‑1 (MQ=255)
                                                                aaaaGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGGACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGGGTGCCCCGGCGGGGGGACGTcc                                                                                     >  2:606892/1‑149 (MQ=255)
                                                                aaaaGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGGTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCgg                                                                                   >  2:764248/1‑151 (MQ=255)
                                                                    ggATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCAt                                                                               <  2:1668671/151‑1 (MQ=255)
                                                                      aTTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATcc                                                                             <  2:1591581/151‑1 (MQ=255)
                                                                      aTTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATcc                                                                             >  1:1723780/1‑151 (MQ=255)
                                                                      aTTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGGACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATcc                                                                             >  1:2149581/1‑151 (MQ=255)
                                                                        tttACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCgg                                                                           >  2:725191/1‑151 (MQ=255)
                                                                          tACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGt                                                                          <  2:1329560/150‑1 (MQ=255)
                                                                          tACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGGGTGCCCCGGCGGCGGGACGTCCGGCCATCCgg                                                                           >  2:1804025/1‑149 (MQ=255)
                                                                          tACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGGCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTTGGGGTGGCCCGGCGGAGTGACGTCCGGCCATCCgg                                                                           >  1:627137/1‑149 (MQ=255)
                                                                           aCGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGCAGCATTGTTTAGGGGTTAGTCGGGGTACGGATTGGTGCCCCCCCCGGCGTGGCCCGTTCGGGGTGCCCCGGGGGGGCGACGTCCCGCCATCCGGTc                                                                         >  2:1984277/1‑150 (MQ=11)
                                                                            cGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCCTCCGGTCa                                                                        >  1:968714/1‑150 (MQ=255)
                                                                            cGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGGACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCa                                                                        >  1:533709/1‑150 (MQ=255)
                                                                            cGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGGTCAGTCGGGGGACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGGGTGCCCCGGCGGAGTGACGGCCGGCCATcc                                                                             >  1:942683/1‑145 (MQ=255)
                                                                             gATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCaa                                                                       >  1:75550/1‑150 (MQ=255)
                                                                             gATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGGACGTATTGGTGCCCCCACCGGCGTGGCCTGTTTGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCaa                                                                       >  1:1439173/1‑150 (MQ=255)
                                                                                gAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGGTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATccc                                                                   >  1:110397/1‑151 (MQ=255)
                                                                                  gCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCt                                                                 <  2:1975338/151‑1 (MQ=38)
                                                                                   cGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTg                                                                <  1:1182635/151‑1 (MQ=38)
                                                                                   cGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGGGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGGCAATcccc                                                                  >  1:1781643/1‑149 (MQ=255)
                                                                                    ggAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGc                                                               <  2:1830296/151‑1 (MQ=38)
                                                                                      aaaaCCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTc                                                             >  1:2005207/1‑151 (MQ=38)
                                                                                      aaaaCCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTc                                                             <  2:1080461/151‑1 (MQ=38)
                                                                                        aaCCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCgg                                                           <  2:1316368/151‑1 (MQ=38)
                                                                                            ggcggcGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCa                                                       <  2:653691/151‑1 (MQ=38)
                                                                                              cggcgCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCaaa                                                     <  2:780822/151‑1 (MQ=37)
                                                                                                   cgTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGGGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTcc                                                >  2:1744506/1‑151 (MQ=25)
                                                                                                   cgTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCCGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTcc                                                >  1:1613093/1‑151 (MQ=25)
                                                                                                       aGGGTGACATGACGCAGTATGGCGGTAGCATTTTTCAGGGGTCAGTCGGGGTACGTATTTGTGTCCCCACCCGCGTGGCCTGTTCGGTGTGCCCCCGCGGAGTGACGGCCCGCCATCCCGGCAATCCCCCGCTTGGTGCAAAGGACCTTcc                                            >  1:1016542/1‑151 (MQ=11)
                                                                                                       aGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGGCAATCCCCTGCTCGGGGCAAAGGGCCTTcc                                            >  1:2133428/1‑151 (MQ=255)
                                                                                                         ggTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCg                                          <  2:2078638/151‑1 (MQ=25)
                                                                                                            gACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGt                                        <  2:436252/150‑1 (MQ=18)
                                                                                                            gACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTg                                       <  1:1321992/151‑1 (MQ=18)
                                                                                                                 gACGCAGTATGGCGGTAGCATTGTTCAGGGGTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGGGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAAcc                                  >  2:1889544/1‑151 (MQ=16)
                                                                                                                   cGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGa                                <  1:660719/151‑1 (MQ=18)
                                                                                                                    gCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGAc                               <  2:118872/151‑1 (MQ=18)
                                                                                                                       gTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATc                            <  1:1151612/151‑1 (MQ=18)
                                                                                                                            gCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCt                       <  1:94036/151‑1 (MQ=18)
                                                                                                                              ggTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTgc                     <  1:906672/151‑1 (MQ=18)
                                                                                                                                tAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTgccc                   >  2:1844463/1‑151 (MQ=18)
                                                                                                                                 aGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTgcccg                  <  2:834809/151‑1 (MQ=18)
                                                                                                                                 aGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTgccc                   <  2:311130/150‑1 (MQ=18)
                                                                                                                                 aGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTgccc                   <  2:195642/150‑1 (MQ=18)
                                                                                                                                   cATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTgcccggc                <  1:90498/151‑1 (MQ=18)
                                                                                                                                      tGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTgcccggccc              <  2:621327/150‑1 (MQ=18)
                                                                                                                                        ttCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGc            >  1:1380053/1‑150 (MQ=18)
                                                                                                                                           aGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGcc        >  2:1662805/1‑151 (MQ=18)
                                                                                                                                                tCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTc    >  1:1095941/1‑150 (MQ=18)
                                                                                                                                                    |                                                                                                                                                   
ATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCAT  >  NC_000913/3619114‑3619409

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