Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,518,082 C→G 6.4% noncoding (35/76 nt) alaX ← tRNA‑Ala

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,518,0820CG 6.4% 221.2 / 2.5 69noncoding (35/76 nt)alaXtRNA‑Ala
Reads supporting (aligned to +/- strand):  ref base C (14/50);  new base G (2/3);  total (16/53)
Fisher's exact test for biased strand distribution p-value = 3.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAG  >  NC_000913/2517933‑2518221
                                                                                                                                                     |                                                                                                                                           
tCCAATGATGCTTTCAGCAGGCGGCACTAACCCGGACGTTCACTGTATTACAAAAGCGGCAAAAAGCAGACACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCCAACCGCTGACCTCTTGcatgc                                                                                                                                             <  2:1225666/150‑1 (MQ=38)
  cAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcatgcca                                                                                                                                           <  2:984161/150‑1 (MQ=255)
    atgatgCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGTTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcatgccatg                                                                                                                                         <  2:1368969/150‑1 (MQ=255)
    atgatgCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcatgccatg                                                                                                                                         <  2:1114399/150‑1 (MQ=255)
        tgCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAg                                                                                                                                     <  1:604635/150‑1 (MQ=255)
          cTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAgcg                                                                                                                                   >  1:1232528/1‑150 (MQ=255)
           tttCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGGCATGCAAgagg                                                                                                                                  <  1:382330/150‑4 (MQ=255)
           tttCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAgcgc                                                                                                                                  <  1:1101363/150‑1 (MQ=255)
              cACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCAGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCAATGCAAGCGctct                                                                                                                               >  2:865073/1‑150 (MQ=255)
               aCCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCtctc                                                                                                                              >  1:1510210/1‑150 (MQ=255)
               aCCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAAACCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGAAAGCGAtctc                                                                                                                              >  2:752430/1‑150 (MQ=255)
                 cAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTccc                                                                                                                            <  1:76872/150‑1 (MQ=255)
                   ggcggcACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGTAAATTTGTTGGAGCTAAGCGGGATCTAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAg                                                                                                                          <  2:648174/150‑1 (MQ=255)
                    gcggcACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGc                                                                                                                         <  1:1067529/150‑1 (MQ=255)
                     cggcACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCt                                                                                                                        >  1:830908/1‑150 (MQ=255)
                        cACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAg                                                                                                                     <  2:8914/150‑1 (MQ=255)
                        cACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAg                                                                                                                     <  1:1342724/150‑1 (MQ=255)
                        cACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAg                                                                                                                     >  1:764235/1‑150 (MQ=255)
                          cTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCt                                                                                                                   <  2:1476878/150‑1 (MQ=255)
                            aaCCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAACAACCCCCGCTTTGCAGCGAGGGTTGGAAATTGGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCtat                                                                                                                 <  2:693160/150‑1 (MQ=255)
                             aCCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCtata                                                                                                                <  1:1526074/150‑1 (MQ=255)
                               ccGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGc                                                                                                              <  1:1477468/150‑1 (MQ=255)
                                cGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGcc                                                                                                             <  1:840129/150‑1 (MQ=255)
                                cGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGcc                                                                                                             <  1:321593/150‑1 (MQ=255)
                                  gACGTAGACTGTATTACAAAAGCGGCAAAACGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGcccc                                                                                                           <  2:330990/150‑1 (MQ=255)
                                  gACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGcccc                                                                                                           >  1:757921/1‑150 (MQ=255)
                                    cGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCAc                                                                                                         <  2:357826/150‑1 (MQ=255)
                                     gTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACg                                                                                                        >  2:1273205/1‑150 (MQ=255)
                                          cTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCg                                                                                                   <  1:1273205/150‑1 (MQ=255)
                                            gTATTACAAAAGCGGCAAAACGCAGAGACAAAAACCCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGtt                                                                                                 <  2:968907/150‑1 (MQ=255)
                                            gTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGtt                                                                                                 <  1:1301676/150‑1 (MQ=255)
                                               ttACAAAAGCGGCAAAAAGCAGAGACAAAAAAGCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTAc                                                                                              <  1:706469/150‑1 (MQ=255)
                                                tACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACg                                                                                             <  1:1484455/150‑1 (MQ=255)
                                                tACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACg                                                                                             >  1:1078449/1‑150 (MQ=255)
                                                tACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACg                                                                                             >  2:1074877/1‑150 (MQ=255)
                                                  cAAAAGCCGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTa                                                                                           <  2:1466575/150‑1 (MQ=255)
                                                        cGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGt                                                                                     <  2:296424/150‑1 (MQ=255)
                                                            aaaaaGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGc                                                                                 <  2:830908/150‑1 (MQ=255)
                                                                aGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCGACGATGCGGTGACGTACCAAGTGTGCTggg                                                                             >  1:536553/1‑150 (MQ=255)
                                                                aGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTggg                                                                             <  1:852847/150‑1 (MQ=255)
                                                                aGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTggg                                                                             <  1:362120/150‑1 (MQ=255)
                                                                aGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGAAGAGTTTACGTACCAAGTTTGGTgcg                                                                             >  1:1221349/1‑148 (MQ=255)
                                                                   agagACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGc                                                                          <  2:764235/150‑1 (MQ=255)
                                                                   agagACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGc                                                                          <  1:834172/150‑1 (MQ=255)
                                                                      gaCAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCaaa                                                                       >  1:1344043/1‑150 (MQ=255)
                                                                            aaaCCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTtgg                                                                 <  1:865073/150‑1 (MQ=255)
                                                                             aaCCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTtggt                                                                <  2:1078449/150‑1 (MQ=255)
                                                                              aCCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTtggtg                                                               <  1:1074877/150‑1 (MQ=255)
                                                                              aCCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTtggtg                                                               <  2:1344043/150‑1 (MQ=255)
                                                                              aCCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTtggtg                                                               <  1:862813/150‑1 (MQ=255)
                                                                               cccccGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTtggtgg                                                              <  2:1136352/150‑1 (MQ=255)
                                                                                ccccGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGa                                                             <  2:443658/150‑1 (MQ=255)
                                                                                ccccGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTAAGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGa                                                             >  1:1095784/1‑150 (MQ=255)
                                                                                ccccGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCT‑CCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAg                                                            <  2:67099/150‑1 (MQ=255)
                                                                                  ccGCTTTCCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGc                                                           <  2:1370452/150‑1 (MQ=255)
                                                                                   cGCTTTGCAGCGAGGGTTTGAAATTTGGTGGAGCTAAGCGGGATCGAACCTCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCTAAATTTGGTGGAGCt                                                          <  2:287983/150‑1 (MQ=255)
                                                                                   cGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCTAAATTTGGTGGAGCt                                                          <  2:1457719/150‑1 (MQ=255)
                                                                                   cGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCt                                                          <  2:342598/150‑1 (MQ=255)
                                                                                    gCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTa                                                         <  1:1248669/150‑1 (MQ=255)
                                                                                    gCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTa                                                         <  1:1254071/150‑1 (MQ=255)
                                                                                             cGAGGGTTTGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCGGAGCTATAGCCCCACGATGCGTTGACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGagc                                                <  1:328780/150‑3 (MQ=255)
                                                                                             cGAGGGTTGGAAATTTGGTTGATCTAATCTGGATCGAACCGCTGACCTCTTCCATGCCATGATAGCGCTCTCCCATCTGAGCTATATCCACACCATGCTTTTACGTACCATGTTTGCTGGGTTCTAAATTTGGTGGAGCTAAGCGGGATc                                                <  2:1270743/150‑1 (MQ=21)
                                                                                              gAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCg                                               <  1:1026653/150‑1 (MQ=255)
                                                                                                 ggTTGGAAATTTGGTGGATCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATCCGTTTACGTACCAAGCTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAAc                                            <  1:1001678/150‑1 (MQ=255)
                                                                                                        aaTTTTGTGGATCTAAGCGGGATCGAACCTCTGAGCTCTTTCATGCAATGCAAGCGCTCTCCAAGCTGAGCTATAGCCCCACTATTCGTTTACGTACCAAGTTTGGTGGGTGGTAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGAc                                     <  2:525272/150‑1 (MQ=37)
                                                                                                                    cTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcatgcc                         <  1:9810/150‑1 (MQ=255)
                                                                                                                           ggATCGAACCTCGGAGCTCTTGCATGGAATGAAAGCGCTCTCCAAGCTGAGCTATAGCCCCACGATGCGTTTACGCACCAAGCTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAg                  <  1:1140631/150‑1 (MQ=25)
                                                                                                                              tCGAACCGCTGTCCTCTTTCATGGCATTCAAGCGCTCTACCAGCTGAGCTATAGCCCCACTATGCGTTTACGTACGAATTTTGCTGGTTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGGCATGCAAgagg               <  2:528670/150‑4 (MQ=25)
                                                                                                                              tCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCATCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGGCATGCAAgagg               <  1:752430/150‑4 (MQ=255)
                                                                                                                               cGAACCGCTGACCTCTTGCATGCCATGAAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGCTGCGTTTACGTACCAAGTTGGCTGGGTGCAAAATTTGCTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGGCATGCAAgaggt              <  1:446971/150‑5 (MQ=37)
                                                                                                                                   ccccgcaCCTCTAGCACGCAATGAAAGCGCTCTCCAAGCGGAGCTATATCCACACTATTCGTTCACGTACCAAGTATGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTcc          <  1:969792/144‑1 (MQ=18)
                                                                                                                                           cgctTGCATGGCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTTTTGCATGCCATGCAAGCGCTCTCCCAGCTGAg  >  2:1209733/3‑150 (MQ=255)
                                                                                                                                           cgctTGCATGGCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAg  >  1:1374702/3‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                           
TCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAG  >  NC_000913/2517933‑2518221

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