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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,1650TA 5.9% 273.4 / 5.5 85L332Q (CTG→CAG) codAcytosine/isoguanine deaminase
Reads supporting (aligned to +/- strand):  ref base T (40/40);  new base A (2/3);  total (42/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.80e-01

CTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAGCGCCAACCTGATTATCCTGCCGG  >  NC_000913/357017‑357314
                                                                                                                                                    |                                                                                                                                                     
cTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGtgctg                                                                                                                                                      <  1:1615485/150‑1 (MQ=255)
cTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGtgctg                                                                                                                                                      <  2:1303357/150‑1 (MQ=255)
cTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGtgctg                                                                                                                                                      <  1:785657/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:439088/150‑5 (MQ=255)
         cGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCAGCACTTGCAGCatattcgc                                                                                                                                             <  1:975928/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                                                                                                                                   >  2:1632639/1‑150 (MQ=255)
                   cGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGttt                                                                                                                                   >  1:1394027/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:1394027/150‑1 (MQ=255)
                                      cATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACgg                                                                                                                >  2:1003861/1‑150 (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                                                                       >  1:1048487/1‑150 (MQ=255)
                                                                               tCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAgcg                                                                       >  2:1324230/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                                                                                    <  1:670830/134‑1 (MQ=255)
                                                                                  gCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATc                                                                                    >  2:670830/1‑134 (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                                                    <  2:1248782/150‑1 (MQ=255)
                                                                                                  tgTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCa                                                    >  1:121806/1‑150 (MQ=255)
                                                                                                   gTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAg                                                   <  1:86713/150‑1 (MQ=255)
                                                                                                   gTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAg                                                   <  1:1518518/150‑1 (MQ=255)
                                                                                                    tCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAgg                                                  >  1:643819/1‑150 (MQ=255)
                                                                                                    tCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAgg                                                  >  1:1221714/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            <  1:288907/150‑1 (MQ=255)
                                                                                                                                          ctgcAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAGCGCCAACCTGAtt            >  2:288907/1‑150 (MQ=255)
                                                                                                                                              aaGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAGCGCCAACCTGATTATcc        <  1:289169/150‑1 (MQ=255)
                                                                                                                                                    tgcATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGCCGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAGCGCCAACCTGATTATCCTGCCgg  >  2:1634041/1‑150 (MQ=255)
                                                                                                                                                    tgcATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGc                                          >  1:390081/1‑110 (MQ=255)
                                                                                                                                                    tgcATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGc                                          <  2:390081/110‑1 (MQ=255)
                                                                                                                                                    tgcATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAGCGCCAACCTGATTATCCTGCCgg  >  1:1196927/1‑150 (MQ=255)
                                                                                                                                                    tgcATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAGCGCCAACCTGATTATCCTGCCgg  >  2:1100944/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAGCGCCAACCTGATTATCCTGCCGG  >  NC_000913/357017‑357314

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: