Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,874,505 T→C 15.6% 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,5050TC15.6% 130.4 / 3.4 66G305G (GGA→GGGcysNsulfate adenylyltransferase, subunit 1
Reads supporting (aligned to +/- strand):  ref base T (40/11);  new base C (4/7);  total (45/21)
Fisher's exact test for biased strand distribution p-value = 9.73e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

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

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