Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,032,336 G→C 8.1% intergenic (‑19/+48) dcm ← / ← yedJ DNA‑cytosine methyltransferase/putative HD superfamily phosphohydrolase YedJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,032,3360GC 8.1% 285.0 / 13.1 86intergenic (‑19/+48)dcm/yedJDNA‑cytosine methyltransferase/putative HD superfamily phosphohydrolase YedJ
Reads supporting (aligned to +/- strand):  ref base G (40/39);  new base C (2/5);  total (42/44)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

TGATTCTCACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAAT  >  NC_000913/2032187‑2032484
                                                                                                                                                     |                                                                                                                                                    
tGATTCTCACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTg                                                                                                                                                      <  1:1286323/150‑1 (MQ=255)
   ttCTCACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTgaga                                                                                                                                                   >  2:779990/1‑150 (MQ=255)
    tctcACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGAt                                                                                                                                                  <  1:1398555/150‑1 (MQ=255)
     ctcACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATg                                                                                                                                                 >  1:1487266/1‑150 (MQ=255)
     ctcACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATg                                                                                                                                                 <  1:1390226/150‑1 (MQ=255)
       cACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCCGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGGTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGc                                                                                                                                               >  2:1119398/1‑150 (MQ=255)
       cACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGc                                                                                                                                               >  1:200599/1‑150 (MQ=255)
       cACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGc                                                                                                                                               >  1:207551/1‑150 (MQ=255)
       cACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGc                                                                                                                                               >  2:941073/1‑150 (MQ=255)
       cACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCAGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGc                                                                                                                                               >  2:185919/1‑150 (MQ=255)
         ccTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCg                                                                                                                                             >  1:559489/1‑150 (MQ=255)
           tACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGAc                                                                                                                                           <  2:783106/150‑1 (MQ=255)
           tACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCggt                                                                                                                                           <  2:721719/150‑3 (MQ=255)
           tACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCggt                                                                                                                                           <  2:381589/150‑3 (MQ=255)
            acacCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCggtg                                                                                                                                          <  1:111626/150‑4 (MQ=255)
             cacCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCggtga                                                                                                                                         <  1:400879/150‑5 (MQ=255)
               ccATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGt                                                                                                                                       >  1:232612/1‑150 (MQ=255)
               ccATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGt                                                                                                                                       >  1:1700725/1‑150 (MQ=255)
                cATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGtt                                                                                                                                      <  1:299217/150‑1 (MQ=255)
                    aaGCTGCGCCACCAACGTTTTAACACCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTg                                                                                                                                                      <  1:986235/130‑1 (MQ=255)
                    aaGCTGCGCCACCAACGTTTTAACAACATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTg                                                                                                                                                      >  2:986235/1‑130 (MQ=255)
                         gcgcCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAAtt                                                                                                                             <  2:1846314/150‑1 (MQ=255)
                            ccaccaACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACa                                                                                                                          >  2:1132507/1‑150 (MQ=255)
                             caccaACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAg                                                                                                                         >  2:1341808/1‑150 (MQ=255)
                             caccaACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAg                                                                                                                         <  2:559489/150‑1 (MQ=255)
                              accaACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGc                                                                                                                        <  2:201132/150‑1 (MQ=255)
                                caACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCa                                                                                                                      >  1:920225/1‑150 (MQ=255)
                                    gTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCGTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGCCAGGTTAACATTATTACAGCCATTTa                                                                                                                  >  2:1434864/1‑150 (MQ=255)
                                        tAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGc                                                                                                              >  1:954580/1‑150 (MQ=255)
                                          aCATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAAGCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCa                                                                                                            <  2:1295791/150‑1 (MQ=255)
                                           catcatAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCACTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCAc                                                                                                           >  1:1335626/1‑150 (MQ=255)
                                             tcatAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCGCCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCAcgc                                                                                                         >  1:559455/1‑150 (MQ=255)
                                              catAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCAcgcg                                                                                                        >  1:1432122/1‑150 (MQ=255)
                                                        cagcagTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTcagccagc                                                                                              >  1:390519/1‑150 (MQ=255)
                                                        cagcagTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTcagccagc                                                                                              >  1:1949884/1‑150 (MQ=255)
                                                         agcagTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCg                                                                                             <  1:555635/150‑1 (MQ=255)
                                                          gcagTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCgg                                                                                            <  2:207551/150‑1 (MQ=255)
                                                          gcagTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCgg                                                                                            >  1:348085/1‑150 (MQ=255)
                                                          gcagTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCgg                                                                                            <  2:200599/150‑1 (MQ=255)
                                                            agTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGAt                                                                                          <  1:1126380/150‑1 (MQ=255)
                                                                tctcCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGaaa                                                                                      <  1:746084/150‑1 (MQ=255)
                                                                   ccAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACg                                                                                   <  1:779990/150‑1 (MQ=255)
                                                                          tGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTAt                                                                            >  2:1805640/1‑150 (MQ=255)
                                                                           gCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATc                                                                           >  1:1142530/1‑150 (MQ=255)
                                                                             cTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCAc                                                                         >  2:275572/1‑150 (MQ=255)
                                                                              tgtgCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCAcc                                                                        <  1:1132507/150‑1 (MQ=255)
                                                                                 gCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAggg                                                                                                               >  2:1481873/1‑108 (MQ=255)
                                                                                 gCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAggg                                                                                                               <  1:1481873/108‑1 (MQ=255)
                                                                                 gCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCACTCAGCCAGCGGATGAAAACGCATCTATCACCTtg                                                                     <  1:1129292/150‑1 (MQ=255)
                                                                                  cGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTtgt                                                                    <  1:1341808/150‑1 (MQ=255)
                                                                                     cATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCGTTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGAt                                                                 >  1:1662269/1‑150 (MQ=255)
                                                                                     cATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGAt                                                                 <  2:1487266/150‑1 (MQ=255)
                                                                                      aTTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATc                                                                <  1:275572/150‑1 (MQ=255)
                                                                                       ttCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCTGCCATCTGAGATGGTCTACATGTTAACATAATTACAGCGATTTAGGGCCACTCTCTCAGCCAGCTGATGAATATGGATGTATCAGCTTGTGTTCg                                                               <  2:678927/150‑1 (MQ=255)
                                                                                        tCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGa                                                              <  1:1833253/150‑1 (MQ=255)
                                                                                          ccGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGaca                                                            <  2:1142530/150‑1 (MQ=255)
                                                                                          ccGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGaca                                                            >  2:564555/1‑150 (MQ=255)
                                                                                            ggTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACAcc                                                          <  1:264208/150‑1 (MQ=255)
                                                                                                gtGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCa                                                      <  2:1962541/149‑1 (MQ=255)
                                                                                                   tATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATtc                                                   <  1:683868/150‑1 (MQ=255)
                                                                                                      gAATCGGTTACTGATATATTTTCCTGGATAGATTTCACCGGCCAGCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTcc                                                <  2:1452858/150‑1 (MQ=255)
                                                                                                        aTCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATtc                                                   >  2:537774/1‑145 (MQ=255)
                                                                                                        aTCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATtc                                                   <  1:537774/145‑1 (MQ=255)
                                                                                                          cGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCg                                            <  1:1872476/150‑1 (MQ=255)
                                                                                                              tACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAg                                        >  2:731565/1‑150 (MQ=255)
                                                                                                              tACGGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAg                                        <  1:595319/150‑1 (MQ=255)
                                                                                                                 tGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTc                                     <  2:348085/150‑1 (MQ=255)
                                                                                                                  gATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCg                                    >  1:702480/1‑150 (MQ=255)
                                                                                                                    tatatTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGc                                  <  2:1335626/150‑1 (MQ=255)
                                                                                                                     atatTTTCCTGCATAGATTTCGCCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCa                                 <  2:559455/150‑1 (MQ=255)
                                                                                                                     atatTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGa                                                              <  1:801968/121‑1 (MQ=255)
                                                                                                                     atatTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGa                                                              >  2:801968/1‑121 (MQ=255)
                                                                                                                     atatTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCa                                 <  2:1949884/150‑1 (MQ=255)
                                                                                                                      tatTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCAc                                >  1:1265770/1‑150 (MQ=255)
                                                                                                                         tttCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGGTGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGa                             >  2:1838653/1‑150 (MQ=255)
                                                                                                                         tttCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGa                             <  1:185919/150‑1 (MQ=255)
                                                                                                                          ttCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAg                            <  1:155939/150‑1 (MQ=255)
                                                                                                                           tCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGc                           <  2:28610/150‑1 (MQ=255)
                                                                                                                                cATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGc                      >  2:50340/1‑150 (MQ=255)
                                                                                                                                 aTAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGcc                     >  1:1395495/1‑150 (MQ=255)
                                                                                                                                  tAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGGGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGACa                    >  2:1646149/1‑150 (MQ=255)
                                                                                                                                       ttCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAAc               <  1:161763/150‑1 (MQ=255)
                                                                                                                                          gtcGGCCATCTCAGATGGCCGGCAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCg            >  1:328076/3‑150 (MQ=255)
                                                                                                                                          gtcGGCCATCTCAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCa                    >  2:1579730/3‑142 (MQ=255)
                                                                                                                                          gtcGGCCATCTCAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCa                    <  1:1579730/140‑1 (MQ=255)
                                                                                                                                          aCCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCg            >  1:1867953/1‑150 (MQ=255)
                                                                                                                                             ggCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACt         >  1:816572/1‑150 (MQ=255)
                                                                                                                                                   cTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTaaaaaa   >  1:795782/1‑150 (MQ=255)
                                                                                                                                                    tGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAAt  >  2:1095184/1‑150 (MQ=255)
                                                                                                                                                    tGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAAt  >  1:975963/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGATTCTCACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAAT  >  NC_000913/2032187‑2032484

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: ‑