Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,202,009 A→T 6.7% Q40L (CAG→CTG)  cca → fused tRNA nucleotidyl transferase/2'3'‑cyclic phosphodiesterase/2'nucleotidase and phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,202,0090AT 6.7% 254.2 / 3.1 93Q40L (CAG→CTG) ccafused tRNA nucleotidyl transferase/2'3'‑cyclic phosphodiesterase/2'nucleotidase and phosphatase
Reads supporting (aligned to +/- strand):  ref base A (37/47);  new base T (2/5);  total (40/53)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGTATTTGACGAAGGGATGATGGCGTGAAGATTTATCTGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACCGGTTTTACTTGCTATGCCGCACCGGATGTCACGCTGGAAGATGATCTTAAGC  >  NC_000913/3201866‑3202158
                                                                                                                                               |                                                                                                                                                     
gtgtATTTGACGAAGGGATGATGGCGTGAAGATTTATCTGGTCGGTGGTGATGTTCGGTATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGAGGGCTACCTGCAGGt                                                                                                                                                 >  1:691132/1‑150 (MQ=255)
ggggaTTTGACGACGGGATGATGGCGTGAAGATTTATCTGGTCGGTGGTGCTGTTCGGGATTCATTGTTAGGGCTACCGGTCAAAGCCAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGt                                                                                                                                                 <  2:1632831/146‑1 (MQ=255)
 tgtATTTGACGAAGGGCTGATGGCGTGAAGAGTTATCTGCTCGCTGGTGCTGCTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTa                                                                                                                                                <  1:100885/150‑1 (MQ=255)
 tgtATTTGACGAAGGGATGATGGCGTGAAGATTTATCTGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTa                                                                                                                                                >  2:1194765/1‑150 (MQ=255)
    aTTTGACGAAGGGATGATGGCGTGAAGATTTAACTGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACTCCACAGGAGATGCTCGACTCGGGCTACCAGCAGGTAGGc                                                                                                                                             >  1:1141110/1‑150 (MQ=255)
      ttGACGACGGGAGGATGGCGTGAAGATTTATCTGGTCGGTGGGGCTGTGCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGCCTACCTGCTGGTAGCCcg                                                                                                                                           <  1:1076287/150‑1 (MQ=255)
       tGACGAAGGGATGATGGCGTGAAGATTTATCTGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGGGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCcgc                                                                                                                                          >  1:1521853/1‑150 (MQ=255)
       tGACGAAGGGATGATGGCGTGAAGATTTATCTGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACTCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCcgc                                                                                                                                          >  1:1819709/1‑150 (MQ=255)
       tGACGAAGGGATGAGGGCGTGAAGATTTATCTGGTCGGTGGTGCGGTTCGGGATGCATTGTGATGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGCCTACCTGCTGGTAGCCcgc                                                                                                                                          <  1:154756/150‑1 (MQ=255)
       tGACGAAGGGAGGATGGCGTGAAGATTTATCAGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGTGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGCCTACCTGCTGGTAGCCcgc                                                                                                                                          <  2:1233368/150‑1 (MQ=255)
        gACGAAGGGATGATGGCGTGAAGATTTATCTGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAATGAAAGATATGGGGTGTTGTTCGGCATTACGCAAAAGAAGATCAACGACCCGGCTTACATGCAGTTAGCCcgcg                                                                                                                                         >  1:214311/1‑150 (MQ=255)
          cGAAGGGATGATGGCGTGAAGATTTATCTGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGAt                                                                                                                                       <  1:528862/150‑1 (MQ=255)
          cGAAGGGATGAGGGCGTGAAGATTTATCTGGTCGGTGGCGCTGTTCGGGATGCATTGTTAGGGCTACCGGGCAAAGACAGAGATTGGGTGGTGGGCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGAt                                                                                                                                       <  1:365862/150‑1 (MQ=255)
            aaGGGAGGATGGCGTGAAGATTTATCTGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGAttt                                                                                                                                     <  2:811241/150‑1 (MQ=255)
              gggATGATGGCGTGAAGATTTATCTGGTCGGTGGGGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGCCTACCTGCTGGTAGCCcgcgtcgagc                                                                                                                                   <  1:1285924/150‑7 (MQ=255)
               gtccgatgGAGTGCAGATTTATCTGGTCGGTGGTGCTCTACGGCATGCAGTGTATGGTCTACCGGTCCAAGACAGAGATTGGGGGGTGGTCGGCAGTACGCCACAGTAGATGCTAGACGCGGGCTACCAGCAGGTAGGCCGCGATGTTcc                                                                                                                                  <  1:338845/146‑1 (MQ=255)
                   gatgGCGTGAAGATTTATCTGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGGCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCtgtg                                                                                                                              <  1:1020994/150‑1 (MQ=255)
                     tgGCGGGAAGATTTATCTGGTCGGTGGGGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGtt                                                                                                                            <  2:372290/150‑1 (MQ=255)
                          tGAAGATTTATCTGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGc                                                                                                                       >  1:1570351/1‑150 (MQ=255)
                              gATTTATCTGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATcc                                                                                                                   >  1:795869/1‑150 (MQ=255)
                              gATTTATCTGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGGTTCTGCATcc                                                                                                                   >  2:64163/1‑150 (MQ=255)
                                   agcTGGTCGGTGCTGCTGTGCGGGATGCCTTGTCAGGGCGACCGGTCAAAGACAGAGATTGGGGGGGGGGCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCaaa                                                                                                              <  2:249500/148‑1 (MQ=255)
                                   aTCTGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCaaa                                                                                                              <  2:595639/150‑1 (MQ=255)
                                     cTGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGAATGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACg                                                                                                            <  2:1441212/150‑1 (MQ=255)
                                      tGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAACACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGc                                                                                                           >  2:175453/1‑150 (MQ=255)
                                       ggTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCa                                                                                                          >  1:181076/1‑150 (MQ=255)
                                       ggTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCa                                                                                                          >  1:581700/1‑150 (MQ=255)
                                       ggTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGGCGGCAGTACGCCACTGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCa                                                                                                          <  1:516343/150‑1 (MQ=255)
                                         ggggtgggGCTGATCGGGATGCTTTGTTAGGGCGACAGGGGAAAGACTGTGAGTGGGGGGTGGTCGGCAGTACGCCACAGGAGCTGCTCGACGCGGGCTACCAGCATGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATg                                                                                                        <  2:97548/148‑1 (MQ=255)
                                          cGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGa                                                                                                       >  1:631699/1‑150 (MQ=255)
                                              tgtgCTGTTCGGGATGGATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTTCCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAagag                                                                                                   <  1:1737362/149‑1 (MQ=255)
                                              ggtgCTGTTCGGGATGCTTTGTCAGGGCTACCGGACAAACACCGCGATGGGGGGGTGGTCGGCTGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCCGTGTTTCTGCATCCGCAAACGCATGAagag                                                                                                   <  2:1653398/150‑1 (MQ=255)
                                              ggtgCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGGGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTTTTTCTGCATCCGCAAACGCATGAagag                                                                                                   <  1:563923/150‑1 (MQ=255)
                                               gtgCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGt                                                                                                  >  2:1547832/1‑150 (MQ=255)
                                               gtgCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGAGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGt                                                                                                  >  2:962368/1‑150 (MQ=255)
                                                ggCTGTTCGGGCTGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTa                                                                                                 <  1:408557/149‑1 (MQ=255)
                                                  cTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACTCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATg                                                                                               <  1:467683/150‑1 (MQ=255)
                                                    gTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATgcg                                                                                             >  2:333429/1‑150 (MQ=255)
                                                       cGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGATGAGTATGCGCTg                                                                                          >  2:912091/1‑150 (MQ=255)
                                                       cGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTg                                                                                          >  1:851025/1‑150 (MQ=255)
                                                           acgCATTGTTAGGTCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGGATGTATGCCGCGATTTTCCTGTGTTTCTGCATCCTCAAACTCATGAAGAGTATGCGCTGGCAc                                                                                      <  2:647894/148‑1 (MQ=255)
                                                              cATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTa                                                                                   >  2:678730/1‑150 (MQ=255)
                                                               aTTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTAc                                                                                  >  1:1535398/1‑150 (MQ=255)
                                                               aTTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGAGCAGGCACGTAc                                                                                  >  1:1535405/1‑150 (MQ=255)
                                                                    tAGGGCTTCCGGTCAAAGACAGAGATTGGGTGGTGGTCGTCAGTACGCCACAGGAGTTGCTCGTCGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAAc                                                                             <  1:1262631/150‑1 (MQ=255)
                                                                    tAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAAc                                                                             >  1:541064/1‑150 (MQ=255)
                                                                     aGGGCTACCGGTCAAAGACAGAGATTGGGTGGGGGGCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACg                                                                            <  2:1060938/150‑1 (MQ=255)
                                                                        gCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGaa                                                                         >  1:631383/1‑150 (MQ=255)
                                                                         cTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGaaa                                                                        >  1:811996/1‑150 (MQ=255)
                                                                          tACCGGTCAACGACAGAGATGGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAt                                                                       <  1:1792849/150‑1 (MQ=255)
                                                                           aCCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCTGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACGGAACGGAAAtc                                                                      <  1:318180/150‑1 (MQ=255)
                                                                              ggTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccgg                                                                   >  1:1034116/1‑150 (MQ=255)
                                                                              ggTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccgg                                                                   >  1:1377498/1‑150 (MQ=255)
                                                                               ggcAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggt                                                                  <  2:670104/148‑1 (MQ=255)
                                                                                tCAAAGCCAGAGATTGGGGGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggtt                                                                 <  1:1738080/150‑1 (MQ=255)
                                                                                tCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggtt                                                                 >  2:1640764/1‑150 (MQ=255)
                                                                                 cAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggttc                                                                >  2:510248/1‑150 (MQ=255)
                                                                                  aaaGCCAGAGATTGGGGGGTGGGCGGCAGATCCGCACAGGAGACGCCCGACGCGGGCTCCCTGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTCTGCGCTGGCACGTACCGAACGGAAAtccggttcc                                                               <  1:1324386/150‑1 (MQ=255)
                                                                                  aaaGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggttcc                                                               <  2:214311/150‑1 (MQ=255)
                                                                                    aGACAGTGTTTGTGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggttccgg                                                             <  2:1312381/150‑1 (MQ=255)
                                                                                    aGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggttccgg                                                             <  1:1325803/150‑1 (MQ=255)
                                                                                       cAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTa                                                          >  2:1210135/1‑150 (MQ=255)
                                                                                       cAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCAAGTACCGAACGGAAATCCGGTTCCGGTTa                                                          >  2:1552556/1‑150 (MQ=255)
                                                                                        ttagaTTGGGTGGTGGTCGGCAGTACGCCTCAGTCGATGCTCTACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCTAACGGAAATCCGGTTCCGGTTac                                                         <  2:857668/148‑1 (MQ=255)
                                                                                           gaTTGGGTGGTGGTCGGCAGTACGCCACAGGAGCTGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGATTCTGCATCCGCAAACGCAGGAAGAGTATGCGCTGGCACGTGCCGAACGGTAATACGGTTCCGGTTTCAcc                                                      >  2:534576/1‑150 (MQ=255)
                                                                                              tGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACCGGt                                                   >  2:1452417/1‑150 (MQ=255)
                                                                                              tGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCGGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACCGGt                                                   <  2:309039/150‑1 (MQ=255)
                                                                                               gggTGGTGGTCGGCAGTACGCCACTGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAATCGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACCGGtt                                                  <  2:1346678/150‑1 (MQ=255)
                                                                                                   ggtggtCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACCGGTTTTAc                                              >  1:1431158/1‑150 (MQ=255)
                                                                                                   ggtggtCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACCGGTTTTAc                                              <  2:1819709/150‑1 (MQ=255)
                                                                                                     tggtCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACCGGTTTTACtt                                            >  1:655084/1‑150 (MQ=255)
                                                                                                         cGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACCGGTTTTACTTGCTa                                        <  2:830051/150‑1 (MQ=255)
                                                                                                         cGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACCGGTTTTACTTGCTa                                        <  2:1671192/150‑1 (MQ=255)
                                                                                                          ggCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACCGGTTTTACTTGCTAt                                       >  2:1738862/1‑150 (MQ=255)
                                                                                                                   ccACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACCGGTTTTACTTGCTATGCCGCACCg                              <  1:149420/150‑1 (MQ=255)
                                                                                                                   ccACAGGAGAGGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATTCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACCGGTTTTACTTGCTATGCCGCACCg                              <  2:1553330/150‑1 (MQ=255)
                                                                                                                    cacaGGAGATGCTCGACGCGGGCTACCCGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGTCCGGAAATCCGGTTCCGGTTACACCGGTTTTACTTGCTATGCCGCACCgg                             >  1:312837/1‑150 (MQ=255)
                                                                                                                       aGGAGATGCGCGACGCTGCCTACCCGCCGGGAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACCGGTTTTACTTGCTATGCCGCACCGTATg                          <  2:901430/150‑1 (MQ=255)
                                                                                                                           gaTGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACCGGTTTTACTTGCTATGCCGCACCGGATGTCAc                      <  1:1181300/150‑1 (MQ=255)
                                                                                                                           gaTGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACCGGTTTTACTTGCTATGCCGCACCGGATGTCAc                      <  1:903514/150‑1 (MQ=255)
                                                                                                                               cTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACCGGTTTTACTTGCTATGCCGCACCGGATGTCACGCTg                  >  2:278779/1‑150 (MQ=255)
                                                                                                                               cTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACCGGTTTTACTTGCTATGCCGCACCGGATGTCACGCTg                  >  2:1390813/1‑150 (MQ=255)
                                                                                                                                 cGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACCGGTTTTACTTGCTATGCCGCACCGGATGTCACGCTGGa                <  1:1243227/150‑1 (MQ=255)
                                                                                                                                 cGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACCGGTTTTACTTGCTATGCCGCACCGGATGTCACGCTGGa                <  1:678730/150‑1 (MQ=255)
                                                                                                                                  gACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTGCCGAACGGAAATCCGGTTCCGGTTACACCGGTTTTACTTGCTATGCCGCACCGGATGTCACGCTGGaa               >  2:242992/1‑150 (MQ=255)
                                                                                                                                  gACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACCGGTTTTACTTGCTATGCCGCACCGGATGTCACGCTGGaa               <  1:215517/150‑1 (MQ=255)
                                                                                                                                    cgcgGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCTCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCAGTTTCCGGTTACACCGGTTTTACTTGCTATGCCGCACCGGATGTCACGCTGGAAga             <  1:373008/150‑1 (MQ=255)
                                                                                                                                    cgcgGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACCGGTTTTACTTGCTATGCCGCACCGGATGTCACGCTGGAAga             <  1:525635/150‑1 (MQ=255)
                                                                                                                                      cgGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGTACGGAAATCCGGTTCCGGTTACACCGGTTTTACTTGCTATGCCGCACCGGATGTCACGCTGGAAgatg           <  1:1632009/150‑1 (MQ=255)
                                                                                                                                          cTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACCGGTTTTACTTGCTATGCCGCACCGGATGTCACGCTGGAAGATGATCt       <  2:1535398/150‑1 (MQ=255)
                                                                                                                                          atACCAGCAGGTAGGCCGCGATTTTCATGTGTTTCTGAATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACCGGTTTTACTTGCTATGCCGCACCGGATGTCACGCTGGAAGATTATCt       <  2:1535405/149‑1 (MQ=255)
                                                                                                                                             ccAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACCGGTTTTACTTGCTATGCCGCACCGGATGTCACGCTGGAAGATGATCTTaa    >  1:878612/1‑150 (MQ=255)
                                                                                                                                             ccAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACCGGTTTTACTTGCTATGCCGCACCGGATGTCACGCTGGAAGATGATCTTaa    <  1:591857/150‑1 (MQ=255)
                                                                                                                                               agcagGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACCGGTTTTACTTGCTATGCCGCACCGGATGTCACGCTGGAAGATGATCTTAAGc  >  1:1418443/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
GTGTATTTGACGAAGGGATGATGGCGTGAAGATTTATCTGGTCGGTGGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACCGGTTTTACTTGCTATGCCGCACCGGATGTCACGCTGGAAGATGATCTTAAGC  >  NC_000913/3201866‑3202158

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