Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 357,160 A→C 5.7% Q330H (CAA→CAC codA → cytosine/isoguanine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,1600AC 5.7% 483.4 / 17.8 159Q330H (CAA→CACcodAcytosine/isoguanine deaminase
Reads supporting (aligned to +/- strand):  ref base A (75/75);  new base C (6/3);  total (81/78)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.62e-01

Only 100 of 161 total aligned reads displayed.

TTCATCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAGCGCCAACCTGATTATCC  >  NC_000913/357012‑357308
                                                                                                                                                    |                                                                                                                                                    
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)
  cATCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGtg                                                                                                                                                   >  1:1412060/1‑150 (MQ=255)
   aTCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTGTCCGCTGGGAACGGCGAATATGCTGCAAGtgc                                                                                                                                                  >  1:1528115/1‑150 (MQ=255)
      tGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGtgctgc                                                                                                                                               <  2:1146679/150‑1 (MQ=255)
        cAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCat                                                                                                                                             >  1:1785892/1‑150 (MQ=255)
         aaGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGGCCGGCATTAACGTCTGCTTTGGTCAGGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCata                                                                                                                                            <  1:1673770/150‑1 (MQ=255)
             aCGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCAGCACTTGCAGCatattcg                                                                                                                                        <  2:1084909/150‑4 (MQ=255)
              cGTTTCGCCACGTATCCAAAACGTCGCGGCATCAAGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCAGCACTTGCAGCatattcgc                                                                                                                                       <  1:2241490/150‑5 (MQ=255)
              cGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATgggg                                                                                                                                       >  2:1546218/1‑150 (MQ=255)
                 ttCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTg                                                                                                                                    >  1:803961/1‑150 (MQ=255)
                   cGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCa                                                                                                                                  <  1:1135255/150‑1 (MQ=255)
                   cGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCa                                                                                                                                  <  2:888567/150‑1 (MQ=255)
                   cGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCa                                                                                                                                  <  1:643551/150‑1 (MQ=255)
                     aTACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGATATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATg                                                                                                                                <  2:2072622/150‑1 (MQ=255)
                      tACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGt                                                                                                                               >  1:1752948/1‑150 (MQ=255)
                           aTCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGcc                                                                                                                          <  2:1689427/150‑1 (MQ=255)
                             ccAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAg                                                                                                                        >  1:58570/1‑150 (MQ=255)
                             ccAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAg                                                                                                                        <  2:1814860/150‑1 (MQ=255)
                              cAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGt                                                                                                                       >  1:1727318/1‑150 (MQ=255)
                              cAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGt                                                                                                                       >  1:1097751/1‑150 (MQ=255)
                              cAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGt                                                                                                                       <  1:855590/150‑1 (MQ=255)
                              cAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGt                                                                                                                       >  2:2090944/1‑150 (MQ=255)
                              cAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGt                                                                                                                       >  2:1344484/1‑150 (MQ=255)
                                 aaCGTCGCGGCATCACGTGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGa                                                                                                                    <  1:1413986/150‑1 (MQ=255)
                                 aaCGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGa                                                                                                                    <  1:1074788/150‑1 (MQ=255)
                                        cgGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTa                                                                                                             >  1:732696/1‑150 (MQ=255)
                                        cgGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTa                                                                                                             >  2:2312343/1‑150 (MQ=255)
                                         ggCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTAc                                                                                                            <  1:980693/150‑1 (MQ=255)
                                         ggCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTAc                                                                                                            <  2:105387/150‑1 (MQ=255)
                                          gCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTATTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACg                                                                                                           <  1:2467270/150‑1 (MQ=255)
                                               ccgcgcgTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAg                                                                                                      <  2:976030/149‑1 (MQ=255)
                                               aCGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAg                                                                                                      >  1:1177445/1‑150 (MQ=255)
                                                          agagATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATgg                                                                                           <  2:1187905/150‑1 (MQ=255)
                                                          agagATGCTGGAGTCCGGCATTAACGTCTGCTTTGGGCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATgg                                                                                           <  2:2218224/150‑1 (MQ=255)
                                                            agaTGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGcc                                                                                         >  1:2247213/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)
                                                                 cTGGAGTCCGTCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTAGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAAt                                                                                    <  2:2005787/150‑1 (MQ=255)
                                                                   ggAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAAtt                                                                                   >  1:1847163/1‑149 (MQ=255)
                                                                   ggAGTCCGGCATTAACGGCTTCTTTGTTCACGTTGATGTCTTCTATCGGTGGTATCGGCGGGGAACGTCGAATATGGTGCAAGTGCTGGATATGTGGCTGCATGTTTGGCAGTTGATGGTCTACGGGCAGATTAACGATGGCCTGAAttt                                                                                  <  2:825864/150‑1 (MQ=255)
                                                                    gAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTa                                                                                 <  1:1238892/150‑1 (MQ=255)
                                                                      gTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAAt                                                                               <  1:2312343/150‑1 (MQ=255)
                                                                       gccGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATc                                                                              <  1:1434985/149‑1 (MQ=255)
                                                                        ccGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATca                                                                             >  1:903484/1‑150 (MQ=255)
                                                                          ggCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATcacc                                                                           >  1:954088/1‑150 (MQ=255)
                                                                           gCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATcaccc                                                                          >  1:1061332/1‑150 (MQ=255)
                                                                           gCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATcaccc                                                                          <  1:498663/150‑1 (MQ=255)
                                                                                aaCGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTCATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCcaccac                                                                     >  1:407569/1‑150 (MQ=255)
                                                                                  cGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAg                                                                   >  1:1355619/1‑150 (MQ=255)
                                                                                  cGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAg                                                                   >  2:129239/1‑150 (MQ=255)
                                                                                     cTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCTGATTAACGATGGCCTGAATTTAATCACCCACCACAgcgc                                                                <  2:1030898/150‑1 (MQ=255)
                                                                                      tGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCa                                                               <  2:1622051/150‑1 (MQ=255)
                                                                                       gCTTTGTTCACGATGATGTGTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGTATGTTTGCCAGTTGATGGGCTACTGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCaa                                                              <  2:2339518/150‑1 (MQ=255)
                                                                                       gCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCaa                                                              >  2:186395/1‑150 (MQ=255)
                                                                                         tttGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAgg                                                            <  2:1412060/150‑1 (MQ=255)
                                                                                               cACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCTAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTg                                                      <  2:1398413/150‑1 (MQ=255)
                                                                                               cACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTg                                                      <  1:2427789/150‑1 (MQ=255)
                                                                                                  gatgatGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAAt                                                   <  1:740697/150‑1 (MQ=255)
                                                                                                   atgatgTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAAtt                                                  <  2:1713624/150‑1 (MQ=255)
                                                                                                    tgatgTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAAttt                                                 <  2:1779795/150‑1 (MQ=255)
                                                                                                       tgTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCa                                              <  2:1962066/150‑1 (MQ=255)
                                                                                                       tgTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCa                                              >  2:1089454/1‑150 (MQ=255)
                                                                                                        gTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACCGCGCAAGGACGTTGAATTTGCAg                                             >  2:2181977/1‑150 (MQ=255)
                                                                                                        gTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAg                                             <  1:1344484/150‑1 (MQ=255)
                                                                                                          cTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGa                                           >  2:334138/1‑150 (MQ=255)
                                                                                                          cTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGa                                           >  1:193296/1‑150 (MQ=255)
                                                                                                          cTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGa                                           >  2:657359/1‑150 (MQ=255)
                                                                                                          cTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGa                                           >  2:1287722/1‑150 (MQ=255)
                                                                                                            tCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGAtt                                         >  2:1104501/1‑150 (MQ=255)
                                                                                                            tCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGAtt                                         >  2:1189028/1‑150 (MQ=255)
                                                                                                                  cGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCa                                   >  1:1946677/1‑150 (MQ=255)
                                                                                                                  cGTGGTAACCGCTCGGAACGGCGAATATGCTGCAAGAGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCAAAGCGCAAGCACGTTGAATTTGCAGGATTACggaa                                   >  2:1907311/1‑148 (MQ=255)
                                                                                                                   gTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCAt                                  <  1:1546218/150‑1 (MQ=255)
                                                                                                                        aTCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTgccg                             <  2:1587092/150‑1 (MQ=255)
                                                                                                                            gCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCgg                         >  2:2095544/1‑150 (MQ=255)
                                                                                                                              tGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGaa                       >  1:2409309/1‑150 (MQ=255)
                                                                                                                                ggAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAAc                     >  2:1760303/1‑150 (MQ=255)
                                                                                                                                ggAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAAc                     >  2:1637039/1‑150 (MQ=255)
                                                                                                                                    cGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAgcg                 <  1:2186210/150‑1 (MQ=255)
                                                                                                                                     ggCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAgcgc                <  2:2209304/150‑1 (MQ=255)
                                                                                                                                     gccccatatGCAGCACTTGCAGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAgcgc                >  1:1664691/6‑150 (MQ=255)
                                                                                                                                      gCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAGCGcc               >  2:2412213/1‑150 (MQ=255)
                                                                                                                                      ccccatatGCAGCACTTGCAGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAGCGcc               >  1:678338/5‑150 (MQ=255)
                                                                                                                                      ccccatatGCAGCACTTGCAGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTAAGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAGCGcc               >  2:1399578/5‑150 (MQ=255)
                                                                                                                                         catatGCAGCACTTGCAGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAGCGCCAAc            >  1:1139480/2‑150 (MQ=255)
                                                                                                                                         aaTATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAGCGCCAAc            >  1:1146984/1‑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)
                                                                                                                                                tgcAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAGCGCCAACCTGATTa     >  2:238868/1‑150 (MQ=255)
                                                                                                                                                tgcAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAGCGCCAACCTGATTa     <  1:1108373/150‑1 (MQ=255)
                                                                                                                                                tgcAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCGACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGGTGAATTTGAAGGATTACGGCATTGCCGCCGGAAACAGCGCCAACCTGATTa     <  2:760105/150‑1 (MQ=255)
                                                                                                                                                   aaGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAGCGCCAACCTGATTATcc  <  2:2247213/150‑1 (MQ=255)
                                                                                                                                                   aaGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAGCGCCAACCTGATTATcc  <  2:942493/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTCATCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAGCGCCAACCTGATTATCC  >  NC_000913/357012‑357308

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: