Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 433,518 C→A 6.0% H22N (CAT→AAT)  ribD → fused diaminohydroxyphosphoribosylaminopyrimidine deaminase/5‑amino‑6‑(5‑phosphoribosylamino)uracil reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913433,5180CA 6.0% 261.9 / 3.8 74H22N (CAT→AAT) ribDfused diaminohydroxyphosphoribosylaminopyrimidine deaminase/5‑amino‑6‑(5‑phosphoribosylamino)uracil reductase
Reads supporting (aligned to +/- strand):  ref base C (34/35);  new base A (3/2);  total (37/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01

GAAGTGCCGAGCAAGATGATTGGCAATCTGGTGATGGAGCAATTGAAAAAGCTCGATAAAGTCGCCTATATCCGTTTTGCCTCTGTCTACCGCAGTTTCGAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTG  >  NC_000913/433305‑433779
                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                      
gAAGTGCCGAGCAAGATGATTGGCAATCTGGTGATGGAGCAATTGAAAAAGCTCGATAAAGTCGCCTATATCCGTTTTGCCTCTGTCTACCGCAGTTTCGAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTCCGACATTGTCAAACAAGGCGAAAATGGAGGTGAAgg                                                                                                                                                                                                              >  2:854882/1‑272 (MQ=255)
           acaGATGTTTGGCAATCTGTTGAGTGAGCAATTGAAAAAGCTCGATAAAGTCGCCTATCTCCGTTTTGCCTCTGGCTACCGCAGTTTCGAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACGAAGCCGGGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAAc                                                                                                                                                                                                   <  1:1066503/270‑1 (MQ=255)
            aaGATTATTGGCAATCTGGTGATGGAGCAATTGAAAAAGCTCGATAAAGTCGCCTATATCCGTTTTGCCTCTGTCTACCGCAGTTTCGAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAAc                                                                                                                                                                                                   <  1:657012/271‑1 (MQ=255)
                 gatTGGCAATCTGGTGATGGAGCAATTGAAAAAGCTCGATAAAGTCGCCTATATCCGTTTTGCCTCTGTCTACCGCAGTTTCGAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAAGTGTCGGTGAAGGTTACCAACAACggtcgg                                                                                                                                                                                             >  2:307575/1‑267 (MQ=255)
                             ggTGATGGAGCAATTGAAAAAGCTCGATAAAGTCGCCTATATCCGTTTTGCCTCTGTCTACCGCAGTTTCGAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCACAGATGGCGAAATTGTCGGTGAAGGGATCCACCAAAGGGCggggggaacaacc                                                                                                                                                                                  >  2:1066502/1‑261 (MQ=255)
                                  tGGAGCAATTGAAAAAGCTCGATTAAGTCGCCTATATCCGTTTTGCCTCTGTCTACCGCAGTTTCGAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGaa                                                                                                                                                                            <  1:854883/272‑1 (MQ=255)
                                        aaTTGAAAAAGCTCGATAAAGTCGCCTATATCCGTTTTGCCTCTGTCTACCGCAGTTTCGAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTacac                                                                                                                                                                      <  1:307576/272‑1 (MQ=255)
                                                                          ttttGCCTCTGTCTACCGCAGTTTCGAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGc                                                                                                                                                                                                                                                 >  2:618384/1‑163 (MQ=255)
                                                                           tttGCCTCTGTCTACCGCAGTTTCGAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCg                                                                                                                                                                                                                                                <  1:618385/163‑1 (MQ=255)
                                                                            ttGCCTCTGTCTACCGCAGTTTCGAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCa                                                                                                                                                                                                                                                                       >  1:802502/1‑139 (MQ=255)
                                                                             tGCCTCTGTCTACCGCAGTTTCGAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCAt                                                                                                                                                                                                                                                                      <  2:802501/139‑1 (MQ=255)
                                                                                             attttCGAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGACTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACt                                                                                                                <  1:31735/270‑1 (MQ=255)
                                                                                                  cGAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCa                                                                                                                                                                                                                                                                       >  1:362330/1‑117 (MQ=255)
                                                                                                  cGAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTg                                                                                                                                                                                          >  1:647044/1‑194 (MQ=255)
                                                                                                  cGAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCCCATGCCGAAGTCAACGCGTTCCGTATGGCGGGTGAAAAAGACCAAGGGAGCAAAGCCCATGTCCCACTCGCAc                                                                                                           >  2:633208/1‑272 (MQ=255)
                                                                                                  cGAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACAt                                                                                                                                                                                  >  2:191685/1‑202 (MQ=255)
                                                                                                  cGAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACAt                                                                                                                                                                                  >  1:886785/1‑202 (MQ=255)
                                                                                                   gAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGTACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAAc                                                                                                                                                                                                   <  1:303925/184‑1 (MQ=255)
                                                                                                   gAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCAt                                                                                                                                                                                                                                                                      <  2:362329/117‑1 (MQ=255)
                                                                                                   gAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATg                                                                                                                                                                                 <  1:191685/202‑1 (MQ=255)
                                                                                                   gAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATg                                                                                                                                                                                 <  2:886784/202‑1 (MQ=255)
                                                                                                     aGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCt                                                                                                        <  1:633209/272‑1 (MQ=255)
                                                                                                     aGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCt                                                                                                        <  1:410545/272‑1 (MQ=255)
                                                                                                            aaaGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGGTACCACCAACGTGCGGGGGAAACACATGCCGAAGTAAACGCGGTGCGTATGGCGGGGGAAAACGACAAAGATGAGACCGACTATGTCACACTCGaaacaagtaaca                                                                                                  >  1:870262/1‑261 (MQ=255)
                                                                                                                                   gcgcCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGt                                                                                                                                                                                                                                                   >  2:562277/1‑104 (MQ=255)
                                                                                                                                    cgcCGGGAGGATTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGATAAAGCTGGCGAAACGAGGACTTTTTACCACCCATCCCAACCCGAATGTCGGGTg                                                                                                                                                                                                                                                  <  1:562278/104‑1 (MQ=255)
                                                                                                                                                  aGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTaca                                                                                                                                                                       >  2:160122/1‑165 (MQ=255)
                                                                                                                                                   gCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTacac                                                                                                                                                                      <  1:160122/165‑1 (MQ=255)
                                                                                                                                                          aGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAAc                                                                                                                                                                                                   >  1:446459/1‑129 (MQ=255)
                                                                                                                                                          aGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGaa                                                                                                                                                                            >  2:35594/1‑152 (MQ=255)
                                                                                                                                                          aGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTg                                                                                                                                                   >  1:942736/1‑177 (MQ=255)
                                                                                                                                                          aGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGt                                                                                          >  1:296234/1‑233 (MQ=255)
                                                                                                                                                           ggACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACg                                                                                                                                                                                                  <  2:446458/129‑1 (MQ=255)
                                                                                                                                                           ggACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAg                                                                                                                                                                           <  1:35594/152‑1 (MQ=255)
                                                                                                                                                           ggACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGa                                                                                                                                                  <  2:942735/177‑1 (MQ=255)
                                                                                                                                                                 gTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCaa                                                                                                                                                                                                    >  1:706961/1‑121 (MQ=255)
                                                                                                                                                                  tATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAAc                                                                                                                                                                                                   <  2:706960/121‑1 (MQ=255)
                                                                                                                                                                    ttACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACCCGAACCCTGTAGcc                                                                                                  >  2:189769/1‑215 (MQ=255)
                                                                                                                                                                    ttACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGGCACACTCGAACCCTGTAGCCATCATGGTCGGACGCCACCGGGCTGTTACCCCCACACCGCCCCCGGccgaccccccg                                         >  2:351167/1‑262 (MQ=255)
                                                                                                                                                                     tACAGGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCca                                                                                                 <  1:189769/215‑1 (MQ=255)
                                                                                                                                                                             gcgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACAt                                                                                                                                                                                  >  1:1084814/1‑127 (MQ=255)
                                                                                                                                                                             gcgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTAt                                                                                                                        >  1:3317/1‑184 (MQ=255)
                                                                                                                                                                             gcgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGAATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAAc                                                                                                                                                                                       >  2:925925/1‑122 (MQ=255)
                                                                                                                                                                              cgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATg                                                                                                                                                                                 <  2:1084813/127‑1 (MQ=255)
                                                                                                                                                                              cgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGAATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAAcc                                                                                                                                                                                      <  1:925926/122‑1 (MQ=255)
                                                                                                                                                                               gggCGCTAAAGCTGGCGCAACGAGGACGTTTTACTACGAATCCCAACCCGAATGTCGGGTGCGTCATTGGCAAAAATGGCGAAATTGTCGGTTAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGGACACGCGTTGCGTATGGCGGGTGAAAAAGACAAAGGTACGa                                                                                                                                >  1:1110707/1‑174 (MQ=255)
                                                                                                                                                                               gggCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAAc                                                                                                                                                                                       >  1:108452/1‑120 (MQ=255)
                                                                                                                                                                               gggCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACAt                                                                                                                                                                                  >  1:145409/1‑125 (MQ=255)
                                                                                                                                                                               gggCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTaca                                                                                                                                                                       >  2:906647/1‑136 (MQ=255)
                                                                                                                                                                               gggCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTg                                                                                                                                   >  2:437176/1‑171 (MQ=255)
                                                                                                                                                                               gggCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTTCGACCGCCTAt                                                                                                                        >  1:508271/1‑182 (MQ=255)
                                                                                                                                                                               gggCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCt                                                                                                        <  1:604243/198‑1 (MQ=255)
                                                                                                                                                                               gggCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGAATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCg                                                                                                                                                                                              >  2:104551/1‑113 (MQ=255)
                                                                                                                                                                                ggCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAAcc                                                                                                                                                                                      <  2:108452/120‑1 (MQ=255)
                                                                                                                                                                                ggCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATg                                                                                                                                                                                 <  2:145409/125‑1 (MQ=255)
                                                                                                                                                                                ggCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTacac                                                                                                                                                                      <  1:906648/136‑1 (MQ=255)
                                                                                                                                                                                ggCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATg                                                                                                                       <  2:508270/182‑1 (MQ=255)
                                                                                                                                                                                ggCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGAATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCgg                                                                                                                                                                                             <  1:104551/113‑1 (MQ=255)
                                                                                                                                                                                ggCGCTAAAGCTGGCGAAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGc                                                                                                                                  <  1:437177/171‑1 (MQ=255)
                                                                                                                                                                                 gcgcTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGAt                            <  1:351168/272‑1 (MQ=255)
                                                                                                                                                                                   gcTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTaca                                                                                                                                                                       >  1:739257/1‑132 (MQ=255)
                                                                                                                                                                                    cTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTacac                                                                                                                                                                      <  2:739256/132‑1 (MQ=255)
                                                                                                                                                                                                 aaCGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGt                                                                                                                                    >  1:125732/1‑152 (MQ=255)
                                                                                                                                                                                                  aCGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTg                                                                                                                                   <  2:125732/152‑1 (MQ=255)
                                                                                                                                                                                                      ggACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTaca                                                                                                                                                                       >  1:408416/1‑113 (MQ=255)
                                                                                                                                                                                                       gACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTacac                                                                                                                                                                      <  2:408415/113‑1 (MQ=255)
                                                                                                                                                                                                          gTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCt                                                                                                                          >  1:38666/1‑153 (MQ=255)
                                                                                                                                                                                                           ttttACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTa                                                                                                                         <  2:38666/153‑1 (MQ=255)
                                                                                                                                                                                                           tttcccaccGCATCACACCCCGATAGTCGTGTGCGTAATTGTAAGAGATGGCGAAATGGTCGGTGAAGGTTACCACAAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTGGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTg  <  2:431723/264‑1 (MQ=255)
                                                                                                                                                                                                              tACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGGCGCt   >  1:431724/1‑268 (MQ=255)
                                                                                                                                                                                                                ccACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCg                                                                                                              >  2:1106357/1‑159 (MQ=255)
                                                                                                                                                                                                                 cACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGa                                                                                                             <  1:1106358/159‑1 (MQ=255)
                                                                                                                                                                                                                  aCGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTcac                                                                                                                   >  1:205260/1‑152 (MQ=255)
                                                                                                                                                                                                                   cGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTcaca                                                                                                                  <  2:205260/152‑1 (MQ=255)
                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                      
GAAGTGCCGAGCAAGATGATTGGCAATCTGGTGATGGAGCAATTGAAAAAGCTCGATAAAGTCGCCTATATCCGTTTTGCCTCTGTCTACCGCAGTTTCGAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTG  >  NC_000913/433305‑433779

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