Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,178,883 C→T A2V (GCT→GTT)  rplA → 50S ribosomal subunit protein L1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,178,8830CT97.9% 154.5 / NA 48A2V (GCT→GTT) rplA50S ribosomal subunit protein L1
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (26/21);  minor base G (0/1);  total (26/22)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

ACCAAGACCCCGCCGGCAGCAGTTCTGCTGAAAAAAGCGGCTGGTATCAAGTCTGGTTCCGGTAAGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGCGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAA  >  NC_000913/4178657‑4179126
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   
aCCAAGACCCCGCCGGCAGCAGTTCTGCTGAAAAAAGCGGCTGGTATCAAGTCTGGTTCCGGTCAGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCgtgt                                                                                                                                                                                                                             >  2:32951/1‑251 (MQ=255)
aCCAAGACCCCGCCGGCAGCAGTTCTGCTGAAAAAAGCGGCTGGTATCAAGTCTGGTTCCGGTAAGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCgtgt                                                                                                                                                                                                                             <  1:270697/251‑1 (MQ=255)
   aaGACCCCGCCGGCAGCAGTTCTGCTGAAAAAAGCGGCTGGTATCAAGTCTGGTTCCGGTAAGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTAt                                                                                                                                                                                                                          <  2:171453/251‑1 (MQ=255)
      aCCCCGGCGGCAGCAGTACTGCTGAAAAAAGCGGCTGGTATCAAGTCTGGTTCCGGTAAGCCGAACAAAGACAAAGTGGGTAAAATTGCCCGCGCTCAGCTGCAGGAAACCGCGCAGACCAAAGCTGCCGACATTACTGGTGCCGACATTGAAGCGATGGCGCGCTCCAGCGAAGTTACTGCACGTTCCATGGGCCTGGTGGTGGAGGACTAAGAAATGGGTAAACTGACCAATCGCATGCGTGTTATCcg                                                                                                                                                                                                                       <  2:204165/251‑1 (MQ=255)
        cccGCCGGCAGCAGTTCTGCTGAAAAAAGCGGCTGGTATCAAGTCTGGTTCCGGTAAGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCcgcg                                                                                                                                                                                                                     <  1:94301/251‑1 (MQ=255)
               gcagcaGTTCTGCTGAAAAAAGCGGCTGGTATCAAGTCTGGTTCCGGTAAGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACGGGTGCCGACATTGAAGCGATGACTCGCTCCATCTAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATTCGTGTTATCCGCGATAAAGt                                                                                                                                                                                                              <  2:163711/251‑1 (MQ=255)
                  gcagTTCTGCTGAAAAAAGCGGCTGGTATCAAGTCTGGTTCCGGTAAGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTg                                                                                                                                                                                                            >  1:61808/1‑250 (MQ=255)
                              aaaaaaGCGGCTGGTATCAAGTCTGGTTCCGGTAAGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAAc                                                                                                                                                                                                >  2:48282/1‑250 (MQ=255)
                                 aaaGCGGCTGGTATCAAGTCTGGTTCCGGTAAGCCGAACAAAGACAAAGTGGGTAACATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCATGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGt                                                                                                                                                                                             >  1:210950/1‑250 (MQ=255)
                                  aaGCGGCTGGTATCAAGTCTGGTTCCGGTAAGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTAc                                                                                                                                                                                           >  1:52613/1‑251 (MQ=255)
                                   aGCGGCTGGTATCAAGTCTGGTTCCGGTAAGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACg                                                                                                                                                                                          <  2:121818/251‑1 (MQ=255)
                                      ggCTGGTATCAAGTCTGGTTCCGGTAAGCCGAACAAAGACAAAGTGGGTAACATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGAcc                                                                                                                                                                                       >  1:129557/1‑250 (MQ=255)
                                         tGGTATCAAGTCTGGTTCCGGAAAGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAAATGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATc                                                                                                                                                                                     >  2:25215/1‑250 (MQ=255)
                                          ggTATCAAGTCTGGTTCCGGTAAGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAAACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTAACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCGTCCAACCAGTACTACATCa                                                                                                                                                                                    <  2:192133/250‑1 (MQ=255)
                                                                        aaagacaaagTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGc                                                                                                                                                     <  2:182781/251‑1 (MQ=255)
                                                                              aaagTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTaaa                                                                                                                                                                                                                                                 <  1:219633/153‑1 (MQ=255)
                                                                              aaagTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTaaa                                                                                                                                                                                                                                                 >  2:219633/1‑153 (MQ=255)
                                                                               aagTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATc                                                                                                                                                                                     >  1:4274/1‑212 (MQ=255)
                                                                               aagTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATc                                                                                                                                                                                     <  2:4274/212‑1 (MQ=255)
                                                                                     gaaaaaTTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTaa                                                                                                                                            >  1:31065/3‑247 (MQ=255)
                                                                                     gTAAAATTTCCCTCGCTCAGCTGCATGAAATCGCGCAGACCAAAGCTGCCGACATGACTGCTGCCGACATTGAAGCTATGACTCTCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTaa                                                                                                                                            <  2:31065/247‑1 (MQ=255)
                                                                                               ccGCGCTCAGCGACAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGTCATTGAATCGATGACTCGCTCCATCGAAGGTACTGCACGGTCCATGGGCCTGGTATTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATTCAACCAAACAGTACGACATCAACGAAGCTAGCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAgtg                                                                                                                              <  1:184418/251‑3 (MQ=255)
                                                                                                   gcTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTgg                                                                                                                           >  1:54843/1‑250 (MQ=255)
                                                                                                   gcTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTgg                                                                                                                           >  1:57361/1‑250 (MQ=255)
                                                                                                   gcTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGa                                                                                                                          <  1:20406/251‑1 (MQ=255)
                                                                                                        gCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAg                                                                                                                     <  2:176543/251‑1 (MQ=255)
                                                                                                          tGCAGGCAATCGCCCAGACCAAAGCTGCCGACATGACTGATGCCGACATTAAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCt                                                                                                                   >  2:125419/1‑251 (MQ=255)
                                                                                                           gCAGGAAAGCGCGCAGACCAAAGCTGCCGACATGACGGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACAGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTTGCt                                                                                                                   <  2:210950/250‑1 (MQ=255)
                                                                                                            cAGGAAATCGCGCAGACCAAAGCTGCCGACATGCCTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGt                                                                                                                 >  2:148894/1‑251 (MQ=255)
                                                                                                                                    gCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTaa                                                                                         <  1:231599/251‑1 (MQ=255)
                                                                                                                                       gACATGACTGGTGACGACATTGAAGAGATGACTCGCAAAATCGAAGGTAATGCACGTTCCATGGGCCTGGAAGTGGAGGACTAAGAAATGGTTAAAATGAACAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGAATGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGAATCGACGCTAGTAAAt                                                                                       >  1:286071/1‑250 (MQ=255)
                                                                                                                                         cATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGACAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTg                                                                                    >  1:17803/1‑251 (MQ=255)
                                                                                                                                                           tGAAGCGATGACTCGCTCCATCGCAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCCACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGCGTGGACGTAGCTGTTAACATCGGCATCGACGCTCGTAAATCTGACCAGAACGCACgcggtg                                                                  >  2:3498/1‑251 (MQ=255)
                                                                                                                                                              aGCGATGACTCGCTCCATCGAAGGTACTGCACGTTACATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCa                                                                >  2:189403/1‑250 (MQ=255)
                                                                                                                                                                  aTGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGt                                                           >  1:174665/1‑251 (MQ=255)
                                                                                                                                                                         gCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCg                                                                                                   >  2:19241/1‑204 (MQ=255)
                                                                                                                                                                         gCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTGATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCg                                                                                                   <  1:19241/204‑1 (MQ=255)
                                                                                                                                                                             cATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACg                                                <  2:156278/251‑1 (MQ=255)
                                                                                                                                                                              atcgacgcgaCTGCACGTTCCATGGCCCTGGGCGTGGAGGACTAAGACATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTCGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACgg                                               >  2:21921/10‑251 (MQ=255)
                                                                                                                                                                                     gTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGc                                        <  1:216037/251‑1 (MQ=255)
                                                                                                                                                                                      tACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGcc                                       >  2:114126/1‑251 (MQ=255)
                                                                                                                                                                                         tGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGccgtt                                    >  1:38143/1‑251 (MQ=255)
                                                                                                                                                                                              gTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGccgttccgtt                               <  2:1959/251‑1 (MQ=255)
                                                                                                                                                                                              gTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGAACTGGccgttccgtt                               >  1:233172/1‑251 (MQ=255)
                                                                                                                                                                                                                   aggaCTCAGAAATGGTTAAACTGACCACGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCa           >  1:73357/1‑250 (MQ=255)
                                                                                                                                                                                                                   aggaCTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCaa          <  1:237840/251‑1 (MQ=255)
                                                                                                                                                                                                                    ggaCTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCa           >  1:203700/1‑249 (MQ=255)
                                                                                                                                                                                                                           aaaaTGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCaaa  <  1:162585/250‑1 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   
ACCAAGACCCCGCCGGCAGCAGTTCTGCTGAAAAAAGCGGCTGGTATCAAGTCTGGTTCCGGTAAGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGCGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAA  >  NC_000913/4178657‑4179126

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