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,8830CT98.2% 173.2 / ‑6.1 55A2V (GCT→GTT) rplA50S ribosomal subunit protein L1
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (25/28);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

TTCTTTCACTTTCGTTACCAAGACCCCGCCGGCAGCAGTTCTGCTGAAAAAAGCGGCTGGTATCAAGTCTGGTTCCGGTAAGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGCGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGC  >  NC_000913/4178641‑4179123
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                
ttCTTTCACTTTCGTTACCAAGACCCCGCCGGCAGCAGTTCTGCTGAAAACAGCGGCTGGTATCAAGTCTGGTTCCGGTACGCCGAACAAAGCCAAAGTGGGAAAAATTTCCCGCGCTCAGCTGCAAGAAATCTCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACTTTCCATGGGCCTGGTAGTGGAGGACTGACAAATTGATAAACTg                                                                                                                                                                                                                                           >  2:479475/1‑250 (MQ=255)
 tCTTTCACTTTCGTTACCAAGACCCCGCCGGCAGCAGTTCTGCTGAAAAAAGCGGCTGGTATCAAGTCTGGTTCCGGTAAGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGa                                                                                                                                                                                                                                          <  1:418307/250‑1 (MQ=255)
      cACTTTCGTTACCAAGACCCCGCCGGCAGCAGTTCTGCTGAAAAAAGCGGCTGGTATCAAGTCTGGTTCCGGTAAGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAgc                                                                                                                                                                                                                                    <  1:423845/251‑1 (MQ=255)
      cACTGTCGTTACCAAGACCCCGCCGGCAGCAGTTCTGCTGAAAAAAGCGGCTGGTATCAAGTCTGGTTCCGGTAAGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAgc                                                                                                                                                                                                                                    <  2:532002/251‑1 (MQ=255)
       acattCGTTTCCAAGACCCCTCCGGCAGCAGTTCTGCTGAAAAAAGCGGCTGGTATCAAGTCTGGTTCCGGTAAGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGTAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATTACTCGCTCCATCGAATGTACTGCACGTTCCATGGGCCGGGTAGTGTAGGACTAAGAAATAGTTAAACGGGCCAAgcg                                                                                                                                                                                                                                   <  2:531335/248‑1 (MQ=255)
                                 agcagTTCTGCTGAAAAAAGCGGCTGGTATCAAGTCTGGTTCCGGTAAGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTg                                                                                                                                                                                                         <  2:178636/251‑1 (MQ=255)
                                 agcagTTCTGCTGAAAAAAGCGGCTGGTATCAAGTCTGGTTCCGGTAAGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACGCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTg                                                                                                                                                                                                         <  2:296501/251‑1 (MQ=255)
                                              aaaaaaGCGGCTGGTATCAAGTCTGGTTCCGGTAAGCCGAACAAAGACAAAGGGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGATGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTACACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTATTATCCTCGAGAAAGGTGATGCAACCAAAc                                                                                                                                                                                             >  1:582838/1‑250 (MQ=255)
                                                   aGCGGCTGGTATCAAGTCTGGTTCCGGTAAGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGGCCAAGCGCATGCGTGTTAGCCGCGAGAAAGTTGATGCAACCAAACAGTACg                                                                                                                                                                                       <  1:195659/251‑1 (MQ=255)
                                                      ggCTGGTATCAAGTCTGGTTCCGGTAAGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGTACTGGTTAAACTGACCAAGCGCCTGCGTGTTATCCGCGAGAAAGTTGATGCAACCCAACCGTACGAc                                                                                                                                                                                     >  2:484211/1‑250 (MQ=255)
                                                       gCTGGTATCAAGTCTGGTTCCGGTAAGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACAt                                                                                                                                                                                   <  2:183435/251‑1 (MQ=255)
                                                       gCTGGTATCAAGTCTGGTTCCGGTAAGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACAt                                                                                                                                                                                   >  1:435602/1‑251 (MQ=255)
                                                                   tCTGGTTCCGGTAAGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTAt                                                                                                                                                                       >  1:176676/1‑251 (MQ=255)
                                                                      ggTTCCGGTAAGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCg                                                                                                                                                                     <  1:566144/250‑1 (MQ=255)
                                                                           cGGTAAGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCActgc                                                                                                                                                               >  1:202803/1‑251 (MQ=255)
                                                                             gTAAGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCActgctg                                                                                                                                                             >  1:108479/1‑251 (MQ=255)
                                                                                aGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGCTGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACAAAGCTATCGCACTGCTGaa                                                                                                                                                           >  1:49960/1‑250 (MQ=255)
                                                                                             caaagTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTg                                                                                                                                             <  2:335316/251‑1 (MQ=255)
                                                                                               aagTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCt                                                                                                                                           >  1:324101/1‑251 (MQ=255)
                                                                                                    ggTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAAtt                                                                                                                                      <  1:243783/251‑1 (MQ=255)
                                                                                                     gTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTc                                                                                                                                     >  1:311420/1‑251 (MQ=255)
                                                                                                                cgcgCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAgtgt                                                                                                                          <  2:165756/251‑4 (MQ=255)
                                                                                                                 gcgcTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGTCCAAGCGCATGCTTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGtgtg                                                                                                                         <  1:369296/251‑1 (MQ=255)
                                                                                                                  cgcTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTAGTAGAAAGTGTgg                                                                                                                        >  2:22717/1‑251 (MQ=255)
                                                                                                                           gCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCt                                                                                                                <  2:221285/250‑1 (MQ=255)
                                                                                                                                 aaTCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAAcc                                                                                                         >  1:270062/1‑251 (MQ=255)
                                                                                                                                          gACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCg                                                                                                <  1:146451/251‑1 (MQ=255)
                                                                                                                                                           tGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGAc                                                                               >  1:301364/1‑251 (MQ=255)
                                                                                                                                                            gACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGAc                                                                               <  2:171539/250‑1 (MQ=255)
                                                                                                                                                                    ccGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGCCTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCActgct                                                                                                                                                              >  2:477501/1‑163 (MQ=255)
                                                                                                                                                                    ccGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCActgct                                                                                                                                                              <  1:477501/163‑1 (MQ=255)
                                                                                                                                                                     cGACATTGAAGCGATGACTCGCTCCATCGAAGGTACAGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGACAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAacgtac                                                                     >  2:421663/1‑251 (MQ=255)
                                                                                                                                                                       aCATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAacgtacgt                                                                   >  1:432793/1‑251 (MQ=255)
                                                                                                                                                                       aCATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAacgtacgt                                                                   >  2:324876/1‑251 (MQ=255)
                                                                                                                                                                         aTTGAAGCGATGACTCGCTCCAGCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCAGCGACGCTCGTAAATCTGACCAGAacgtacgtag                                                                 <  2:482948/251‑3 (MQ=255)
                                                                                                                                                                          ttGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACgtggg                                                                >  1:96033/1‑250 (MQ=255)
                                                                                                                                                                                  aTGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGAGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGt                                                        >  1:23655/1‑251 (MQ=255)
                                                                                                                                                                                     aCTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACt                                                     <  2:441237/251‑1 (MQ=255)
                                                                                                                                                                                     aCTCGCTCCATCGAAGGTACTGCACGGTCCATTGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCTGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGTTTCAACTGTACt                                                     <  1:479475/251‑1 (MQ=255)
                                                                                                                                                                                            ccATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGc                                                <  1:283166/249‑1 (MQ=255)
                                                                                                                                                                                             cATCGAATGTACTGCACGTTCCATGGGCCTTGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGAGGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGCCTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGTCCATAACGTACGTGGTGCAACTGTACTGCCGctcg                                             <  1:457417/251‑4 (MQ=255)
                                                                                                                                                                                             cATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACg                                             >  1:441174/1‑251 (MQ=255)
                                                                                                                                                                                             cATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGGGGTGCAACTGTACTGCCGCACg                                             <  1:169804/251‑1 (MQ=255)
                                                                                                                                                                                                cGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGt                                           >  1:550904/1‑250 (MQ=255)
                                                                                                                                                                                                   aGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCACCTGTACTGCCGCACGGTACTg                                       >  2:106279/1‑251 (MQ=255)
                                                                                                                                                                                                    ggTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTgg                                      <  2:437043/251‑1 (MQ=255)
                                                                                                                                                                                                           cACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGccgttcc                               <  2:108479/251‑1 (MQ=255)
                                                                                                                                                                                                            aCGTTCGATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACTACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGccgttccg                              <  2:530013/251‑1 (MQ=255)
                                                                                                                                                                                                                tCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACTTGGTGCAACTGTACTGCCGCACGGTACTGGCcgttccgttc                           <  2:486379/250‑1 (MQ=255)
                                                                                                                                                                                                                  cATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGAGTGTTATCCGCGAGAAAGATGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGTCGAATGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTAATGCCGCACGGTACTGGCCGGTCCGTTcgcg                        >  2:454647/1‑251 (MQ=255)
                                                                                                                                                                                                                    tGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGt                       >  1:215436/1‑250 (MQ=255)
                                                                                                                                                                                                                              aGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACCTCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTAttt             >  1:216776/1‑250 (MQ=255)
                                                                                                                                                                                                                                     gaCTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAgg     <  2:432793/251‑1 (MQ=255)
                                                                                                                                                                                                                                       cTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCACGGTg   >  2:541518/1‑251 (MQ=255)
                                                                                                                                                                                                                                         aaGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGc  >  1:548956/1‑250 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                
TTCTTTCACTTTCGTTACCAAGACCCCGCCGGCAGCAGTTCTGCTGAAAAAAGCGGCTGGTATCAAGTCTGGTTCCGGTAAGCCGAACAAAGACAAAGTGGGTAAAATTTCCCGCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGCGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGC  >  NC_000913/4178641‑4179123

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: