Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,449,799 C→G 5.8% T517T (ACG→ACC tynA ← tyramine oxidase, copper‑requiring

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,449,7990CG 5.8% 365.0 / 12.8 121T517T (ACG→ACCtynAtyramine oxidase, copper‑requiring
Reads supporting (aligned to +/- strand):  ref base C (57/56);  new base G (3/4);  total (61/60)
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.48e-01

Only 100 of 123 total aligned reads displayed.

CTACATCCAGATCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGATGCTTA  >  NC_000913/1449651‑1449946
                                                                                                                                                    |                                                                                                                                                   
cTACATCCAGATCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCg                                                                                                                                                    >  1:854073/1‑150 (MQ=255)
cTACATCCAGATCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCg                                                                                                                                                    <  1:1618290/150‑1 (MQ=255)
cTACATCCAGATCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCg                                                                                                                                                    >  2:1295652/1‑150 (MQ=255)
 tACATCCAGATCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGt                                                                                                                                                   >  2:1474531/1‑150 (MQ=255)
  aCATCCAGATCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTa                                                                                                                                                  >  2:1714368/1‑150 (MQ=255)
    aTCCAGATCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACCCCTTTCACCGCTTTGATGCCCGTAGc                                                                                                                                                >  2:1804449/1‑150 (MQ=255)
    aTCCAGATCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGc                                                                                                                                                <  1:31248/150‑1 (MQ=255)
     tCCAGATCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGATTTAACACCTTTCACCGCTTCGATGCCCGTAGCa                                                                                                                                               >  1:1254197/1‑150 (MQ=255)
      ccAGATCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCAc                                                                                                                                              >  1:1092444/1‑150 (MQ=255)
        aGATCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCg                                                                                                                                            <  2:1119243/150‑1 (MQ=255)
           tagtggCGGAAATTATATATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGGCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGTCa                                                                                                                                         <  2:1250295/146‑1 (MQ=255)
           tCGTGGCTGATAGTATTAATTTGTTGGTGTGTAGTACAGAAGATATTGTGTTCGTTTAGAGTTCGGTAGCGCGTGTCATATTTCGCCGTCTCATCGTGCATGGTTTTCGCTGTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCa                                                                                                                                         <  1:349262/150‑1 (MQ=255)
           tCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCa                                                                                                                                         <  1:1398846/150‑1 (MQ=255)
             gAGGCTGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCGGTGCTACGGGCATc                                                                                                                                       <  2:1515795/150‑1 (MQ=255)
             gAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATc                                                                                                                                       <  1:1725879/150‑1 (MQ=255)
              aGGCGGAAATTATAAATCTGTTGGTGTGTAGTACCCACGATATTGTGATCGATAATCTTGACGTAGCTCGTGTCATCTTTCGCCGTCTCATAGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCg                                                                                                                                      <  2:314254/150‑1 (MQ=255)
                gCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGAt                                                                                                                                    >  2:258420/1‑150 (MQ=255)
                 cGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCATCGCTTCGATGCCCGTAGCACCGGCATCGATg                                                                                                                                   >  2:1381956/1‑150 (MQ=255)
                 cGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCGGTGCTACGGGCATCGaag                                                                                                                                   <  2:1560011/150‑3 (MQ=255)
                  ggAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGc                                                                                                                                  >  1:256472/1‑150 (MQ=255)
                  ggAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGAGCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCATCGATGCCGGTGCTACGGGCATCGAAGc                                                                                                                                  <  1:1313436/150‑1 (MQ=255)
                        tataAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCTGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAg                                                                                                                            >  2:1825915/1‑150 (MQ=255)
                        tataAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAg                                                                                                                            <  1:1194893/150‑1 (MQ=255)
                        tataAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTGTCGCTTTAACACCTTTCACCGATTCGATGCCCGTAGCACCGGCATCGATGCCAATAg                                                                                                                            >  1:347734/1‑150 (MQ=255)
                          taAATTTGGTGGTGTGTAGTACCCACGATATTGTGTTCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTGATCGTGCATGGTTTTCGGTTTAACACCTTGCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTg                                                                                                                          <  2:669405/150‑1 (MQ=255)
                           aaaTATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCCCCGCTTCGATGCCCGTCGCACCTGCGGCGATGCCAATAGGGc                                                                                                                         >  2:949835/1‑150 (MQ=255)
                           aaaTATGTTGGTGGGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGAACCGGCATCGATGCCAATAGTGc                                                                                                                         <  1:139673/150‑1 (MQ=255)
                            aaTATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTGTCGCTTTAACACCTTTCCCCGCTTCGATGCCAGTAGCACCGGCATCGATGCCAATAGTGcc                                                                                                                        >  1:340091/1‑150 (MQ=255)
                                     gtgtgtAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAAGACCTTTCACCGCTTCGTTGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATg                                                                                                               <  2:28696/150‑1 (MQ=255)
                                     gtgtgtAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATg                                                                                                               <  1:1900569/150‑1 (MQ=255)
                                     gtgtgtAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATg                                                                                                               <  1:419704/150‑1 (MQ=255)
                                     gggtgtAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCCTCTCAGCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATg                                                                                                               <  2:342960/148‑1 (MQ=255)
                                      tgtgtAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATgg                                                                                                              <  2:1903657/150‑1 (MQ=255)
                                      ggtgtAGTACCCACGATATTGTGATCGATAAGCGGGCCGTAGCGCGGGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATgg                                                                                                              <  2:1387843/149‑1 (MQ=255)
                                            gtaCCCACGATATTGTGATCGATAAGCGGGCCGTAGCGCGTGGCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATc                                                                                                        <  2:1245346/150‑1 (MQ=255)
                                             taCACACGATATTGTGATCGATAAGCTTGCCGAAGCGCGGGTAATCTATCGCCGTCTCATCGTGCATGGGTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGTCATCGATGCCAAAAGTGCCGTTTTCATGGAAGATcc                                                                                                       <  1:1573966/150‑1 (MQ=255)
                                              aCCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCa                                                                                                      <  1:1197072/150‑1 (MQ=255)
                                               cccACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAg                                                                                                     >  2:1558186/1‑150 (MQ=255)
                                                         tgtgATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGt                                                                                           <  2:256472/150‑1 (MQ=255)
                                                         tgtgATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGt                                                                                           >  2:191259/1‑150 (MQ=255)
                                                          gtgATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTa                                                                                          >  2:91258/1‑150 (MQ=255)
                                                             aTCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGACCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTc                                                                                       <  2:1771022/150‑1 (MQ=255)
                                                                  tAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGt                                                                                  <  1:1571553/150‑1 (MQ=255)
                                                                    aGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTa                                                                                >  1:289609/1‑150 (MQ=255)
                                                                    aGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGACGTTTTCATGGAAGATCAAGTCAAAAATGTAGTCATAGTTa                                                                                >  2:1225200/1‑150 (MQ=255)
                                                                     gCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGATGATCCAGTCAAAAATGTAGTCATAGTTAc                                                                               >  1:1019887/1‑150 (MQ=255)
                                                                       gTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTAccc                                                                             <  2:1092444/150‑1 (MQ=255)
                                                                         gCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCAc                                                                           <  1:1295652/150‑1 (MQ=255)
                                                                           cGTAGCGCGTGTCATCTTTAGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTg                                                                         <  1:1770284/150‑1 (MQ=255)
                                                                              aGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTAc                                                                      <  2:966186/150‑1 (MQ=255)
                                                                              aGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTAc                                                                      <  1:3482/150‑1 (MQ=255)
                                                                                cgcgTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTCCCCACAGTACTg                                                                    >  2:145545/1‑150 (MQ=255)
                                                                                cgcgTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTg                                                                    >  1:1968175/1‑150 (MQ=255)
                                                                                    tgtCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATcc                                                                <  1:599406/150‑1 (MQ=255)
                                                                                    tgtCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATcc                                                                >  1:887520/1‑150 (MQ=255)
                                                                                    tgtCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATcc                                                                <  2:241622/150‑1 (MQ=255)
                                                                                      tcttCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAg                                                              <  2:1993173/147‑1 (MQ=255)
                                                                                      tCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAg                                                              >  2:501440/1‑150 (MQ=255)
                                                                                       cATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAgc                                                             >  2:1830633/1‑150 (MQ=255)
                                                                                         tCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAgcgc                                                           <  1:1410471/150‑1 (MQ=255)
                                                                                         acTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGATTCGATGCCCGTAGCACCGGCATAGATGCCAATAGTGCCGTTTTCATGGAAGATACAGTCAAAAATGTAGTCATAGTTACCGACTGTACTGATCCAgcgc                                                           >  1:71158/2‑150 (MQ=255)
                                                                                            ttCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTTCTGATCCAGCGcacc                                                        >  1:1279116/1‑150 (MQ=255)
                                                                                            ttCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGcacc                                                        >  1:805162/1‑150 (MQ=255)
                                                                                            ttCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGcacc                                                        >  1:797947/1‑150 (MQ=255)
                                                                                            ttCGCCGTCTCATCGAGCATGGTATTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTATCACCGGGATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGcacc                                                        <  2:315632/150‑1 (MQ=255)
                                                                                                ccGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGGTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTa                                                    <  2:334089/150‑1 (MQ=255)
                                                                                                  gTCTCAGCGTGCATTGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCAACGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAAc                                                  <  1:1015768/150‑1 (MQ=255)
                                                                                                      cATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGCTGCCCGTAGCACCGGCATCGCTGCCAATAGTGCCGTTTTCATGGAAGCTCCAGGCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTccc                                              >  1:1228614/1‑150 (MQ=255)
                                                                                                       aTCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCg                                             <  1:1585168/150‑1 (MQ=255)
                                                                                                            gCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGtt                                        >  1:605358/1‑150 (MQ=255)
                                                                                                              aTGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCg                                      <  1:1027872/150‑1 (MQ=255)
                                                                                                               tGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCgg                                     <  1:1714368/150‑1 (MQ=255)
                                                                                                               tGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCgg                                     <  1:608697/150‑1 (MQ=255)
                                                                                                                  ttttCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTAc                                  <  1:91258/150‑1 (MQ=255)
                                                                                                                  ttttCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTAc                                  >  1:869846/1‑150 (MQ=255)
                                                                                                                     tCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGa                               <  1:1676704/150‑1 (MQ=255)
                                                                                                                      cGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTGCGGTACTGAc                              >  2:1981489/1‑150 (MQ=255)
                                                                                                                      cGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATACAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGTTACTGAc                              >  1:569428/1‑150 (MQ=255)
                                                                                                                       gCTTTAACACCTTTCACCGCTTCGATGCACGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCAAGCGCACCACTAACTCCCGGCGATCGGTACTGTCg                             >  1:105093/1‑150 (MQ=255)
                                                                                                                          ttAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTg                          >  1:1035144/1‑150 (MQ=255)
                                                                                                                          ttAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTg                          >  1:1911348/1‑150 (MQ=255)
                                                                                                                             acacCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGGTCGGTACTGACGTTGGGc                       >  1:874282/1‑150 (MQ=255)
                                                                                                                                ccTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTgg                    >  2:69508/1‑150 (MQ=255)
                                                                                                                                  tttCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGcc                  <  2:805162/150‑1 (MQ=255)
                                                                                                                                  tttCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGcc                  <  2:797947/150‑1 (MQ=255)
                                                                                                                                   ttCACCACTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTAGTGACGTTGGGCTGGccc                 >  2:1170259/1‑150 (MQ=255)
                                                                                                                                     cACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCAt               >  1:915745/1‑150 (MQ=255)
                                                                                                                                      aCCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCCCCCCTAACTCCCGGCGTGCGGTCCTGACGTTGGGCTGGCCCAtt              >  2:1704366/1‑150 (MQ=255)
                                                                                                                                         gCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTcc           <  1:1474531/150‑1 (MQ=255)
                                                                                                                                          catCGATGCCGGTGCTACGGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCt          >  2:1624597/3‑150 (MQ=255)
                                                                                                                                          cTTCGATGCCCGTAGCACCGTCATCGATGCCAATAGTGCCGTTTTAATGGAAGATCCAGTCAAATATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCAGCACTAACTCCCGGCGTTCGGTACTGACGTTTGGCTGGCCCATTTTCt          <  1:301297/150‑1 (MQ=255)
                                                                                                                                          cTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCt          <  2:964447/150‑1 (MQ=255)
                                                                                                                                          cTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCt          <  2:1742364/150‑1 (MQ=255)
                                                                                                                                          cTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCt          <  2:1943048/150‑1 (MQ=255)
                                                                                                                                           atCGATGCCGGTGCTACGGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTGCTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTTCTGACGTTGGGCTGGCCCATTTCCTg         >  1:1300324/2‑150 (MQ=255)
                                                                                                                                             cGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTATTCATGGAAGATCCAGTCAAAAATGTACTCATAGTTAACCACAGTACTGATCAAGCGCACCACTAACTCCCTGCGTTCGGAACTAAAGATGGACTGTACCATTTACTGAt       >  2:1593830/1‑150 (MQ=255)
                                                                                                                                             cGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGGTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGCCGTTGGGCTGGCCCATGTCCTGAt       >  1:598983/1‑150 (MQ=255)
                                                                                                                                              gATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGATg      >  2:1727932/1‑150 (MQ=255)
                                                                                                                                                tGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGATGCt    <  2:1039852/150‑1 (MQ=255)
                                                                                                                                                  cccGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGATGCTta  >  1:1653636/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CTACATCCAGATCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGATGCTTA  >  NC_000913/1449651‑1449946

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: