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,2620CT95.5% 170.6 / ‑4.8 66A24V (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 (33/30);  minor base G (3/0);  total (36/30)
Fisher's exact test for biased strand distribution p-value = 2.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

ATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCAT  >  NC_000913/3619114‑3619409
                                                                                                                                                    |                                                                                                                                                   
aTTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCg                                                                                                                                                   >  1:961379/1‑151 (MQ=255)
aTTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCg                                                                                                                                                   <  2:320670/151‑1 (MQ=255)
 ttCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCgg                                                                                                                                                  >  1:711772/1‑151 (MQ=255)
    aTTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGt                                                                                                                                               <  2:593910/151‑1 (MQ=255)
           ccACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTAtt                                                                                                                                        <  1:59637/151‑1 (MQ=255)
            cacaAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGGTCAGTCGGGGGACGTAtt                                                                                                                                        >  2:243777/1‑150 (MQ=255)
              caAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGt                                                                                                                                     <  2:195604/151‑1 (MQ=255)
               aaGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTgg                                                                                                                                      >  2:349287/1‑149 (MQ=255)
                         tcGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCAcc                                                                                                                            <  1:581797/149‑1 (MQ=255)
                         tcGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGGACGTATTGGTGCCCCCACCgg                                                                                                                          >  2:137725/1‑151 (MQ=255)
                              ttttCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCCGGGTTCAGTCGGGGGACGTATTGGTGCCCCCACCggggggg                                                                                                                     >  2:969739/1‑146 (MQ=255)
                              ttttCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTgg                                                                                                                     >  1:592332/1‑151 (MQ=255)
                               tttCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTgg                                                                                                                     <  1:915677/150‑1 (MQ=255)
                                  cGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGAGACATGACGCAGTATGGCGGTAGCGTTGTTCAGGGTTTAGTCGGGGTTCGTATTTGTGACCCCACCGGCGt                                                                                                                       >  1:951032/1‑145 (MQ=21)
                                       gAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCg                                                                                                             <  1:313020/150‑1 (MQ=255)
                                          gATACTGTCATGAAAATAATAGAAAAGGATTGTACGATGAGCGGAAAACCGGCGGCGCGGCAGGGTGACATGACGCAGTATGGAGGTAGCATTGTTCAGGGTTCAGTCGAGGTACGTATTGGTGCCGCCACCGGCGTGGCCTGTTCGgtgt                                                                                                         <  2:4361/151‑1 (MQ=18)
                                           aTACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGtgtg                                                                                                        <  1:594357/151‑1 (MQ=255)
                                                  cATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCcggc                                                                                                 <  1:166625/151‑1 (MQ=255)
                                                     tAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGGACGTATTGGTGCCCCCCCCGGCGGGGCCTGTTTGGGGTGCCCcggcggg                                                                                              >  2:957417/1‑150 (MQ=18)
                                                      aaaaTAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGa                                                                                              <  2:342048/150‑1 (MQ=255)
                                                         ataataGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGGACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGGGTGCCCcggcgg                                                                                               >  2:624354/1‑146 (MQ=255)
                                                            ataGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGAGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGGTCAGGCGGGGTAAGGATTTGGGCCCCCCCCGGCGTGGCCTGTTCGGGGGGCCCcggccgggggac                                                                                         >  2:217265/1‑141 (MQ=14)
                                                             taGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGGACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGGGCCCCGGCGGGGTGAc                                                                                         >  2:401091/1‑148 (MQ=255)
                                                               gAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGGGACGTcc                                                                                     >  1:630064/1‑150 (MQ=255)
                                                               gAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCg                                                                                    <  1:929802/151‑1 (MQ=255)
                                                                  aaGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGcc                                                                                 >  1:1050455/1‑151 (MQ=255)
                                                                  aaGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGcc                                                                                 <  1:2369/151‑1 (MQ=255)
                                                                   aGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCa                                                                                <  1:677120/151‑1 (MQ=255)
                                                                    ggATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGGTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGGGTGCCCCGGCGGAGGGACGTCCGGcc                                                                                 >  2:455010/1‑149 (MQ=255)
                                                                          tACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGAAGCATTGTTCAGGGCTCAATAGGGGGACGGATTGGTGCCCCCACCCGCGTGGCCTGTTTGGGGTGCCCCCGCGGAGTGAAGTCCCGCCCTCCCGt                                                                          >  2:255232/1‑150 (MQ=14)
                                                                           aCGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCa                                                                        <  1:711104/151‑1 (MQ=255)
                                                                            cGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGGATTGTTCCGGGTTCAGGCGGGGGACGGAATTGTGCCCCCCCCGGCGTGGCCTGTTTGGGGGGCCCCCGGGGGGTGGCGTCCGGCCCTCCGGGCa                                                                        >  2:267515/1‑150 (MQ=11)
                                                                             gATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCCGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCaa                                                                       >  1:1022285/1‑150 (MQ=255)
                                                                              aTGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCaa                                                                       >  1:249875/1‑149 (MQ=38)
                                                                                  gCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCCGGGGTCAGTCGGGGGACGTATTGGTGCCCCCACCGGCGGGGCCTGGTCGGTGTGCCCCGGGGGGGTGAGGTCCGGCCATcccggccaccccc                                                                  >  2:518911/1‑139 (MQ=11)
                                                                                   cGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGGTCAGTCGGGGGACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGGGTGCCCCCGCGGGGGGACGTCCGGCCCTCCCGTCAATcccc                                                                  >  2:976866/1‑149 (MQ=11)
                                                                                    ggAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGGTCAGTCGGGGGACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGGGTGCCCCGGCGGAGGGACGGCCGGCCATCCGGTCAATcccc                                                                  >  2:151337/1‑148 (MQ=14)
                                                                                     gAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCCATGTTCAGGGGTCAGGCGGGGGACGTATTGGGGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGGGACGTCCGGCCCTCCCGTCCATcccc                                                                  >  2:248854/1‑147 (MQ=11)
                                                                                         aCCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGGTCAGTCGGGGGACGTATTGGTGCCCCCACCGGCGGGGCCTGTTCGGGGTGCCCCCGCGGGGGGACGTCCGGCCCTCCGGTCAATCCCCCGCTCgg                                                           >  2:1047563/1‑150 (MQ=14)
                                                                                            ggcggcGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGGTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGGGGCCTGTTCGGTGTGCCCCGGCGGAGGGACGTCCCGCCATCCGGTCAATCCCCTGCTCGGTg                                                         >  1:225650/1‑149 (MQ=255)
                                                                                              cggcgCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGGGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGGGCaa                                                      >  2:890238/1‑150 (MQ=255)
                                                                                                  gcgTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTACGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGGCTGTTTGGTGTGCCCCGGCGGAGTGACGTTCGGCCATCCGGTCAATCCCCTGCTCGGTGGaaaag                                                   >  2:92338/1‑147 (MQ=11)
                                                                                                   cgTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGGTCAGTCGGGGTACGGATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCCGCCATCCGGTCAATCCCCCGCTCGGGGCAAAGGGcc                                                >  2:297767/1‑151 (MQ=11)
                                                                                                     tCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCCGTGTGCCCCGGCGGAGGGACGGCCGGCCATCCGGTCAATCCCCCGCTCGGTGCAAAGGGCCtt                                              >  2:950491/1‑151 (MQ=11)
                                                                                                       aGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTcc                                            <  2:257127/151‑1 (MQ=25)
                                                                                                         ggTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGGCCTTccc                                           >  1:759298/1‑150 (MQ=25)
                                                                                                            gACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTg                                       <  1:569678/151‑1 (MQ=18)
                                                                                                             aCATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTGCCCGGTGa                                      <  1:1028894/151‑1 (MQ=18)
                                                                                                                   cGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCg                                 <  2:508130/150‑1 (MQ=18)
                                                                                                                   cGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGa                                <  2:615915/151‑1 (MQ=18)
                                                                                                                     cAGTATGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACa                              <  1:680127/151‑1 (MQ=18)
                                                                                                                         aTGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCg                           <  1:802313/150‑1 (MQ=18)
                                                                                                                         aTGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCg                           <  1:837619/150‑1 (MQ=18)
                                                                                                                         aTGGCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGc                          <  2:423743/151‑1 (MQ=18)
                                                                                                                            gCGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCt                       >  1:880869/1‑151 (MQ=18)
                                                                                                                             cGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTg                      >  1:531308/1‑151 (MQ=18)
                                                                                                                             cGGTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTg                      <  1:393032/151‑1 (MQ=18)
                                                                                                                              ggTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTgc                     <  2:868008/151‑1 (MQ=18)
                                                                                                                              ggTAGCATTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCACTgc                     >  2:581037/1‑151 (MQ=18)
                                                                                                                                    aTTGTTCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTgcccggcc               >  2:250359/1‑151 (MQ=18)
                                                                                                                                         tCAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCt           <  2:963304/150‑1 (MQ=18)
                                                                                                                                          cAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGc         >  2:610333/1‑151 (MQ=18)
                                                                                                                                          cAGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGc         >  1:402108/1‑151 (MQ=18)
                                                                                                                                           aGGGTTCAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCACTGCCCGGCCCGCTGcc        <  1:581052/151‑1 (MQ=18)
                                                                                                                                                 cAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCa   <  1:638099/150‑1 (MQ=18)
                                                                                                                                                 cAGTCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCAt  <  2:863508/151‑1 (MQ=18)
                                                                                                                                                    |                                                                                                                                                   
ATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCAT  >  NC_000913/3619114‑3619409

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