Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,449,779 G→A S35S (AGC→AGT rpsC ← 30S ribosomal subunit protein S3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,449,7790GA97.3% 262.9 / ‑6.6 74S35S (AGC→AGTrpsC30S ribosomal subunit protein S3
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (32/40);  minor base C (2/0);  total (34/40)
Fisher's exact test for biased strand distribution p-value = 2.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCAACCAGTTTTGCGTCCAGTTCAGGCTTACGAACTTCGGCGATGTTGATCTGTGCAGGAACGCCAGCGATGTCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCGCTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTA  >  NC_000913/3449533‑3450021
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  
aGCAACCAGTTTTGCGTCCAGTTCAGGCTTACGAACTTCGGCGATGTTGATCTGTGCAGGAACGCCAGCGATGTCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGt                                                                                                                                                                                                                                                >  1:166411/1‑251 (MQ=255)
       aGTTTTGCGTCCAGTTCAGGCTTACGAACTTCGGCGATGTTGATCTGTGCAGGAACGCCAGCGATGGCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCGTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCATGTACTGACGTACTTTAAAATCACTGTCCAGGtt                                                                                                                                                                                                                                         <  2:279864/251‑1 (MQ=255)
       aGTTTTGCGTCCAGTTCAGGCTTACGAACGTCGGCGATGTTGATCTGTGCAGGAACGCCAGCGATGTCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCGTTAGCCAGTTCCTTAGTCATGTACTGACGTACTTTAAAATCACTGTCCAGGtt                                                                                                                                                                                                                                         <  2:327986/251‑1 (MQ=255)
        gTTTTGCGTCCAGTTCAGGCTTACGAACTTCGGCGATGTTGATCTGTGCAGGAACGCCAGCGATGTCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAACATCACTGTCCAGGTTg                                                                                                                                                                                                                                        >  1:465430/1‑251 (MQ=255)
           ttGCGTCCAGTTCAGGCTTACGAACTTCGGCGATGTTGATCTGTGCAGGAACGCCAGCGATGGCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATTCCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTc                                                                                                                                                                                                                                      <  1:430655/250‑1 (MQ=255)
              cGTCCAGTTCAGGCTTACGAACTTCGGCGATGTTGATCTGTGCAGGAACGCCAGCGATGTCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCg                                                                                                                                                                                                                                  <  1:283917/251‑1 (MQ=255)
                  cAGTTCAGGCTTACGAACTTCGGCGATGTTGATCTGTGCAGGAACGCCAGCGATGTCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGAGAACGATACGAGATACTGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGGCAGCGAAtt                                                                                                                                                                                                                              <  2:209571/251‑1 (MQ=255)
                       cAGGCTTACGAACTTCGGCGATGTTGATCTGTGCAGGACCGCCAGAGATGTCCGCTACGACCTTCCGCAGTTTTTCTACGTCATCACCTTTTTTACCGCTAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTATTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTg                                                                                                                                                                                                                         >  1:447738/1‑251 (MQ=255)
                         ggCTTACGAAATTCGGCGATGTTGATCTGTGCAGGAACGCCAGCGCTGTCCGCTACGTCCTTCCGCAGTTTTTATACGTCTTCACCTTTTTTACCGATAACGCTACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGCTACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCCGCGAATTCTTTGgt                                                                                                                                                                                                                       >  2:318149/1‑251 (MQ=255)
                                     tCGGCGATGTTGATCTGTGCAGGAACGCCAGCGATGTCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAAc                                                                                                                                                                                                             <  1:425972/249‑1 (MQ=255)
                                      cGGCGATGTTGATCTGTGCAGGAACGCCAGCGATGTCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTGAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTGCGCAAACCa                                                                                                                                                                                                           <  2:442092/250‑1 (MQ=255)
                                       agcgctGTTGATCTGTGCAGGAACGACGGCGATGTCCACTCCGACCTTACGAAGTGTTTCTACGTCTTCACCTTTTTTACCGATAACGATACGCGGCCGAGCAGTGTGAATGGTTACACGGATGTTCTTAGCCGGACGCTCGATCACGATACTAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGCACTGACGTACTTTAAAATCCCTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAgg                                                                                                                                                                                                         >  2:36712/6‑251 (MQ=255)
                                          gATGTTGATCTGTGCAGGAACGCCAGCGATGTCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGCTAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTCGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGt                                                                                                                                                                                                        >  1:45837/1‑249 (MQ=255)
                                          gATGTTGATCTGTGCAGGAACGCCAGCGATGTCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTag                                                                                                                                                                                                      <  1:82807/251‑1 (MQ=255)
                                               tGATCTGTGCAGGAACGCCAGCGATGTCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGtt                                                                                                                                                                                                  <  1:234928/250‑1 (MQ=255)
                                                     gtgCAGGAACGCCAGCGATGTCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATg                                                                                                                                                                                             <  2:383493/249‑1 (MQ=255)
                                                       gCAGGAACGCCAGCGATGTCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCTATAACGATACGAGATACGGACGCTTTAGCCAGTGCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGtt                                                                                                                                                                                          <  1:297787/250‑1 (MQ=255)
                                                        cAGGAACGCCAGCGATGTCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGtttt                                                                                                                                                                                        >  2:157911/1‑251 (MQ=255)
                                                        cAGGAACGCCAGCGATGTCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGtttt                                                                                                                                                                                        <  2:394275/251‑1 (MQ=255)
                                                                ccAGCGATGTCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATAcc                                                                                                                                                                                <  1:381388/251‑1 (MQ=255)
                                                                  aGCGATGTCCGCTACGACCTTACGCAGTTTTTCTACGTTTTCACCTTTTTTATCGATAACGATAACCGGGCGAGCAGTGTGCATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGCTACGCACGCTTTAGCCAGTTCCTTAGTCAGGTCCTGAAGTACTTTAAACTCACTGTCCAGGTTGTCAGCGAATTCTTTGGCGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATAccc                                                                                                                                                                               >  1:447538/1‑250 (MQ=255)
                                                                     gATGTCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGc                                                                                                                                                                           <  2:300144/251‑1 (MQ=255)
                                                                     gATGTCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGc                                                                                                                                                                           >  1:291824/1‑251 (MQ=255)
                                                                        gTCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGa                                                                                                                                                                         >  1:208987/1‑250 (MQ=255)
                                                                        gTCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGa                                                                                                                                                                         >  1:151375/1‑250 (MQ=255)
                                                                         tCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAAt                                                                                                                                                                       >  2:262854/1‑251 (MQ=255)
                                                                            gCTACGGCCTTACGCAGTTTTTCTACGTCTTCCCCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGGTTTACAATACCCAGGCGAATAcc                                                                                                                                                                    <  2:323390/251‑1 (MQ=255)
                                                                                 gACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATTCCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAg                                                                                                                                                               <  1:261841/251‑1 (MQ=255)
                                                                                 gACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAg                                                                                                                                                               <  1:82958/251‑1 (MQ=255)
                                                                                      tACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGGTCCTTAGTCAGGTACTGACGTGCTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGAt                                                                                                                                                            <  1:28091/249‑1 (MQ=255)
                                                                                       actcAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGt                                                                                                                                                          >  1:361600/4‑250 (MQ=255)
                                                                                       aagCAGTGTTTCTAAGTCTGCACCTGTGTTACCGAGAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGAGGATCTGAGCCGGACGCTCGATAACGAGACGGGAGACGGACGCTTTAGACAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAAGTCGTTGGTGTGCGCAAACCAGGTAGAGGGCCATGGTGTTACAATACCGAGGCGAAGGCCATTAGGATGt                                                                                                                                                          <  2:43720/248‑1 (MQ=255)
                                                                                       aCGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGt                                                                                                                                                          <  2:351566/250‑1 (MQ=255)
                                                                                          cAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACtt                                                                                                                                                      <  1:355579/251‑1 (MQ=255)
                                                                                                  cTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGt                                                                                                                                                                                                        <  2:228662/193‑1 (MQ=255)
                                                                                                  cTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGt                                                                                                                                                                                                        >  1:228662/1‑193 (MQ=255)
                                                                                                      gTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTg                                                                                                                                          <  1:262854/251‑1 (MQ=255)
                                                                                                          tCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTa                                                                                                                                       <  1:448034/250‑1 (MQ=255)
                                                                                                                  tttACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATGCCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCag                                                                                                                               <  1:84145/250‑1 (MQ=255)
                                                                                                                    tACCGATAACGATACCCGGGCGAGCAGTGTGACTGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATAAGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTCAAATCACTGTCAAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGt                                                                                                                            >  2:179361/1‑251 (MQ=255)
                                                                                                                     aCCGATAACGATAACGGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATACCGATCCGAGCTACGGACGCTTTAGACAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGCTGTCAGCGACTCCTTTGGTGTTCGCAATCCAGGTAGAGTTCCCTGGTTTTACAATACCCAGGCGAATACCATTAGGATGTAATTTCTGACCCATTGCTAGTCTcccgcgtc                                                                                                                           >  2:284132/1‑245 (MQ=255)
                                                                                                                        gATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGtctc                                                                                                                         >  1:139808/1‑250 (MQ=255)
                                                                                                                             cGATACCCGGGCGAGCAGTGTGAATGGTTACACGGAGGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGAt                                                                                                                   <  2:390873/251‑1 (MQ=255)
                                                                                                                                tACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAGATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGCTCGCAAACCAGGTAGAGATCCATGGTTTTACAATACCCAGGCGAATACCAGTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATccg                                                                                                                >  2:347336/1‑249 (MQ=255)
                                                                                                                                     gggCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGacac                                                                                                            >  2:142930/1‑250 (MQ=255)
                                                                                                                                           gCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACCCAACcac                                                                                                      >  2:26499/1‑250 (MQ=255)
                                                                                                                                             aGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGGTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGt                                                                                                   <  2:45837/251‑1 (MQ=255)
                                                                                                                                                      ggTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCt                                                                                          <  1:459855/251‑1 (MQ=255)
                                                                                                                                                      ggTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCt                                                                                          >  1:153897/1‑251 (MQ=255)
                                                                                                                                                       gTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTg                                                                                         <  1:236112/251‑1 (MQ=255)
                                                                                                                                                                aTGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTc                                                                                >  2:262107/1‑251 (MQ=255)
                                                                                                                                                                  gCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAg                                                                              <  2:301779/251‑1 (MQ=255)
                                                                                                                                                                      ttAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAgggtg                                                                          >  2:44579/1‑248 (MQ=255)
                                                                                                                                                                        aGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGc                                                                         <  2:430627/250‑1 (MQ=255)
                                                                                                                                                                            ggACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCt                                                                    >  1:181726/1‑251 (MQ=255)
                                                                                                                                                                            ggACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCt                                                                    >  1:100253/1‑251 (MQ=255)
                                                                                                                                                                                gCTCGATAACGAGACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTTTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTGCAGTATGCGATCTGCAc                                                                <  2:425840/251‑1 (MQ=255)
                                                                                                                                                                                    gATAAAGATACGAGCTACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGCCACAACCCCAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGAcc                                                            >  2:137290/1‑251 (MQ=255)
                                                                                                                                                                                       aaCGATACGCGATACGGACGCTTTAGCCAGTTCCTTAGACAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGAAACCCAGGTAGAGTTCCATGGTTTTACAATACACAGGCGAATACCATTAGGCTGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCCCTTCAGGATGCGATCTGCACGAACttt                                                         >  1:420470/1‑251 (MQ=255)
                                                                                                                                                                                        aCGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCtttt                                                        <  1:465102/251‑1 (MQ=255)
                                                                                                                                                                                          gATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGc                                                      >  1:148283/1‑251 (MQ=255)
                                                                                                                                                                                              cGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCAcgc                                                  >  1:71506/1‑251 (MQ=255)
                                                                                                                                                                                               gagaTACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCAcgcg                                                 <  2:186143/251‑1 (MQ=255)
                                                                                                                                                                                                    aCGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCAAGCGGCAt                                             >  1:326459/1‑250 (MQ=255)
                                                                                                                                                                                                      ggACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAAt                                          <  1:458341/251‑1 (MQ=255)
                                                                                                                                                                                                               aGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCAt                                 <  1:157911/251‑1 (MQ=255)
                                                                                                                                                                                                                           aGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTCGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGGGGCTGGTGCGACTCAGGATGCGATCTGCCCGACCTTTTGCCCGCGGCCTAATTCGCTTCATGCCCGGGccctc                     >  2:287664/1‑248 (MQ=255)
                                                                                                                                                                                                                            gTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTc                                                                                >  2:343795/1‑191 (MQ=255)
                                                                                                                                                                                                                            gTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTc                                                                                <  1:343795/191‑1 (MQ=255)
                                                                                                                                                                                                                             tCAGGTACTGACGTACTTTAAAATCCCTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCCGGATGCGATCTGCACGACCTTTTGCACGCGGGATAATGCGCTTCCTGCTCGGGCCTTCGt                   >  2:439449/1‑251 (MQ=255)
                                                                                                                                                                                                                                       aCGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCg                                                                        <  1:88237/188‑1 (MQ=255)
                                                                                                                                                                                                                                       aCGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCg                                                                        >  2:88237/1‑188 (MQ=255)
                                                                                                                                                                                                                                               tAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTa  <  2:153897/250‑1 (MQ=255)
                                                                                                                                                                                                                                                aaaaTCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCAt                                             <  1:36080/206‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  
AGCAACCAGTTTTGCGTCCAGTTCAGGCTTACGAACTTCGGCGATGTTGATCTGTGCAGGAACGCCAGCGATGTCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCGCTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTA  >  NC_000913/3449533‑3450021

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