Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,846,923 C→G 5.3% R13P (CGC→CCC)  topB ← DNA topoisomerase III

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,846,9230CG 5.3% 256.7 / 4.2 76R13P (CGC→CCC) topBDNA topoisomerase III
Reads supporting (aligned to +/- strand):  ref base C (34/38);  new base G (2/2);  total (36/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

TTCCAGCGCGCATAGCGGCTGTCGTAGGCGTCTGGCTGCGCCTGCTCAAGCAGGTGACCGATACACCAGGTCACCACCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCAATTGCCGTCAGTTCAATGCCAAACTCGGCGGCTGTAGCTTTGACCTCATTTTCT  >  NC_000913/1846776‑1847070
                                                                                                                                                   |                                                                                                                                                   
ttCCAGCGCGCATAGCGGCTGTCGTAGGCGTCTGGCTGCGCCTGCTCAAGCAGGTGACCGATACACCAGGTCCCCACCTGACCATTACCGCACTCGATAAAGCCAGCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGgcgcg                                                                                                                                                    >  3:254366/1‑149 (MQ=255)
ttCCAGCGCGCATAGCGGCTGTCGTAGGCGTCTGGCTGCGCCTGCTCAAGCAGGTGACCGATACACCAGGTCACCACCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGgcgcg                                                                                                                                                    >  5:93124/1‑149 (MQ=255)
ttCCAGCGCGCATAGCGGCTGTCGTAGGCGTCTGGCTGCGCCTGCTCAAGCAGGTGACCGATACACCAGGTCACCACCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGgcgcg                                                                                                                                                    >  2:9977/1‑149 (MQ=255)
     gcgcgcATAGCGGCTGTCGTAGGCGTCTGGCTGCGCCTGCTCAAGCAGGTGACCGATACACCAGGTCACCACCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCa                                                                                                                                               <  2:8630/149‑1 (MQ=255)
             aGCGGCTGTCGTAGGCGTCTGGCTGCGCCTGCTCAAGCAGGTGACCGATACACCAGGTCACCACCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGt                                                                                                                                       >  7:14341/1‑149 (MQ=255)
                 gCTGTCGTAGGCGTCTGGCTGCGCCTGCTCAAGCAGGTGACCGATACACCAGGTCACCACCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGttttt                                                                                                                                   >  4:127958/1‑149 (MQ=255)
                  cggTCGTAGGCGTCTGGCTGCGCCTGCTCAAGCAGGTGACCGATACACCAGGTCACCACCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTc                                                                                                                                  <  6:146762/147‑1 (MQ=255)
                  cTGTCGTAGGCGTCTGGCTGCGCCTGCTAAAGCAGGTGACCGATACACCAGGTCACCACCTGACCATTACCGCACTCGATAAAGCCCTCGCCTTTCCGGTGCGGTTTGGGCAGGACATCCGCAATGGCGCGCGCCAGACTCGGTTTTTc                                                                                                                                  >  6:205987/1‑149 (MQ=255)
                    gTCGTAGGCGTCTGGCTGCGCCTGCTCAAGCAGGTGACCGATACACCAGGTCACCACCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCgg                                                                                                                                >  1:216944/1‑149 (MQ=255)
                     tCGTAGGCGTCTGGCTGCGCCTGCTCAAGCAGGTGACCGATACACCAGGTCACCACCTGACCATTACCGCACTCGATAAAACCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCg                                                                                                                                 <  2:150037/147‑1 (MQ=255)
                           gCGTCTGGCTGCGCCTGCTCAAGCAGGTGACCGATACACCAGGTCACCACCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGGTTTTCGGCAATaaa                                                                                                                         >  2:13748/1‑149 (MQ=255)
                            cGTCTGGCTGCGCCTGCTCAAGCAGGTGACCGATACACCAGGTCACCACCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAaac                                                                                                                        <  2:202232/149‑1 (MQ=255)
                                    tGCGCCTACTCAAGCAGGTGACCGATACACCAGGTCACCACCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCAt                                                                                                                >  4:137552/1‑149 (MQ=255)
                                      cgcCTGCTCAAGCAGGTGACCGATACACCAGGTCACCACCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGCTTTTCGGCAATAAACAACCGCAt                                                                                                                <  1:13748/147‑1 (MQ=255)
                                       gcCTGCTCAAGCAGGTGACCGATACACCAGGTCACCACCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCACGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTg                                                                                                              >  6:289872/1‑148 (MQ=255)
                                        ccTGCTCAAGCAGGTGACCGATACACCAGGTCACCACCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGaa                                                                                                            <  7:384054/149‑1 (MQ=255)
                                            cTCAAGCAGGTGACCGATACACCAGGTCACCACCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAAtt                                                                                                          <  3:232737/147‑1 (MQ=255)
                                            cTCAAGCAGGTGACCGATACACCAGGTCACCACCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAAtt                                                                                                          <  5:118200/147‑1 (MQ=255)
                                               aaGCAGGTGACCGATACACCAGGTCACCACCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATaa                                                                                                                          >  4:352111/1‑128 (MQ=255)
                                               aaGCAGGTGACCGATACACCAGGTCACCACCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAAc                                                                                                       <  6:93124/147‑1 (MQ=255)
                                                aGCAGGTGACCGATACACCAGGGCACCACCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGaa                                                                                                    <  1:177722/149‑1 (MQ=255)
                                                  cAGGTGACCGATACACCAGGTCACCACCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAAtc                                                                                                  <  8:14341/149‑1 (MQ=255)
                                                   aGGTGACCGATACACCAGGTCACCACCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAAtct                                                                                                 <  1:9977/149‑1 (MQ=255)
                                                   aGGTGACCGATACACCAGGTCACCACCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAAtct                                                                                                 >  5:389695/1‑149 (MQ=255)
                                                        aCCGATACACCAGGTCACCACCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAAc                                                                                             >  4:376337/1‑148 (MQ=255)
                                                        aCCGATACACCAGGTCACCACCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAAGTAACGAATCTCAAc                                                                                             >  1:333247/1‑148 (MQ=255)
                                                          cGATACACCAGGTCACCACCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCATTAAACAACCGCATTGAATTAACGAATCTCAACCAt                                                                                          <  2:216944/149‑1 (MQ=255)
                                                          cGATACACCAGGTCACCACCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCAt                                                                                          >  6:180848/1‑149 (MQ=255)
                                                           gATACACCAGGTCACCACCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATg                                                                                         <  3:193934/149‑1 (MQ=255)
                                                                acCAGGTCACCACCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACTAATCTCAACCATGGCAc                                                                                     >  8:190150/1‑148 (MQ=255)
                                                                acCAGGTCACCACCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCAc                                                                                     >  2:87931/1‑148 (MQ=255)
                                                                 ccAGGTCACCACCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGa                                                                                   <  5:180848/149‑1 (MQ=255)
                                                                  cAGGTCACCACCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGc                                                                                       <  6:375255/144‑1 (MQ=255)
                                                                  cAGGTCACCACCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGc                                                                                       >  5:375255/1‑144 (MQ=255)
                                                                  cAGGTCACCACCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGAc                                                                                  >  2:405028/1‑149 (MQ=255)
                                                                                aCCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTgg                                                                     <  6:389695/148‑1 (MQ=255)
                                                                                 ccATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAAACGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCa                                                                   >  5:44993/1‑149 (MQ=255)
                                                                                  cATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCACGCGCTGGCAc                                                                  >  7:288862/1‑149 (MQ=255)
                                                                                       ccGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGtt                                                             >  6:324467/1‑149 (MQ=255)
                                                                                       ccGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGtt                                                             <  3:137552/149‑1 (MQ=255)
                                                                                       ccGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGtc                                                             >  4:161930/1‑148 (MQ=255)
                                                                                       ccGCAATCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATTGCGCGCGCAAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGtt                                                             >  2:297842/1‑149 (MQ=255)
                                                                                        cGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTc                                                            >  5:408620/1‑149 (MQ=255)
                                                                                               gATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCACGCGCTGGCACCAGTTCGCCAAt                                                      <  8:288862/148‑1 (MQ=255)
                                                                                                 tAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGGGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCAATTGc                                                   >  4:295060/1‑149 (MQ=255)
                                                                                                 tAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGGGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCAATTGc                                                   >  4:289213/1‑149 (MQ=255)
                                                                                                   aaGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCAATTGcc                                                  <  4:254366/148‑1 (MQ=255)
                                                                                                   aaGCCATCGCCTTTCCGGTGCGGTTTGGGCACGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCAATTGcc                                                  <  5:289872/148‑1 (MQ=255)
                                                                                                      ccATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGGGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCAATTGCCGTCa                                              <  3:289213/149‑1 (MQ=255)
                                                                                                      ccATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGGGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCAATTGCCGTCa                                              <  3:295060/149‑1 (MQ=255)
                                                                                                      ccATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCAATTGCCGTCa                                              <  5:324467/149‑1 (MQ=255)
                                                                                                      ccATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCAATTGCCGTCa                                              <  1:87931/149‑1 (MQ=255)
                                                                                                        aTCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCAATTGCCGTCAGt                                            >  7:183596/1‑149 (MQ=255)
                                                                                                         tCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCAATTGCCGTCAGtt                                           >  5:33784/1‑149 (MQ=255)
                                                                                                           gCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTg                                                                      >  1:464706/1‑120 (MQ=255)
                                                                                                           gCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTg                                                                      <  2:464706/120‑1 (MQ=255)
                                                                                                            ccTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCAATTGCCGTCAGTTCaa                                        <  7:194626/149‑1 (MQ=255)
                                                                                                            ccTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCAATTGCCGTCAGTTCaa                                        <  8:183596/149‑1 (MQ=255)
                                                                                                             cTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCAATTGCCGTCAGTTCAAt                                       <  1:405028/149‑1 (MQ=255)
                                                                                                                 ccGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCTGTTTTTCGGCAATAaacaa                                                                                                                      <  6:7305/66‑1 (MQ=255)
                                                                                                                 ccGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAaacaa                                                                                                                      >  5:7305/1‑66 (MQ=255)
                                                                                                                  cGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCAATTGCCGTCAGTTCAATGCCa                                   >  5:208151/1‑148 (MQ=255)
                                                                                                                    gTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCAATTGCCGTCAGTTCAATGCCAAAc                                >  5:347321/1‑149 (MQ=255)
                                                                                                                         gTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCAATTGCCGTCAGTTCAATGCCAAACTcggc                           >  5:53276/1‑149 (MQ=255)
                                                                                                                         gTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCAATTGCCGTCAGTTCAATGCCAAACTcggc                           <  7:190150/149‑1 (MQ=255)
                                                                                                                              ggCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCAATTGCCGTCAGTTCAATGCCAAACTCGGCGGCt                       <  6:347321/148‑1 (MQ=255)
                                                                                                                              ggCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCAATTGCCGTCAGTTCAATGCCAAACTCGGCGGCTg                      >  8:394728/1‑149 (MQ=255)
                                                                                                                               gCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACcgccgc                                                                            >  3:201217/1‑94 (MQ=255)
                                                                                                                               gCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACcgccgc                                                                            <  4:201217/94‑1 (MQ=255)
                                                                                                                                 aGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACCACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCTCGCTGGCACCAGTTCGCCAATTGCCGTCAGTTCAATGCCAAACTCGGCGGCTGTAg                   <  5:205987/149‑1 (MQ=255)
                                                                                                                                   gACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAAtctc                                                                                                <  4:139795/70‑1 (MQ=255)
                                                                                                                                     cATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCAATTGCCGTCAGTTCAATGCCAAACTCGGCGGCTGTAGCttt               <  6:408620/149‑1 (MQ=255)
                                                                                                                                     cATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCAATTGCCGTCAGTTCAATGCCAAACTCGGCGGCTGTAGCttt               <  6:33784/149‑1 (MQ=255)
                                                                                                                                      aTCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCAATTGCCGTCAGTTCAATGCCAAACTCGGCGGCTGTAGCttt               >  2:49113/1‑148 (MQ=255)
                                                                                                                                      aTCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCAATTGCCGTCAGTTCAATGCCAAACTCGGCGGCTGTAGCttt               >  4:187909/1‑148 (MQ=255)
                                                                                                                                       tCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCAATTGCCGTCAGTTCAATGCCAAACTCGGCGGCTGTAGCTTTGa             <  1:49113/149‑1 (MQ=255)
                                                                                                                                         aGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCAATTGCCGTCAGTTCAATGCCAAACTCGGCGGCTGTAGCTTTGAc            <  3:127958/148‑1 (MQ=255)
                                                                                                                                          gCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTc                                                            >  3:112260/1‑99 (MQ=255)
                                                                                                                                          gCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTc                                                            <  4:112260/99‑1 (MQ=255)
                                                                                                                                                  gcgcgcCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCAATTGCCGTCAGTTCAATGCCAAACTCGGCGGCTGTAGCTTTGACCTCATTTTCt  >  3:122226/1‑149 (MQ=255)
                                                                                                                                                  gcgcgcCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCAATTGCCGTCAGTTCAATGCCAAACTCGGCGGCTGTAGCTTTGACCTCATTTTCt  <  3:376337/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TTCCAGCGCGCATAGCGGCTGTCGTAGGCGTCTGGCTGCGCCTGCTCAAGCAGGTGACCGATACACCAGGTCACCACCTGACCATTACCGCACTCGATAAAGCCATCGCCTTTCCGGTGCGGTTTGGGCAGGACATCAGCAATGGCGCGCGCCAGACTCGGTTTTTCGGCAATAAACAACCGCATTGAATTAACGAATCTCAACCATGGCACGACCGCCGCGCGCTGGCACCAGTTCGCCAATTGCCGTCAGTTCAATGCCAAACTCGGCGGCTGTAGCTTTGACCTCATTTTCT  >  NC_000913/1846776‑1847070

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