Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,226,119 C→G 6.1% intergenic (‑39/‑481) minC ← / → ycgJ inhibitor of FtsZ ring polymerization/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,226,1190CG 6.1% 458.2 / 15.3 131intergenic (‑39/‑481)minC/ycgJinhibitor of FtsZ ring polymerization/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (56/67);  new base G (4/4);  total (60/71)
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.25e-01

Only 100 of 136 total aligned reads displayed.

CCAATAACCCGCAAACCGGTTGCCGAAACCGCCTTATGCATCGCTGACCAGTTTACCGGGTCTTCCAGTGCACTGACGTTGAGTACAACGGGGGCATGTTTTAAAAATGCGGGGGCCTGAGCGATTTTGTCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGTCAATCAACCCATTGCCGTCGTACTTTTGATTGTTCTTATTTACGCTTCCTTTTTCCGCACCCTAATCGTGACAGCATCTCCAGACCTGTTTAAATATCCAGTATTGAGGCAAGTGTA  >  NC_000913/1225851‑1226371
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                            
ccAACAACCCGCAAACCGGTGGCCCAAACCGCCTTATGCATCGCTGACCAGTTTACCGGGTCTTCCAGTGCACTGACGTTGAGTACAACGGGGGCATGTTTTAAAAATGCGGGGGCCTGAGCGATTTTGTCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCgatg                                                                                                                                                                                                                                                        <  2:324122/275‑1 (MQ=255)
         cGCAAACCGGTTGCCGAAACCGCCTTATGCATCGCTGACCAGTTTACCGGGTCTTCCAGTGCACTGACGTTGAGTACAACGGGGGCATGTTTTAAAAATGCGGGGGCCTGAGCGATTTTGTCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGa                                                                                                                                                                                                                                               <  1:21452/275‑1 (MQ=255)
          gccaaCCGGTTGCCTAAACCGCCTTATGCATCGCCGACCAGTTTACCGGGTCTTCCAGTGCACTGACGTTGAGTACAACGGGGGCATGTTTTAAAAATGCGGGGGCCTGACCGATTTTGTCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGaa                                                                                                                                                                                                                                              <  4:17309/272‑1 (MQ=255)
                 ggTTGTCGAAACCGCCTTATGCATCGCTGACCAGTTTACCGGGTCTTCCAGTGCACTGACGTTGAGTAGAACGGGGGCATGTTTTAAAAATGCGGGGGCCTGAGCGATTTTGTCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACtata                                                                                                                                                                                                                                       <  2:64636/275‑1 (MQ=255)
                     gCCGAACCCGCCTTAGCCAACGCTTACCAGTTTCCCGGGTCTTCCAGTGCACTGACGTTGATGACAACGGGGGCATGNTTTAAAAATGCGTGGGCCTTAGCGATTTTGTCTTCCAGCGCCTGATGGATACCCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGc                                                                                                                                                                                                                                   <  3:121406/275‑1 (MQ=255)
                     gCCGAAACCGCCTTATGGAGCGCGGACCAGTTTACCGGGTCTTCCGTTGCACTGACGTGGAGTACAACGGGGGCAGGTTTTAAAAAAGCGGGGGCCTGAGCGATTTTGTCTTCCAGCGCCTGATGGTTAACCTTAGGTTCTGCCTCATGGAGATGAACAACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGc                                                                                                                                                                                                                                   <  4:469383/275‑1 (MQ=255)
                           aCCGCCTTAGGGATCGCTGCCCAGTTTTCCGGGTATTCCAGTGAACGGACGTAGAGTACAACGGGGGCATGTTTTAAAAAGTCGGGGGCCTGGGCGATTTTGTCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTta                                                                                                                                                                                                                             <  4:136325/275‑1 (MQ=255)
                               ccTTATGCATCGCTGACCAGTTTACCGGGTCTTCCAGTGCACTGACGTTGAGTACAACGGGGGCATGTTTTAAAAATGCGGGGGCCTGAGCTATTTTGTCTTCCAGCGCCTGATGGATAACTTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCGACTGCCTTTAAGCTAGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGt                                                                                                                                                                                                                         <  1:95670/275‑1 (MQ=255)
                                                          tgttttccgtttgaCTGATGTTGATGACNANGGGGGCATTTTTTAAAAATGCGGGGGCCTGAGCGATGTTGTCTTCNAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTNTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCt                                                                                                                                                                                              <  3:47171/262‑1 (MQ=255)
                                                          ggTCATCCAGGGCACTGGAGGTGAGTACAACGGGGTCATGTTTTCAAAAGGCGGGGGCCTGAGCGCTTGTGGTTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGAGGAACCACAGATAACGTGAAGCTACTGCCTTTAAGCTCGATTGTCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCt                                                                                                                                                                                              <  4:243087/275‑1 (MQ=255)
                                                                     gCACTGACGTTGAGTACAACGGGGGCATGTTTTAAAAATGCGGGGGCCTGAGCGATTTTGTCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCCATTAGCTATTAATCATCGCCAGCGCGCGATTATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGTGGCAAGTTACCCTCCACTTTTTTc                                                                                                                                                                                   >  4:362281/1‑275 (MQ=255)
                                                                          gACGTTGAGTACAACGGGGGCATGTTTTAAAAATGCGGGGGCCTGAGCGATTTTGTCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGtt                                                                                                                                                                                                                              <  2:192985/227‑1 (MQ=255)
                                                                                         gggggAATGTTTTAAAAATTCGGGGGCCTGAGCGATTTTGTCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCtg                                                                                                                                                                <  2:362282/274‑1 (MQ=255)
                                                                                                     tAAAAATGCGGGGGCCTGAGCGATTTTGTCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCa                                                                                                                                                                                                  >  3:135991/1‑228 (MQ=255)
                                                                                                         aaTGCGGGGGCCTGAGCGATTTTGTCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCt                                                                                                                                                                                              <  1:135991/228‑1 (MQ=255)
                                                                                                                  gCCTGAGCGATTTTGTCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCtt                                                                                                                                      <  2:383099/275‑1 (MQ=255)
                                                                                                                  gCCTGAGCGATTTTGTCTTCCAGCGCCTGATAGATAACCGTAGGTTCTGCCTCATGCGGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTTGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGAAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCtt                                                                                                                                      <  2:150128/275‑1 (MQ=255)
                                                                                                                                 tCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCgatgat                                                                                                                                                                                                                                                      <  5:721052/148‑1 (MQ=255)
                                                                                                                                      cAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATa                                                                                                                                          >  4:413075/1‑251 (MQ=255)
                                                                                                                                          gcCTGATGGATAACCTTAGGTTCTGCCTAATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCtt                                                                                                                                      <  2:413076/251‑1 (MQ=255)
                                                                                                                                                      aCCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGTTGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTTTATTGAGGCAAGTCACCCCCCCAATTATTCAAATTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTTTCAATAAACCAATTGCcgtc                                                                                                    >  4:256165/1‑273 (MQ=255)
                                                                                                                                                        cTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGcgcgcg                                                                                                                                                                                                                                                           <  8:937157/120‑1 (MQ=255)
                                                                                                                                                        cTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGcgcgcg                                                                                                                                                                                                                                                           >  7:937157/1‑120 (MQ=255)
                                                                                                                                                          tAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTta                                                                                                                                                                                                                             >  8:715176/1‑148 (MQ=255)
                                                                                                                                                            ggTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCGCGCTGGCgatgat                                                                                                                                                                                                                                                      >  7:367017/1‑121 (MQ=255)
                                                                                                                                                            ggTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCGCGCTGGCgatgat                                                                                                                                                                                                                                                      <  8:367017/121‑1 (MQ=255)
                                                                                                                                                                 tGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCa                                                                                                                                                                                                  >  2:318302/1‑168 (MQ=255)
                                                                                                                                                                 tGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGTCAATCAACCCATTGCCGTCGTACTTTTGATTg                                                                                       >  2:416405/1‑275 (MQ=255)
                                                                                                                                                                  gCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGAc                                                                                                                                                                                                                                           >  2:259161/1‑126 (MQ=255)
                                                                                                                                                                   ccTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGtt                                                                                                                                                                                                                                                   <  6:96481/117‑1 (MQ=255)
                                                                                                                                                                   ccTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGTCAATCAACCCATTGCCGTCGTACTTTTGATTGtt                                                                                     <  4:159250/275‑1 (MQ=255)
                                                                                                                                                                     tCATGCAGATGAACCACAGATAAAGTTAANCTACTGCCNNTAAGCTCGATTGGCGTGTTTGACNTCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCt                                                                                                                                                                                              <  3:30431/168‑1 (MQ=255)
                                                                                                                                                                     tCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCGCGCTGGCgatgatt                                                                                                                                                                                                                                                     >  7:805134/1‑112 (MQ=255)
                                                                                                                                                                     tCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCGCGCTGGCgatgatt                                                                                                                                                                                                                                                     <  8:805134/113‑2 (MQ=255)
                                                                                                                                                                     tCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCt                                                                                                                                                                                              <  4:318301/168‑1 (MQ=255)
                                                                                                                                                                      cATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACtata                                                                                                                                                                                                                                       <  4:259160/126‑1 (MQ=255)
                                                                                                                                                                       aTGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATtat                                                                                                                                                                                                                <  7:715176/148‑1 (MQ=255)
                                                                                                                                                                        tGCAGATGAACCACAGATAAAGTGACGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCATGGCCTTACTCAATTACCTATTAATCATCGCCAGCGCGGGATGATGTTCCGAAGACTATAAGGCATGTTAGAGTCTGGTTTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATACCTTAAAAAATACCATTGTCATTCACACCATTGCCGTCGTACTTTTGATTGTTCtatt                                                                                >  4:78374/1‑272 (MQ=255)
                                                                                                                                                                        tGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATCGCTTAAAAAATACCCTTGTCAATCAACCCATTCCCGTCGTACTTTTGATTTTTCTTAt                                                                                >  2:107986/1‑275 (MQ=255)
                                                                                                                                                                             aTGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGTCAATCAACCCATTGCCGTCGTACTTTTGATTGTTCTTATTTAc                                                                            >  4:382639/1‑274 (MQ=255)
                                                                                                                                                                                  ccACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCa                                                                                                                                                                                                  >  2:547519/1‑151 (MQ=255)
                                                                                                                                                                                      aGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCt                                                                                                                                                                                              <  4:547518/151‑1 (MQ=255)
                                                                                                                                                                                         tAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCt                                                                                                                                                                                              >  7:624410/1‑148 (MQ=255)
                                                                                                                                                                                          aaaGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCga                                                                                                                                                                                                                                                          <  8:76958/87‑1 (MQ=255)
                                                                                                                                                                                          aaaGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCga                                                                                                                                                                                                                                                          >  7:76958/1‑87 (MQ=255)
                                                                                                                                                                                            aGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCt                                                                                                                                       <  4:320303/200‑1 (MQ=255)
                                                                                                                                                                                             gTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCagag                                                                                                                                                                               >  3:60766/1‑159 (MQ=255)
                                                                                                                                                                                                aaGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCAttt                                                                                                                                                                                       >  5:781742/1‑148 (MQ=255)
                                                                                                                                                                                                 aGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTaaa                                                                                                                                                                           <  1:60766/159‑1 (MQ=255)
                                                                                                                                                                                                   cTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATTCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGAGCGCGATGGTGTTCCGa                                                                                                                                                                                                                                               <  7:795360/89‑1 (MQ=255)
                                                                                                                                                                                                   cTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGa                                                                                                                                                                                                                                               >  8:795360/1‑89 (MQ=255)
                                                                                                                                                                                                    tACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTc                                                                                                                                                                                   <  8:624410/148‑1 (MQ=255)
                                                                                                                                                                                                       tGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAg                                                                                                                                                                                                          <  6:775254/122‑1 (MQ=255)
                                                                                                                                                                                                       tGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGTCAATCAACCCATTGCCGTCGTACTTTTGATTGTTCTTATTTACGCTTCCTTTTTCCGCTCCCTAATCGt                                                  >  1:19249/1‑274 (MQ=255)
                                                                                                                                                                                                          cTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGatgatg                                                                                                                                                                                                                                                     <  6:331604/76‑1 (MQ=255)
                                                                                                                                                                                                                     aTTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCGCGCTGGCgatga                                                                                                                                                                                                                                                       >  8:1095758/1‑63 (MQ=37)
                                                                                                                                                                                                                     aTTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCGCGCTGGCgatga                                                                                                                                                                                                                                                       <  7:1095758/63‑1 (MQ=37)
                                                                                                                                                                                                                     aTTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACtata                                                                                                                                                                                                                                       >  5:94895/1‑79 (MQ=255)
                                                                                                                                                                                                                     aTTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACtata                                                                                                                                                                                                                                       <  6:94895/79‑1 (MQ=255)
                                                                                                                                                                                                                      ttGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGa                                                                                                                                                                                                                                               <  6:124997/70‑1 (MQ=255)
                                                                                                                                                                                                                      ttGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGa                                                                                                                                                                                                                                               >  5:124997/1‑70 (MQ=255)
                                                                                                                                                                                                                       tGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGa                                                                                                                                                                                                                    <  8:615188/96‑1 (MQ=255)
                                                                                                                                                                                                                          cGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCGCGCTGGCgatgat                                                                                                                                                                                                                                                      >  7:610747/1‑59 (MQ=25)
                                                                                                                                                                                                                          cGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCtg                                                                                                                                                                >  8:1113418/1‑145 (MQ=255)
                                                                                                                                                                                                                          cGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCtg                                                                                                                                                                <  7:1113418/145‑1 (MQ=255)
                                                                                                                                                                                                                             gtTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCa                                                                                                                                                                                                                                  <  8:561948/76‑1 (MQ=255)
                                                                                                                                                                                                                             gtTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCa                                                                                                                                                                                                                                  >  7:561948/1‑76 (MQ=255)
                                                                                                                                                                                                                              tttGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGTCAATCAACCCATTGCCGTCGTACTTTTGATTGTTCTTAtt                                                                               >  2:187868/1‑222 (MQ=255)
                                                                                                                                                                                                                               ttGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGtt                                                                                                                                                                                                                              >  7:97830/1‑78 (MQ=255)
                                                                                                                                                                                                                               ttGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCAcc                                                                                                                                                                                                >  5:521596/1‑108 (MQ=255)
                                                                                                                                                                                                                               ttGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCAcc                                                                                                                                                                                                <  6:1029854/108‑1 (MQ=255)
                                                                                                                                                                                                                               ttGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCAcc                                                                                                                                                                                                >  5:1029854/1‑108 (MQ=255)
                                                                                                                                                                                                                               ttGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCAcc                                                                                                                                                                                                <  6:521596/108‑1 (MQ=255)
                                                                                                                                                                                                                                  aCATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGTCAATCAACCCATTGCCGTCGTACTTTTGATTGTTCTTATTTACg                                                                           <  4:187867/222‑1 (MQ=255)
                                                                                                                                                                                                                                        tGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTc                                                                                                                                                                                                                        <  5:931417/75‑1 (MQ=255)
                                                                                                                                                                                                                                        tGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTc                                                                                                                                                                                                                        >  6:931417/1‑75 (MQ=255)
                                                                                                                                                                                                                                          gCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTa                                                                                                                                                                                      <  6:437302/107‑1 (MQ=255)
                                                                                                                                                                                                                                          gCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTa                                                                                                                                                                                      >  5:437302/1‑107 (MQ=255)
                                                                                                                                                                                                                                              tACTCAATTAGCTATTAATCATCGCCAGGGCGCGATGATGTTCCGAAGACTATAAg                                                                                                                                                                                                                                     <  8:628240/56‑1 (MQ=255)
                                                                                                                                                                                                                                              tACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAg                                                                                                                                                                                                                                     >  7:628240/1‑56 (MQ=255)
                                                                                                                                                                                                                                               aCTCAATTAGCTTTTAAGCATCGCCAGTGCGCGATGAGGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGTCAATCAACCCATTGCCGTCGTACTTTTGATTGTTCTTATTTACGCTTCCTTTTTCCGCACCCTAATCGTGACAGCATCTCCAGACCTGTTTAAATATCCAGTATTGAGGc         <  4:107985/275‑1 (MQ=255)
                                                                                                                                                                                                                                               aCTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGTCAATCAACCCATTGCCGTCGTACTTTTGATTGTTCTTATTTACGCTTCCTTTTTCCGCACCCTAATCGTGACAGCATCTCCAGACCTGTTTAAATATCCAGTATTGAGGc         >  2:194934/1‑275 (MQ=255)
                                                                                                                                                                                                                                                cTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTa                                                                                                                                     >  2:549750/1‑150 (MQ=255)
                                                                                                                                                                                                                                                    aTTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTaaaaa                                                                                                                                 <  4:549749/150‑1 (MQ=255)
                                                                                                                                                                                                                                                     ttAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGTCAATCAACCCATTGCCGTCGTACTTTTGATTGTTCTTATTTACGCTTCCTTTTTCCGCACCCTAATCGTGACAGCATCTcc                                      >  4:113768/1‑240 (MQ=255)
                                                                                                                                                                                                                                                        gCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGTCAATCAACCCATTGCCGTCGTACTTTTGATTGTTCTTATTTACGCTTCCTTTTTCCGCACCCTAATCGTGACAGCATCTCCAGACCTGTTTAAATATCCAGTATTGAGGCAAGTGTa  <  1:96480/273‑1 (MQ=255)
                                                                                                                                                                                                                                                         cTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCtg                                                                                                                                                                <  8:955883/114‑1 (MQ=255)
                                                                                                                                                                                                                                                         cTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCtg                                                                                                                                                                >  7:955883/1‑114 (MQ=255)
                                                                                                                                                                                                                                                         cTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGTCAATCAACCCATTGCCGTCGTACTTTTGATTGTTCTTATTTACGCTTCCTTTTTCCGCACCCTAATCGTGACAGCATCTCCAGAc                                  <  2:113769/240‑1 (MQ=255)
                                                                                                                                                                                                                                                            ttAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCa                                                                                                                                            >  7:920458/1‑131 (MQ=255)
                                                                                                                                                                                                                                                            ttAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCa                                                                                                                                            <  8:920458/131‑1 (MQ=255)
                                                                                                                                                                                                                                                                tcatcGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCaa                                                                                                                                                                                                      >  7:754403/1‑69 (MQ=255)
                                                                                                                                                                                                                                                                tcatcGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCaa                                                                                                                                                                                                      <  8:754403/69‑1 (MQ=255)
                                                                                                                                                                                                                                                                  atcGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCt                                                                                                                                                                 <  8:1112830/104‑1 (MQ=255)
                                                                                                                                                                                                                                                                  atcGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCt                                                                                                                                                                 >  7:1112830/1‑104 (MQ=255)
                                                                                                                                                                                                                                                                   tcGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAAt                                                                                                                                                        >  6:50943/1‑112 (MQ=255)
                                                                                                                                                                                                                                                                      ccAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGTCAATCAACCCATTGCCGTCGTACTTTTGATTGTTCTTATTTACGCTTCCTTTTTCCGCACCCTAATCGTGACAGCATCTCCAGACCTGt                              >  2:360277/1‑231 (MQ=255)
                                                                                                                                                                                                                                                                       cAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCag                                                                                                                                                                                 <  5:772514/83‑1 (MQ=255)
                                                                                                                                                                                                                                                                       cAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCag                                                                                                                                                                                 >  6:772514/1‑83 (MQ=255)
                                                                                                                                                                                                                                                                         gcgcgcGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAAt                                                                                                                               <  5:424557/131‑1 (MQ=255)
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                            
CCAATAACCCGCAAACCGGTTGCCGAAACCGCCTTATGCATCGCTGACCAGTTTACCGGGTCTTCCAGTGCACTGACGTTGAGTACAACGGGGGCATGTTTTAAAAATGCGGGGGCCTGAGCGATTTTGTCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGTCAATCAACCCATTGCCGTCGTACTTTTGATTGTTCTTATTTACGCTTCCTTTTTCCGCACCCTAATCGTGACAGCATCTCCAGACCTGTTTAAATATCCAGTATTGAGGCAAGTGTA  >  NC_000913/1225851‑1226371

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