Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 357,160 A→C 5.4% 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.4% 288.9 / 7.8 93Q330H (CAA→CACcodAcytosine/isoguanine deaminase
Reads supporting (aligned to +/- strand):  ref base A (43/45);  new base C (2/3);  total (45/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

TTCATCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAGCGCCAACCTGATTATCC  >  NC_000913/357012‑357308
                                                                                                                                                    |                                                                                                                                                    
ttCATCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAg                                                                                                                                                     >  2:594559/1‑150 (MQ=255)
ttCATCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAg                                                                                                                                                     <  1:33354/150‑1 (MQ=255)
ttCATCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAg                                                                                                                                                     <  1:1389319/150‑1 (MQ=255)
  cATCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGtg                                                                                                                                                   >  2:828090/1‑150 (MQ=255)
  cATCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGtg                                                                                                                                                   >  2:528799/1‑150 (MQ=255)
     cTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGtgctg                                                                                                                                                <  2:1303357/150‑1 (MQ=255)
     cTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGtgctg                                                                                                                                                <  1:785657/150‑1 (MQ=255)
     cTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGtgctg                                                                                                                                                <  1:1615485/150‑1 (MQ=255)
      tGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTTGGAACGGCGAATATGCTGCAAGtgctgc                                                                                                                                               <  2:222561/150‑1 (MQ=255)
       gCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCa                                                                                                                                              >  2:590629/1‑150 (MQ=255)
        cAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCat                                                                                                                                             <  1:528799/150‑1 (MQ=255)
        cAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCat                                                                                                                                             >  1:1548009/1‑150 (MQ=255)
        cAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCat                                                                                                                                             <  1:1505351/150‑1 (MQ=255)
            gACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATgg                                                                                                                                         >  2:646602/1‑150 (MQ=255)
              cGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATgggg                                                                                                                                       <  1:594559/150‑1 (MQ=255)
              cGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCAGCACTTGCAGCatattcgc                                                                                                                                       <  1:975928/150‑5 (MQ=255)
              cGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCAGCACTTGCAGCatattcgc                                                                                                                                       <  1:439088/150‑5 (MQ=255)
               gTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAAAGGCGAATATGCAGCACTTGCAGCatattcgcc                                                                                                                                      <  1:87343/150‑6 (MQ=255)
                  tCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGc                                                                                                                                   <  1:431508/150‑1 (MQ=255)
                   cGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCa                                                                                                                                  <  2:696488/150‑1 (MQ=255)
                   cGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCa                                                                                                                                  <  2:994604/150‑1 (MQ=255)
                      tACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGt                                                                                                                               >  2:73618/1‑150 (MQ=255)
                        cGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGttt                                                                                                                             >  1:1394027/1‑150 (MQ=255)
                        cGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGttt                                                                                                                             >  2:1632639/1‑150 (MQ=255)
                         gTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTg                                                                                                                            <  2:1503410/150‑1 (MQ=255)
                          tATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGc                                                                                                                           >  2:684286/1‑150 (MQ=255)
                             ccAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAg                                                                                                                        >  2:639067/1‑150 (MQ=255)
                              cAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTACTGCATATGGGGCTGCATGTTTGCCAGt                                                                                                                       >  1:189006/1‑150 (MQ=255)
                               aaaaCGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGtt                                                                                                                      >  2:80588/1‑150 (MQ=255)
                               aaaaCGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCAAGTTTGCCAGtt                                                                                                                      >  2:949551/1‑150 (MQ=255)
                                    gTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATgg                                                                                                                 <  2:1554187/150‑1 (MQ=255)
                                      cgcgGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGc                                                                                                               <  2:104410/150‑1 (MQ=255)
                                       gcgGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCt                                                                                                              >  1:1282465/1‑150 (MQ=255)
                                          gCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACg                                                                                                           <  1:73618/150‑1 (MQ=255)
                                           cATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACgg                                                                                                          >  2:1003861/1‑150 (MQ=255)
                                           cATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACgg                                                                                                          <  2:1394027/150‑1 (MQ=255)
                                              cACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCa                                                                                                       <  1:639067/150‑1 (MQ=255)
                                               aCGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAg                                                                                                      <  1:1541346/150‑1 (MQ=255)
                                                cgcgcgTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGa                                                                                                     <  1:590629/150‑1 (MQ=255)
                                                 gcgcgTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGAt                                                                                                    <  1:80588/150‑1 (MQ=255)
                                                     gTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTTTCCGCTGTGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGTCTACGGTCATATTAAc                                                                                                <  2:420335/150‑1 (MQ=255)
                                                     gTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAAc                                                                                                <  1:929566/150‑1 (MQ=255)
                                                     gTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTCCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAAc                                                                                                >  1:559399/1‑150 (MQ=255)
                                                        aaaGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGAt                                                                                             <  1:646602/150‑1 (MQ=255)
                                                        aaaGAGATGCTGGAGTAAGGCATTACCGTCTGCTTTGGTCACGAAGATGTCTTAGATCAGTGGTCTCCGCTGGGAACGGCAAATATGCTGCAAGTGCTGCATAAGCGGCTGCATGTTTGACAGTTGATGGGCTAAAGGCAGATTAACCAt                                                                                             >  2:1096870/1‑150 (MQ=255)
                                                         aaGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATg                                                                                            <  2:1373624/150‑1 (MQ=255)
                                                         aaGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATg                                                                                            <  2:1574312/150‑1 (MQ=255)
                                                            cgaTGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGcc                                                                                         <  1:654783/149‑1 (MQ=255)
                                                             gaTGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCt                                                                                        >  2:713878/1‑150 (MQ=255)
                                                              aTGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTg                                                                                       >  2:504861/1‑150 (MQ=255)
                                                                      gTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAAt                                                                               <  1:1632639/150‑1 (MQ=255)
                                                                        ccGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAAttt                                                                                  <  1:1092501/145‑1 (MQ=255)
                                                                        ccGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCCCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAAttt                                                                                  >  2:1092501/1‑145 (MQ=255)
                                                                            cATTAACGTCTGCATTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCCAATATGCTGCAAGTTCTGCATATTGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGAAGGCCTGAATTTAATcaccca                                                                         >  2:1331495/1‑150 (MQ=255)
                                                                                  cGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAg                                                                   >  1:97678/1‑150 (MQ=255)
                                                                                   gTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCAACCACAgc                                                                  <  1:520798/150‑1 (MQ=255)
                                                                                    tCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAgcg                                                                 >  2:1324230/1‑150 (MQ=255)
                                                                                    tCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAgcg                                                                 >  1:1048487/1‑150 (MQ=255)
                                                                                     cTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAgcgc                                                                <  1:457343/150‑1 (MQ=255)
                                                                                       gCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATcac                                                                            <  1:1063961/136‑1 (MQ=255)
                                                                                       gCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATcac                                                                            >  2:1063961/1‑136 (MQ=255)
                                                                                       gCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATc                                                                              >  2:670830/1‑134 (MQ=255)
                                                                                       gCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATc                                                                              <  1:670830/134‑1 (MQ=255)
                                                                                          ttGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGa                                                           >  2:830196/1‑150 (MQ=255)
                                                                                           tGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGAc                                                          >  1:1518086/1‑150 (MQ=255)
                                                                                            ggTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACg                                                         >  2:338462/1‑150 (MQ=255)
                                                                                             gTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGt                                                        >  2:1279562/1‑150 (MQ=255)
                                                                                               cACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGtta                                                      >  2:886040/1‑149 (MQ=255)
                                                                                               cACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTg                                                      <  1:830196/150‑1 (MQ=255)
                                                                                               cACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTg                                                      >  2:1518518/1‑150 (MQ=255)
                                                                                                aCGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGa                                                     >  1:891042/1‑150 (MQ=255)
                                                                                                 cGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTACTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGaa                                                    <  2:189006/150‑1 (MQ=255)
                                                                                                   atgatgTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAAtt                                                  <  2:1537352/150‑1 (MQ=255)
                                                                                                    tgatgTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAAttt                                                 <  1:1312145/150‑1 (MQ=255)
                                                                                                     gatgTCTTCGATCCGTGGTATCCGCTGGGAACGTCGAATATGCTTCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATTTTCTACGGGAAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTg                                                <  2:867582/150‑1 (MQ=255)
                                                                                                       tgTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCa                                              >  1:121806/1‑150 (MQ=255)
                                                                                                       tgTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCa                                              <  2:1248782/150‑1 (MQ=255)
                                                                                                        gTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAg                                             <  1:1518518/150‑1 (MQ=255)
                                                                                                        gTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAg                                             <  1:86713/150‑1 (MQ=255)
                                                                                                         tCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAgg                                            >  1:1221714/1‑150 (MQ=255)
                                                                                                         tCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAgg                                            >  1:643819/1‑150 (MQ=255)
                                                                                                            tCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCATGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGAtt                                         >  2:1118009/1‑150 (MQ=255)
                                                                                                             cGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGTTATACTTTTATATTAATTATGGCCTGAATTTAATCAACCACGACAGAGAAAGGAAGTTTAATTTGCAGGATTa                                        <  2:787578/150‑1 (MQ=255)
                                                                                                                  cGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCa                                   <  2:994406/150‑1 (MQ=255)
                                                                                                                            gCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCgg                         >  1:404073/1‑150 (MQ=255)
                                                                                                                                  aaCGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAg                   >  1:593661/1‑150 (MQ=255)
                                                                                                                                   cagccccatatGCAGCACTTGCAGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAgc                  >  1:1351999/8‑150 (MQ=255)
                                                                                                                                     ggCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAgcgc                >  1:1259552/1‑150 (MQ=255)
                                                                                                                                     gccccatatGCAGCACTTGCAGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAgcgc                >  1:89206/6‑150 (MQ=255)
                                                                                                                                            atGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCCGAGCTCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAGCGCCAACCTg         <  2:559399/150‑1 (MQ=255)
                                                                                                                                             tgctgcAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAGCGCCAACCTGa        >  2:94380/1‑150 (MQ=255)
                                                                                                                                               ctgcAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAGCGCCAACCTGAtt      >  2:288907/1‑150 (MQ=255)
                                                                                                                                               ctgcAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAGCGCCAACCTGAtt      <  1:288907/150‑1 (MQ=255)
                                                                                                                                                   aaGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAGCGCCAACCTGATTATcc  <  1:289169/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: