Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 357,156 T→A 6.6% L329Q (CTG→CAG)  codA → cytosine/isoguanine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,1560TA 6.6% 470.9 / 16.7 151L329Q (CTG→CAG) codAcytosine/isoguanine deaminase
Reads supporting (aligned to +/- strand):  ref base T (69/72);  new base A (6/4);  total (75/76)
Fisher's exact test for biased strand distribution p-value = 5.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.38e-01

Only 100 of 164 total aligned reads displayed.

CAATATTCATCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAGCGCCAACCTGATTA  >  NC_000913/357007‑357305
                                                                                                                                                     |                                                                                                                                                     
cAATATTCATCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATAtgct                                                                                                                                                       >  1:1030898/1‑150 (MQ=255)
  atatTCATCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATAtgctgc                                                                                                                                                     <  1:1853932/150‑1 (MQ=255)
     ttCATCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAg                                                                                                                                                  >  1:220476/1‑150 (MQ=255)
     ttCATCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAg                                                                                                                                                  <  1:883996/150‑1 (MQ=255)
      tCATCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTACAAGt                                                                                                                                                 <  2:812524/150‑1 (MQ=255)
       cATCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGtg                                                                                                                                                >  1:2472852/1‑150 (MQ=255)
        aTCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATgcagcacttgc                                                                                                                                               <  1:195285/150‑6 (MQ=255)
        aTCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGtgc                                                                                                                                               <  1:2446881/150‑1 (MQ=255)
          cTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGtgctg                                                                                                                                             <  1:1481159/150‑1 (MQ=255)
             cAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCat                                                                                                                                          >  1:1785892/1‑150 (MQ=255)
              aaGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGGCCGGCATTAACGTCTGCTTTGGTCAGGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCata                                                                                                                                         <  1:1673770/150‑1 (MQ=255)
                ggACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATg                                                                                                                                       <  2:1271478/150‑1 (MQ=255)
                  aCGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCAGCACTTGCAGCatattcg                                                                                                                                     <  2:1084909/150‑4 (MQ=255)
                   cGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATgggg                                                                                                                                    >  2:1662915/1‑150 (MQ=255)
                   cGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATgggg                                                                                                                                    >  2:1546218/1‑150 (MQ=255)
                   cGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCAGCACTTGCAGCatattcgc                                                                                                                                    <  2:1044123/150‑5 (MQ=255)
                    gTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATGCGCTGGGAAGGGCGAATATGCTGCAAGTGCTGCATATGGGGc                                                                                                                                   <  2:2027142/150‑1 (MQ=255)
                      ttCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTg                                                                                                                                 >  1:803961/1‑150 (MQ=255)
                        cGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCa                                                                                                                               >  2:1382691/1‑150 (MQ=255)
                        cGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCa                                                                                                                               <  1:643551/150‑1 (MQ=255)
                         gATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCAt                                                                                                                              >  1:141814/1‑150 (MQ=255)
                                aTCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGGTGCATGTTTGcc                                                                                                                       <  1:548716/150‑1 (MQ=255)
                                 tCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCa                                                                                                                      >  1:1951991/1‑150 (MQ=255)
                                   cAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGt                                                                                                                    >  1:1745140/1‑150 (MQ=255)
                                   cAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGt                                                                                                                    >  1:1727318/1‑150 (MQ=255)
                                   cAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGt                                                                                                                    >  2:208070/1‑150 (MQ=255)
                                   cAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGt                                                                                                                    <  1:855590/150‑1 (MQ=255)
                                   cAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGt                                                                                                                    >  1:1097751/1‑150 (MQ=255)
                                    aaaaCGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGtt                                                                                                                   >  1:1170128/1‑150 (MQ=255)
                                      aaCGTCGCGGCATCACGTGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGa                                                                                                                 <  1:1413986/150‑1 (MQ=255)
                                            gcgGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCt                                                                                                           <  2:386135/150‑1 (MQ=255)
                                             cgGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTa                                                                                                          >  2:2312343/1‑150 (MQ=255)
                                              ggCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTAc                                                                                                         <  2:105387/150‑1 (MQ=255)
                                              ggCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTAc                                                                                                         <  1:980693/150‑1 (MQ=255)
                                              ggCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCAGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTAc                                                                                                         <  1:1119970/150‑1 (MQ=255)
                                              ggCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGAAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTAc                                                                                                         <  2:2153737/150‑1 (MQ=255)
                                                 aTCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACggg                                                                                                      >  2:1951551/1‑150 (MQ=255)
                                                    ccgcgcgTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAg                                                                                                   <  2:976030/149‑1 (MQ=255)
                                                         cgTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTaa                                                                                              >  1:2209304/1‑150 (MQ=255)
                                                         cgTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTaa                                                                                              >  1:1207573/1‑150 (MQ=255)
                                                               agagATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATgg                                                                                        <  2:803961/150‑1 (MQ=255)
                                                               agagATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATgg                                                                                        <  2:1187905/150‑1 (MQ=255)
                                                               agagATGCTGGAGTCCGGCATTAACGTCTGCTTTGGGCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATgg                                                                                        <  2:2218224/150‑1 (MQ=255)
                                                                gagaTGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGc                                                                                       >  1:113744/1‑150 (MQ=255)
                                                                 agaTGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGcc                                                                                      <  2:917935/150‑1 (MQ=255)
                                                                 agaTGCTGGAGTCCGGCATTAACGGCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACTATGGcc                                                                                      <  2:306137/150‑1 (MQ=255)
                                                                   aTGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTTCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTg                                                                                    <  2:1174208/150‑1 (MQ=255)
                                                                    tGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGa                                                                                   >  1:1082958/1‑150 (MQ=255)
                                                                      cTGGAGTCCGTCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTAGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAAt                                                                                 <  2:2005787/150‑1 (MQ=255)
                                                                        ggAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAAtt                                                                                >  1:1847163/1‑149 (MQ=255)
                                                                        ggAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAAtt                                                                                <  2:1847163/149‑1 (MQ=255)
                                                                           gTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAAt                                                                            <  1:2312343/150‑1 (MQ=255)
                                                                            tCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATc                                                                           <  2:1707285/150‑1 (MQ=255)
                                                                             ccGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATca                                                                          >  1:903484/1‑150 (MQ=255)
                                                                             ccGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATca                                                                          >  1:690405/1‑150 (MQ=255)
                                                                               ggCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATcacc                                                                        >  1:954088/1‑150 (MQ=255)
                                                                               ggCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATcacc                                                                        >  1:1653223/1‑150 (MQ=255)
                                                                                gCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATcaccc                                                                       >  1:1061332/1‑150 (MQ=255)
                                                                                gCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATcaccc                                                                       <  1:498663/150‑1 (MQ=255)
                                                                                 cATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATcaccca                                                                      >  1:1398413/1‑150 (MQ=255)
                                                                                     aaCGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTCATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCcaccac                                                                  >  1:407569/1‑150 (MQ=255)
                                                                                       cGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAg                                                                >  1:1355619/1‑150 (MQ=255)
                                                                                           tGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCa                                                            <  2:1622051/150‑1 (MQ=255)
                                                                                            gCTTTGTTCACGATGATGTGTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGTATGTTTGCCAGTTGATGGGCTACTGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCaa                                                           <  2:2339518/150‑1 (MQ=255)
                                                                                            gCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCaa                                                           >  2:1588241/1‑150 (MQ=255)
                                                                                            gCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCaa                                                           >  1:2319540/1‑150 (MQ=255)
                                                                                              tttGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAgg                                                         <  2:1412060/150‑1 (MQ=255)
                                                                                                  gTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGt                                                     >  2:2186210/1‑150 (MQ=255)
                                                                                                    cACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCTAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTg                                                   <  2:1398413/150‑1 (MQ=255)
                                                                                                    cACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTg                                                   <  1:2427789/150‑1 (MQ=255)
                                                                                                       gatgatGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAAt                                                <  1:740697/150‑1 (MQ=255)
                                                                                                         tgatgTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAAttt                                              <  2:1779795/150‑1 (MQ=255)
                                                                                                            tgTCTTCGATCCGTGGTATCCGCTGGGAACTGGGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCa                                           <  1:224590/150‑1 (MQ=255)
                                                                                                            tgTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCa                                           >  2:1089454/1‑150 (MQ=255)
                                                                                                            tgTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCa                                           <  2:1785892/150‑1 (MQ=255)
                                                                                                             gTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAg                                          >  1:134224/1‑150 (MQ=255)
                                                                                                              gcTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAgg                                         <  2:1302169/149‑1 (MQ=255)
                                                                                                               cTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGa                                        >  1:193296/1‑150 (MQ=255)
                                                                                                               cTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGa                                        >  2:1287722/1‑150 (MQ=255)
                                                                                                               cTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGa                                        >  2:657359/1‑150 (MQ=255)
                                                                                                                 tCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGAtt                                      >  2:1104501/1‑150 (MQ=255)
                                                                                                                    aTCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATTGGCTACGGGAAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACg                                   <  2:2472852/150‑1 (MQ=255)
                                                                                                                     tCCGTGGTATCCGCTGGGAACGGCTAATATTATTCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGTCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACgg                                  <  2:871879/150‑1 (MQ=255)
                                                                                                                       cGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCa                                <  2:1406158/150‑1 (MQ=255)
                                                                                                                       cGTGGTAACCGCTCGGAACGGCGAATATGCTGCAAGAGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCAAAGCGCAAGCACGTTGAATTTGCAGGATTACggaa                                >  2:1907311/1‑148 (MQ=255)
                                                                                                                        gTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCAt                               <  1:1546218/150‑1 (MQ=255)
                                                                                                                             aTCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTgccg                          <  2:1750334/150‑1 (MQ=255)
                                                                                                                                cGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGCGGCTGCATGTTTCCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGccgccg                       >  2:904128/1‑150 (MQ=255)
                                                                                                                                 gCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCgg                      >  2:2095544/1‑150 (MQ=255)
                                                                                                                                   tGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGaa                    <  2:1170128/150‑1 (MQ=255)
                                                                                                                                          ggCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAgcgc             <  2:2209304/150‑1 (MQ=255)
                                                                                                                                           gCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAGCGcc            >  2:2412213/1‑150 (MQ=255)
                                                                                                                                           ccccatatGCAGCACTTGCAGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAGCGcc            >  1:678338/5‑150 (MQ=255)
                                                                                                                                            cccatatGCAGCACTTGCAGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAGCGCCa           >  2:2386313/4‑150 (MQ=255)
                                                                                                                                              catatGCAGCACTTGCAGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAGCGCCAAc         >  1:1139480/2‑150 (MQ=255)
                                                                                                                                              aaTATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAGCGCCAAc         <  2:732696/150‑1 (MQ=255)
                                                                                                                                               atatGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAGCGCCAAcc        <  2:1712534/150‑1 (MQ=255)
                                                                                                                                                   gctgcAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACGACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAGCGCCAACCTGAt    >  1:2445094/1‑150 (MQ=255)
                                                                                                                                                    ctgcAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAGCGCCAACCTGAtt   >  1:1191250/1‑150 (MQ=255)
                                                                                                                                                    ctgcAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAGCGCCAACCTGAtt   >  1:26061/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAATATTCATCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAGCGCCAACCTGATTA  >  NC_000913/357007‑357305

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