Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,031,651 G→A 8.6% L223F (CTC→TTC)  dcm ← DNA cytosine methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,031,6510GA 8.6% 297.3 / 13.0 93L223F (CTC→TTC) dcmDNA cytosine methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (44/41);  new base A (4/4);  total (48/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAATGTTTGCCGTCGATGATTTTCGGATCGTCTGGCCCATTATCTTCTGCATCAGCCACGTCATAGCCCAGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTCGATACGCCAGCCAGCGAAAATGGCTGACAAGGGAAACC  >  NC_000913/2031505‑2031798
                                                                                                                                                  |                                                                                                                                                   
aaaaTGTTTGCCGTCGATGATTTTCGGAACGTCTGGCCCATTATCTTCTGCATCAGCCACGTCATAGCCCAGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGcac                                                                                                                                                  >  1:1453141/1‑150 (MQ=255)
 aaaTGTTTGCCGTCGATGATTTTCGGATCGTCTGGCCCATTATCTTCTGCATCAGCCACGTCATAGCCCAGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGcaca                                                                                                                                                 >  2:659223/1‑150 (MQ=255)
       ttGCCGTCGATGATTTTCGGATCGTCTGGCCCATTATCTTCTGCATCAGCCACGTCATAGCCCAGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATc                                                                                                                                           <  2:578643/150‑1 (MQ=255)
          ccGTCGATGATTTTCGGATCGTCTGGCCCATTATCTTCTGCATCAGCCACGTCATAGCCCAGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGcc                                                                                                                                        <  2:1155785/150‑1 (MQ=255)
             tCGATGATTTTCGGATCGTCTGGCCCATTATCTTCTGCATCAGCCACGTCATAGCCCAGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGga                                                                                                                                     <  1:112711/150‑1 (MQ=255)
                       tCGGATCGTCTGGCCCATTATCTTCTGCATCAGCCACGTCATAGCCCAGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGt                                                                                                                           <  1:1493547/150‑1 (MQ=255)
                        cGGATCGGCTGGCCCATTATCTTCTGCATCAGCCACGTCATAGCCCAGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTc                                                                                                                          <  1:446098/150‑1 (MQ=255)
                         ggATCGTCTGGCCCATTATCTTCTGCATCAGCCACGTCATAGCCCAGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCg                                                                                                                         >  2:427632/1‑150 (MQ=255)
                           aTCGTCTGGCCCATTATCTTCTGCATCAGCCACGTCATAGCCCAGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGAt                                                                                                                       >  2:407245/1‑150 (MQ=255)
                           aTCGTCTGGCCCATTATCTTCTGCATCAGCCACGTCATAGCCCAGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGAt                                                                                                                       >  2:1181560/1‑150 (MQ=255)
                               tCTGGCCCATTATCTTCTGCATCAGCCACGTCATAGCCCAGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATg                                                                                                                   <  2:363313/150‑1 (MQ=255)
                               tCTGGCCCATTATCTTCTGCATCAGCCACGTCATAGCCCAGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATg                                                                                                                   >  1:319439/1‑150 (MQ=255)
                               tCTGGCCCATTATCTTCTGCATCAGCCACGTCAGAGCCCAGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATg                                                                                                                   <  2:351157/150‑1 (MQ=255)
                                 tGGCCCATTATCTTCTGCATCAGCCACGTCATAGCCCAGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCg                                                                                                                 <  1:63186/150‑1 (MQ=255)
                                 tGGCCCATTATCTTCTGCATCAGCCACGTCATAGCCCAGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCg                                                                                                                 <  2:132612/150‑1 (MQ=255)
                                 tGGCCCATTATCTTCTGCATCAGCCACGTCATAGCCCAGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCg                                                                                                                 <  2:725930/150‑1 (MQ=255)
                                  ggCCCATTATCTTCTGCATCAGCCACGTCATAGCCCAGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGt                                                                                                                <  2:754211/150‑1 (MQ=255)
                                  ggCCCATTATCTTCTGCATCAGCCACGTCATAGCCCAGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGt                                                                                                                <  2:570217/150‑1 (MQ=255)
                                    cccATTATCTTCTGCATCAGCCACGTCATAGCCCAGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTAc                                                                                                              >  1:1245894/1‑150 (MQ=255)
                                    cccATTATCTTCTGCATCAGCCACGTCATAGCCCAGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTAc                                                                                                              >  1:430801/1‑150 (MQ=255)
                                       aTTATCTTCTGCATCAGCCACGTCATAGCCCAGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACcac                                                                                                           <  1:1623116/150‑1 (MQ=255)
                                        ttATCTTCTGCATCAGCCACGTCATAGCCCAGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACcaca                                                                                                          >  2:410601/1‑150 (MQ=255)
                                        ttATCTTCTGCATCAGCCACGTCATAGCCCAGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACcaca                                                                                                          >  1:1120147/1‑150 (MQ=255)
                                        ttATCTTCTGCATCAGCCACGTCATAGCCCAGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACcaca                                                                                                          >  1:742301/1‑150 (MQ=255)
                                        ttATCTTCTGCATCAGCCACGTCATAGCCCAGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACcaca                                                                                                          >  1:420197/1‑150 (MQ=255)
                                             ttctGCATCAGCCACGTCATAGCCCAGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCaaa                                                                                                     >  1:528823/1‑150 (MQ=255)
                                             ttctGCATCAGCCACGTCATAGCCCAGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAACACAAACATCGCCGGACGCCGCGCGTCGATAATGCGTACCACATCaaa                                                                                                     >  1:1544589/1‑150 (MQ=255)
                                             ttctGCATCAGCCACGTCATAGCCCAGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAACACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCaaa                                                                                                     >  1:1544632/1‑150 (MQ=255)
                                             ttctGCATCAGCCACGTCATAGCCCAGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAACACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCaaa                                                                                                     >  1:751455/1‑150 (MQ=255)
                                             ttctGCATCAGCCACGTCATAGCCCAGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAACACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCaaa                                                                                                     >  1:751229/1‑150 (MQ=255)
                                              tctGCATCAGCCACGTCATAGCCCAGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCaaaa                                                                                                    <  1:833189/150‑1 (MQ=255)
                                                tGCATCAGCCACGTCATAGCCCAGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCaaaaaa                                                                                                  <  1:1279371/150‑1 (MQ=255)
                                                  cATCAGCCACGTCATAGCCCAGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACa                                                                                                >  1:1149619/1‑150 (MQ=255)
                                                    tCAGCCACGTCATAGCCCAGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGc                                                                                              >  2:149694/1‑150 (MQ=255)
                                                    tCAGCCACGTCATAGCCCAGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGc                                                                                              >  2:1562811/1‑150 (MQ=255)
                                                    tCAGCCACGTCATAGCCCAGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGc                                                                                              >  2:850020/1‑150 (MQ=255)
                                                    tCAGCCACGTCATAGCCCAGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGc                                                                                              >  2:849847/1‑150 (MQ=255)
                                                    tCAGCCACGTCATAGCCCAGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACGTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGc                                                                                              >  2:149692/1‑150 (MQ=255)
                                                             tCATAGCCCAGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTgg                                                                                     >  1:726002/1‑150 (MQ=255)
                                                             tCATAGCCCAGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTgg                                                                                     >  1:725992/1‑150 (MQ=255)
                                                             tCATAGCCCAGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTgg                                                                                     >  1:1423012/1‑150 (MQ=255)
                                                                   cccAGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCg                                                                               >  1:898349/1‑150 (MQ=255)
                                                                      aGTTCGTCCAGCTTCTGCATGATGATGCGGTACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGTCGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATC‑AAAAACAGCGTGCCCTGTGTATCGCAgt                                                                           <  2:977477/150‑2 (MQ=255)
                                                                      aGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAg                                                                            <  1:427632/150‑1 (MQ=255)
                                                                      aGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAg                                                                            <  1:976906/150‑1 (MQ=255)
                                                                       gTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAgg                                                                           <  2:446161/150‑1 (MQ=255)
                                                                        ttCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTTAAGTTTTCGAGCACAAACATCGCCGGTCGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGc                                                                          >  1:464852/1‑150 (MQ=255)
                                                                        ttCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGc                                                                          >  1:946336/1‑150 (MQ=255)
                                                                            tCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCaaaa                                                                      <  1:1672685/150‑1 (MQ=255)
                                                                              cAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAAcc                                                                    <  1:1531708/150‑1 (MQ=255)
                                                                                gCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCCAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGt                                                                  >  2:593898/1‑150 (MQ=255)
                                                                                gCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGt                                                                  <  1:55884/150‑1 (MQ=255)
                                                                                gCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGt                                                                  <  2:260774/150‑1 (MQ=255)
                                                                                   tCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTgcg                                                               <  1:1011072/150‑1 (MQ=255)
                                                                                   tCTGCATGATGATGCGGAACGATTTACCCTGGTCGTGACTATACAGGTTTTTGACGATTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTgcg                                                               <  2:1701352/150‑1 (MQ=255)
                                                                                                gCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCcgagcga                                                  <  2:447145/150‑1 (MQ=255)
                                                                                                gCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCcgagcga                                                  <  2:254875/150‑1 (MQ=255)
                                                                                                          ttACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTc                                        <  2:420197/150‑1 (MQ=255)
                                                                                                          ttACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTc                                        <  2:742301/150‑1 (MQ=255)
                                                                                                          ttACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTc                                        <  2:1120147/150‑1 (MQ=255)
                                                                                                          ttACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTc                                        <  1:410601/150‑1 (MQ=255)
                                                                                                          ttACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTc                                        <  2:1381900/150‑1 (MQ=255)
                                                                                                             cccTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAACACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAACAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTCGAt                                     <  2:1544589/150‑1 (MQ=255)
                                                                                                             cccTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAACACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTCGAt                                     <  2:1544632/150‑1 (MQ=255)
                                                                                                             cccTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAACACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTCGAt                                     <  2:751229/150‑1 (MQ=255)
                                                                                                             cccTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAACACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTCGAt                                     <  2:751455/150‑1 (MQ=255)
                                                                                                              ccTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTCGATa                                    <  1:731580/150‑1 (MQ=255)
                                                                                                              ccTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTCGATa                                    >  2:309026/1‑150 (MQ=255)
                                                                                                                     gTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTCGATACgccagc                             >  2:1260017/1‑150 (MQ=255)
                                                                                                                      tGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTCGATACgccagcc                            >  2:342692/1‑150 (MQ=255)
                                                                                                                      tGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTCGATACgccagcc                            >  1:1217382/1‑150 (MQ=255)
                                                                                                                      tGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTCGATACgccagcc                            >  2:1208861/1‑150 (MQ=255)
                                                                                                                          ttttCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTCGATACGCCAGCCAGCg                        >  2:38872/1‑150 (MQ=255)
                                                                                                                          ttttCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTCGATACGCCAGCCAGCg                        >  2:39067/1‑150 (MQ=255)
                                                                                                                           tttCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTCGATACGCCAGCCAGCGa                       >  2:1551137/1‑150 (MQ=255)
                                                                                                                              cAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTCGATACGCCAGCCAGCGaaaa                    <  1:83968/150‑1 (MQ=255)
                                                                                                                              cAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTCGATACGCCAGCCAGCGaaaa                    <  1:755613/150‑1 (MQ=255)
                                                                                                                              cAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTCGATACGCCAGCCAGCGaaaa                    <  1:755605/150‑1 (MQ=255)
                                                                                                                               aGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTCGATACGCCAGCCAGCGAAAAt                   >  2:983199/1‑150 (MQ=255)
                                                                                                                               aGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTCGATACGCCAGCCAGCGAAAAt                   >  2:107175/1‑150 (MQ=255)
                                                                                                                               aGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTCGATACGCCAGCCAGCGAAAAt                   >  2:203852/1‑150 (MQ=255)
                                                                                                                               aGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTCGATACGCCAGCCAGCGAAAAt                   >  2:1658949/1‑150 (MQ=255)
                                                                                                                                 gTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTCGATACGCCAGCCAGCGAAAATgg                 <  2:474809/150‑1 (MQ=255)
                                                                                                                                 gTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTCGATACGCCAGCCAGCGAAAATgg                 >  1:277534/1‑150 (MQ=255)
                                                                                                                                  tttttGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTCGATACGCCAGCCAGCGAAAATGGc                >  1:636317/1‑150 (MQ=255)
                                                                                                                                  tttttGACGTTTTCGAGCACAAACATCACCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTCGATACGCCAGCCAGCGAAAATGGc                <  2:1435288/150‑1 (MQ=255)
                                                                                                                                        aCGTTTTCGCGCACAAACACCGCCGGACGACGCGCGTAGATAATGCGTACCACATCACAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTCGATACGCCAGCCAGCGAAAATGGCTGACaa          <  1:742596/150‑1 (MQ=255)
                                                                                                                                        aCGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTCGATACGCCAGCCAGCGAAAATGGCTGACaa          <  1:1598263/150‑1 (MQ=255)
                                                                                                                                          gTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTCGATACGCCAGCCAGCGAAAATGGCTGACAAgg        >  2:1376812/1‑150 (MQ=255)
                                                                                                                                          gTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTCGATACGCCAGCCAGCGAAAATGGCTGACAAgg        >  2:911788/1‑150 (MQ=255)
                                                                                                                                           ttttCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGCATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTCGATACGCCAGCCAGCGAAAATGGCTGACAAggg       <  1:267436/150‑1 (MQ=255)
                                                                                                                                                gAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTCGATACGCCAGCCAGCGAAAATGGCTGACAAGGGAAAcc  >  1:744152/1‑150 (MQ=255)
                                                                                                                                                gAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTCGATACGCCAGCCAGCGAAAATGGCTGACAAGGGAAAcc  >  1:1719275/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AAAATGTTTGCCGTCGATGATTTTCGGATCGTCTGGCCCATTATCTTCTGCATCAGCCACGTCATAGCCCAGTTCGTCCAGCGTCTGCATGATGATGCGGAACGTTTTACCCTGGTCGTGACTTTTCAGGTTTTTGACGTTTTCGAGCACAAACATCGCCGGACGACGCGCGTCGATAATGCGTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTCGATACGCCAGCCAGCGAAAATGGCTGACAAGGGAAACC  >  NC_000913/2031505‑2031798

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