Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,874,503 T→C 17.4% G305G (GGA→GGG cysN ← sulfate adenylyltransferase, subunit 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,874,5030TC17.4% 131.8 / 4.4 69G305G (GGA→GGGcysNsulfate adenylyltransferase, subunit 1
Reads supporting (aligned to +/- strand):  ref base T (40/12);  new base C (5/8);  total (46/23)
Fisher's exact test for biased strand distribution p-value = 1.55e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTT  >  NC_000913/2874285‑2874729
                                                                                                                                                                                                                              |                                                                                                                                                                                                                                      
ttaATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGCAGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCCCCACAGCCAGATCGCCCCGGCCGATGTCGATCCCATCCGCCAGCACCAGGGGGAGCGCTTCTCCGGCaaaa                                                                                                                                                                                                                               >  1:127976/1‑227 (MQ=255)
ttaATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGGGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGCCGATCTCATCCGTCAGCACCAGGGTGATCGCTCCTCCCGCaacg                                                                                                                                                                                                                               >  2:325025/1‑226 (MQ=255)
      tCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCTCCGGTAACGCTTCGTCTGCCGCCCGCACCCGATCGCCACGGCAGATGTCGACCCCATCCGTCAGCCCCACGGCGTCGCCCTCTCCGGCAAAGGcttct                                                                                                                                                                                                                         >  1:647102/1‑228 (MQ=255)
       cAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCACCCGATCCCCCCGGCCGCGGTCGATCCCATCCCTCCGCCCCCGGGTGACCGCTTCACCGGCCAAGGcttcat                                                                                                                                                                                                                        >  1:558570/1‑226 (MQ=255)
       cAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACCCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCCCCCCGGCTGATGTCGATCTCATCCGTCAGCACCAGGGCGCACGCTTTCCCCGCAAAGGgttcct                                                                                                                                                                                                                        >  1:436081/1‑226 (MQ=255)
          ccTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCATAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCCGCACCGCCGGTAACCCTTCGCCTGCCGCCAGCAGCCGACCGCCACGACGGACGTCGATCCCATCCGTCACCACCAGGGACACCGCTTCTCACGCAAAGGccacctcaa                                                                                                                                                                                                                     >  1:62187/1‑220 (MQ=255)
              aTAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCCCGGCTGATGTCGCTCTCATCCGTCAGCCCCAGCGGGATCGCTTCTCCGGCAAAGGCTCCTTCCCGATc                                                                                                                                                                                                                 >  1:619178/1‑228 (MQ=255)
              aTAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCTTCAGCACCAGGGGGATCGCTTCTCCGGCAAAGGCttcttcccccta                                                                                                                                                                                                                 >  1:644413/1‑222 (MQ=255)
              aTAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCCCGAt                                                                                                                                                                                                                  >  2:359783/1‑227 (MQ=255)
                    aaTGCCATCCACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCACCCATCCCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTCCCGCCACCCGCAGAACGCCCCCGCTGAGATCGAGCACATCCATCACCACCCGGTGTATCCCTTCCCCGGCaaatgacttttaccggcccccccc                                                                                                                                                                                                           >  1:538626/1‑207 (MQ=255)
                    aaTGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCCTAACTCTGCCCTGCAGAAAGCGGCTCTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTCACGCTTCGTCTGCCGCCAGCAGCAGATCGCCCCGGCTGACTCCGATCTCACCTGCCCGCCCCCTGCTCCTCCCTTATCCGCCAAAGCCCTCTTCCCTACCACCAt                                                                                                                                                                                                            >  1:387136/1‑227 (MQ=255)
                    aaTGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGGCGGAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGGAGGCCTCCAGCAGCACATCGCAACGGCCGATGCAGATGTCATCCGTCAGCCCCACAGGATAGCGTTCTCCGGCACAGGCCGCTCCGCGTACAccca                                                                                                                                                                                                            >  1:448002/1‑225 (MQ=255)
                     aTGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCCGATCGCCACGGCTGATGTCAACCTCTTCCGCCAGCACCACGGTGACCCCTCCTCCGGCAAAGGCTCTTTCCCGCTCACCaac                                                                                                                                                                                                           >  1:655966/1‑225 (MQ=255)
                            cAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGGGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTcc                                                                                                                                                                                                   >  2:370084/1‑227 (MQ=255)
                                     cgcgcgTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGCCGATCTCATCCGTCAGCACCAGGGCGATCACTTCTCCGGCAAAGGCTTCTTCCCGATCACCCTCCAAAGTCACGATCcgc                                                                                                                                                                                           >  2:89041/1‑227 (MQ=255)
                                     cgcgcgTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCACAAGCAGATCGCCACGGCTGATGTCGCTCTCATCCGTCAGCCCCAGGGAGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGACACGATCcgc                                                                                                                                                                                           >  2:655191/1‑227 (MQ=255)
                                       cgcgTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACCGCTGATGTCGATCTCATCCATCAGCACCAGGGAGATCGCCTCCCCGGCAAAGGCTCCTTCGCGATCACCATCCAAAGTCACGATCcccgg                                                                                                                                                                                         >  1:627940/1‑226 (MQ=255)
                                       cgcgTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAACTCGCCACGGCCGATGTCGGTCTCCTCCGTCAGCCCCAGGGTGATCGCTTCCCCGGCCAAGGCCTCTTCCCCATCACCATCAAAAGCCACCATCcgcag                                                                                                                                                                                         >  1:504295/1‑225 (MQ=255)
                                                     taataggagccagccgaaaCTCTGCCCGGGGAAAACGGTCTCTCCCCCCACCCATCCCACACCCACGGCCGCGCTCGGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCca                                                                                                                                                                         <  2:21329/210‑1 (MQ=255)
                                                      gCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCcac                                                                                                                                                                        >  2:13350/1‑227 (MQ=255)
                                                                           tGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCTCTTCTCCGGCAAAGGCTTCTTCGCGTTCACCATCAAAAGTCACGATCCCCGCGACGTTTGATTCCCCACCAAAGGGCAGAACttgaa                                                                                                                                                    >  1:316848/1‑223 (MQ=255)
                                                                            gCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCCCGGCTGATGTCGCTCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCCTCTTCGCCATCACCATCCAAAGTCACCATCCCCGCGACCTTTGATTCACAACCAGAGGGCACACCTTTGac                                                                                                                                                   >  1:380095/1‑226 (MQ=255)
                                                                            gCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCCCGATCACCATCAAAAGTCACGATCCGCCCGACGTTTGATTCCACACCCGCGGGCAGCACTTTGaca                                                                                                                                                  >  2:354343/1‑227 (MQ=255)
                                                                             cccTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCCCGCTCCCCATCAAAAGCCACAGTCCCCCCGACGTTTGATTCCCAACCCAAGGgccgtgctcttaa                                                                                                                                                   >  1:21329/1‑214 (MQ=255)
                                                                                  gagaAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTtct                                                                                                                                                                                                                                        >  2:83474/1‑137 (MQ=255)
                                                                                  gagaAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGCCAGCACTTTGACACGTTGc                                                                                                                                            >  2:246434/1‑227 (MQ=255)
                                                                                    gaAAGCGGCTCTCCCGCCACCCCTCCCACCTCCCCCGCCGCGCTCGCCCCCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGcc                                                                                                                                           <  2:267590/226‑1 (MQ=255)
                                                                                      agcgggCTGTTCCGCCTCCAATCCCCCACCCACCGCCGCGCTCTGCCCCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGAc                                                                                                                                       <  2:619178/224‑1 (MQ=255)
                                                                                       aGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGCTCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCCCCATCAAAAGTCACGATCCCCGCGCCGTTTGCTTCCCCACCACAGGGCCGCACTTTGACCCGTTGCCCCACt                                                                                                                                      >  1:541391/1‑228 (MQ=255)
                                                                                       aGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGCTCACCATCAAAAGTCACGATCCGCCCGACGTTTGATTCCACACCAGAGGGCAGCACCTTGACACAGTGCCCGACt                                                                                                                                      >  2:134029/1‑228 (MQ=255)
                                                                                       aGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCCCACCAGAGGGCAGCACTTTGAAAAGTGGCCCgcc                                                                                                                                       >  1:267590/1‑225 (MQ=255)
                                                                                             cTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGCTGTCGATCTCCTCCGTCAGCCCCAGGGTGACCGCTTCTCCGGAAAAGGCTTCTTCGCGATCACCACCAAAAGTCACCACCCGCCCGACGTTTGATTCCCCACCGAAGGGCAGCACTTTGACCCGTTGCCCCACTTTCaa                                                                                                                                 >  1:649396/1‑227 (MQ=255)
                                                                                             cTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCCCGCCGTTTGTTTCCACACCCGAGGGCAGCCCTTTGACACGTTGCCCGACTTCCa                                                                                                                                  >  2:378264/1‑227 (MQ=255)
                                                                                                                  aTCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGTCAAAGGCTTCTTCGCGATCACCATCACAAGTCACGATCCGCGCGACGTTTGATTCCACACCATAGGTCAGCACTTTGACACGTTGCCCGCCTTCCACGCGCCCGGTGGCCAGCGTTcc                                                                                                            >  2:184213/1‑228 (MQ=255)
                                                                                                                  aTCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGATCCCCGACTTCCACGCGACCGGATTCCACCGTTcc                                                                                                            >  2:347477/1‑228 (MQ=255)
                                                                                                                  aTCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGGGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCCGAGGGAAGCACTTTCACACGTTGCCCGACTTCCACGCGACCGGTAGCCGACGTTcc                                                                                                            >  1:309460/1‑228 (MQ=255)
                                                                                                                  aTCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCACCACCAGGGGGATCGCTTCTCCGGCAAAGGCTTCTTCCCGATCACCATCAAAAGTCACGATCCGCCCGACGTTTGATTCCCACCCCAAGGGCAGCACTTTTACACGTTTCCCGCGTTCCACCCACCCGGAGTCCAGCGTTcc                                                                                                            >  1:596749/1‑228 (MQ=255)
                                                                                                                         gACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGGGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCACGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCCACCGGATGCCCGCGTTCCGGccgaa                                                                                                     >  1:526843/1‑224 (MQ=255)
                                                                                                                            gcgcTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGCTCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCTCCGTTTGATTCCACACCCGAGGGCCGCACTTTGCCCCGTTGCCCGCCTTCCCCGCGACCGATGGCCACCGTTCCGCCGTAcccc                                                                                                  >  2:198785/1‑227 (MQ=255)
                                                                                                                              gttcggccccgccgGTACCGCTTCGGCTGCCGCCAGCAGCAGTTCGCCAAGGCTGATTACGATCCCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACg                                                                                                <  2:316848/220‑1 (MQ=255)
                                                                                                                                         ccgGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGCAAATCGAGAt                                                                                     >  2:88515/1‑228 (MQ=255)
                                                                                                                                          cgGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCAACCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCCCGATCCGCGCGACGTTTGATTCCACACCAGAGGGAAGCACTTTGACACGTTGCCCGACTTCCACCAGACCGGGTGCCAGCGTTCCGGGGCAACACCGAAAAtgcaaat                                                                                     >  1:433251/1‑221 (MQ=255)
                                                                                                                                             tAACGCTTCGTCTGCCGGCAGGAGCAGATCGCCCCGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGAACGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCgg                                                                                 <  1:13350/228‑1 (MQ=255)
                                                                                                                                                  cTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGGCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGAt                                                                            <  2:436081/228‑1 (MQ=255)
                                                                                                                                                    tccgcTGCCTCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTa                                                                          <  2:433251/224‑1 (MQ=255)
                                                                                                                                                    tCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGGGATTCGGGCGAtt                                                                           >  2:199347/1‑227 (MQ=255)
                                                                                                                                                         tGCCGCCAGCCGCAGATCGCCACGGCTGATGGCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTa                                                                     <  2:127976/228‑1 (MQ=255)
                                                                                                                                                         tGCCGCCAGCAGCAGATGGCCCCGGCTGATGCCGATCTCATCCGTCAGCTCCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTa                                                                     <  2:309460/228‑1 (MQ=255)
                                                                                                                                                          gccggcAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCCCCGGCAAAGGCTTCTGCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTa                                                                     <  2:380095/227‑1 (MQ=255)
                                                                                                                                                          gccgccAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTa                                                                     <  2:373452/227‑1 (MQ=255)
                                                                                                                                                              ccAGCAGCAGACCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCa                                                                <  2:196289/228‑1 (MQ=255)
                                                                                                                                                                        gATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCCCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGGGATTAACCGAATGCCCCGGGAAgc                                                       >  1:638406/1‑228 (MQ=255)
                                                                                                                                                                           cGCCACGGCTGATGTCGATTGCATCCGTCAGCACCAGGGTGATCGCTTCGCCGGCAAAGGCTTCCTCGCGATCACCATCAAAAGTCACGATCCGCTCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCGCa                                                   <  2:596749/228‑1 (MQ=255)
                                                                                                                                                                                    ttattcgcattcccgccGGCAGCCCCAGGGTGATCGCTTCCCCGGCAAAGGCTTCCTCGCGATCACCATCAAAAGTCACGATCCGCGCGGCGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGa                                           <  1:88515/212‑1 (MQ=255)
                                                                                                                                                                                         acggatgagatcGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTCACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACTGAAATCGACATTCGGGCGATTAAGGTCCTTGACCGGAAGAGCCATTGGCTGAGCCTcc                                    >  1:374057/12‑228 (MQ=255)
                                                                                                                                                                                             tctcCTCCCTCCGCACCAGGGTGATCGGTTCCTCGGCAAAGGCTTTCTCGCGATCACCAACAAAAGTCACCATCCGGCGCACTTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATccacca                                <  1:359783/228‑1 (MQ=255)
                                                                                                                                                                                             tctcACCCGTCAGCACCAGGGTGAGCGCTTCTCCGGCAAAGGCTGCTTCGCGATCACCAGCAAAAGTCACGAACCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATccacca                                <  1:89041/228‑1 (MQ=255)
                                                                                                                                                                                              cccacccGTCAACACCATGGTTATCGCTCTACCGGAAAAGGCTTTTCCGCTATCACGATCAAAGGTCACGATCCGCGCGCTGTTTGATCCCACACCAGAGGGCCGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATccacca                                <  1:347477/222‑1 (MQ=255)
                                                                                                                                                                                                            ccAGGGTTATCGCTTCCCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCCACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCa                 <  2:644413/228‑1 (MQ=255)
                                                                                                                                                                                                             cAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAACTCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCa                 >  2:122286/1‑227 (MQ=255)
                                                                                                                                                                                                              aGGGTGATCGCTTCCCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACg               <  2:37861/228‑1 (MQ=255)
                                                                                                                                                                                                              aGGGTGATCGCTTCCCCCGCAAAGGCTTCTTCCCGATCCCCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACCCCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACg               <  2:638406/228‑1 (MQ=255)
                                                                                                                                                                                                               gggTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCTGGCGATTAACGTACTGCACCGGGAAGGGCATTGTCCGAGCATCCACCACTCGCGGGAACCCCACgg              >  1:89142/1‑228 (MQ=255)
                                                                                                                                                                                                                    aTCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCTGGCGATTAACGTACTGCACCGGGAATCGCATTGGCTGTGCATCCCCCACTCGCTGGATATCCACGGTTTcc         >  2:21444/1‑228 (MQ=255)
                                                                                                                                                                                                                    aTCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTc          >  2:481791/1‑227 (MQ=255)
                                                                                                                                                                                                                    aTCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGAAATTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACTACCTGCACCCGGAACCGCATTGGCACACATACCACCACTCGCTGGTtccacccgtgttc          >  1:92298/1‑216 (MQ=255)
                                                                                                                                                                                                                     tCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCCCCGCGACGCGCATTGGCTGGGCATCCCCCACTCCCTGACTACCCCCGGTTTccc        >  1:46976/1‑227 (MQ=255)
                                                                                                                                                                                                                      tcccTTCCCCGGCAAAGGTTTCTTCGCGATCCCCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGCCACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCa        <  2:243018/225‑1 (MQ=255)
                                                                                                                                                                                                                      cGCTTCCCCGGCAAAGGCTTCTTCGCGATCCCCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCa        <  2:541391/227‑1 (MQ=255)
                                                                                                                                                                                                                      cGCTTCACCGGCAAAGGCTTCTGCGCGATCACCATCAAAAGCCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCa                 <  2:46976/218‑1 (MQ=255)
                                                                                                                                                                                                                             ctcCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTcc                                                                                                                                                                          >  1:247456/1‑62 (MQ=255)
                                                                                                                                                                                                                             ctcCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGCGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAACTCGCGATTCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCAACCACTCGCTGGATCTCCCCGGGTTcaaaacctc  >  1:551590/1‑220 (MQ=255)
                                                                                                                                                                                                                             ctcCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGCCTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGCTTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCAACCACTCGCTGGATCTCCACGGTTTcacgaactt  >  2:164457/1‑220 (MQ=255)
                                                                                                                                                                                                                              tcCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCca                                                                                                                                                                         <  2:247456/62‑1 (MQ=255)
                                                                                                                                                                                                                              |                                                                                                                                                                                                                                      
TTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAGATTCGGGCGATTAACGTACTGCACCGGGAAGCGCATTGGCTGAGCATCCACCACTCGCTGGATCTCCACGGTTTCCAGCACTT  >  NC_000913/2874285‑2874729

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