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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,032,3360GC10.0% 404.5 / 28.1 119intergenic (‑19/+48)dcm/yedJDNA cytosine methyltransferase/putative HD superfamily phosphohydrolase
Reads supporting (aligned to +/- strand):  ref base G (59/48);  new base C (6/6);  total (65/54)
Fisher's exact test for biased strand distribution p-value = 7.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

Only 100 of 121 total aligned reads displayed.

CACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAAT  >  NC_000913/2032194‑2032484
                                                                                                                                              |                                                                                                                                                    
cACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATgg                                                                                                                                                >  3:131614/1‑149 (MQ=255)
 aCCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGc                                                                                                                                               <  4:131614/149‑1 (MQ=255)
  ccTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGcc                                                                                                                                              >  4:725432/1‑149 (MQ=255)
   cTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCg                                                                                                                                             <  2:165399/149‑1 (MQ=255)
   cTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCg                                                                                                                                             <  3:287030/149‑1 (MQ=255)
     acacCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGAc                                                                                                                                           >  2:46262/1‑149 (MQ=255)
     acacCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCggt                                                                                                                                           <  1:1043756/149‑3 (MQ=255)
      cacCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACa                                                                                                                                          <  3:529332/149‑1 (MQ=255)
        ccATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAgg                                                                                                                                        <  2:187508/149‑1 (MQ=255)
         cATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGt                                                                                                                                       >  3:450905/1‑149 (MQ=255)
         cATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGt                                                                                                                                       >  2:920065/1‑149 (MQ=255)
         cATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGt                                                                                                                                       <  2:85476/149‑1 (MQ=255)
         cATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGt                                                                                                                                       >  1:791973/1‑149 (MQ=255)
          aTTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGtt                                                                                                                                      >  1:205903/1‑149 (MQ=255)
          aTTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGtt                                                                                                                                      >  1:62932/1‑149 (MQ=255)
          aTTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGtt                                                                                                                                      <  1:920065/149‑1 (MQ=255)
          aTTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGtt                                                                                                                                      <  1:970472/149‑1 (MQ=255)
              aGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACa                                                                                                                                  >  1:763673/1‑149 (MQ=255)
               gCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACAt                                                                                                                                 <  4:57058/149‑1 (MQ=255)
                 tGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATaa                                                                                                                               <  1:501895/149‑1 (MQ=255)
                  gcgcCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAAt                                                                                                                              <  1:421732/149‑1 (MQ=255)
                      caccaACGTTTTAACATCATAAATTTTCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGTATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACa                                                                                                                          <  2:65262/149‑1 (MQ=255)
                       accaACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAg                                                                                                                         >  4:952304/1‑149 (MQ=255)
                          aaCGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGAc                                                                                                                                           >  1:604970/1‑128 (MQ=255)
                          aaCGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGAc                                                                                                                                           <  2:604970/128‑1 (MQ=255)
                          aaCGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCa                                                                                                                      >  1:352842/1‑149 (MQ=255)
                             gTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCAttt                                                                                                                   <  3:497014/149‑1 (MQ=255)
                               tttAACATCATAAATTTGCAGCAGTTTCTCTAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAg                                                                                                                 <  3:879124/149‑1 (MQ=255)
                                   aCATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGACATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGcc                                                                                                             >  2:967948/1‑149 (MQ=255)
                                    catcatAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCa                                                                                                            >  1:817006/1‑149 (MQ=255)
                                     atcatAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCggt                                                                                                                                           >  4:105065/1‑115 (MQ=255)
                                      tcatAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAAt                                                                                                                              >  2:334001/1‑129 (MQ=255)
                                      tcatAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCAcg                                                                                                          >  4:524129/1‑149 (MQ=255)
                                        atAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCAcgcg                                                                                                        >  2:120462/1‑149 (MQ=255)
                                         tAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACgcgc                                                                                                       <  4:856056/149‑1 (MQ=255)
                                          aaaTTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCt                                                                                                      <  2:205903/149‑1 (MQ=255)
                                          aaaTTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCt                                                                                                      >  4:474798/1‑149 (MQ=255)
                                                gcagcaGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTcagcca                                                                                                <  2:973171/149‑1 (MQ=255)
                                                   gcagTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCg                                                                                             >  3:404789/1‑149 (MQ=255)
                                                   gcagTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCg                                                                                             >  1:1079477/1‑149 (MQ=255)
                                                   gcagTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGCCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCg                                                                                             <  1:2023/149‑1 (MQ=255)
                                                    cagTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCgg                                                                                            <  2:196825/149‑1 (MQ=255)
                                                     agTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGa                                                                                           <  2:791973/149‑1 (MQ=255)
                                                      gTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGAt                                                                                          <  2:817006/149‑1 (MQ=255)
                                                       tttCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATg                                                                                         >  1:1041366/1‑149 (MQ=255)
                                                           tcCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACGTAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGAATGaaaa                                                                                     >  3:931436/1‑149 (MQ=255)
                                                           tcCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAgg                                                                                                                >  1:954424/1‑122 (MQ=255)
                                                           tcCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAgg                                                                                                                <  2:954424/122‑1 (MQ=255)
                                                           tcCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGaaaa                                                                                     >  4:212790/1‑149 (MQ=255)
                                                           tcCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGaaaa                                                                                     >  3:474631/1‑149 (MQ=255)
                                                           tcCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGaaaa                                                                                     >  3:25360/1‑149 (MQ=255)
                                                           tcCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGaaaa                                                                                     >  2:395431/1‑149 (MQ=255)
                                                              aGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGc                                                                                  >  3:773845/1‑149 (MQ=255)
                                                                cATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCAt                                                                                >  2:681783/1‑149 (MQ=255)
                                                                cATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAAGTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCAt                                                                                <  4:474631/149‑1 (MQ=255)
                                                                 aTTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATc                                                                               <  2:610081/149‑1 (MQ=255)
                                                                             gCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTtgtg                                                                   <  4:549311/149‑1 (MQ=255)
                                                                                ttCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATc                                                                <  3:701218/149‑1 (MQ=255)
                                                                                  cccGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGa                                                              >  2:250222/1‑149 (MQ=255)
                                                                                   ccGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGac                                                             >  2:412432/1‑149 (MQ=255)
                                                                                      gTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACa                                                                                                                          <  1:520422/85‑1 (MQ=255)
                                                                                      gTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACa                                                                                                                          >  2:520422/1‑85 (MQ=255)
                                                                                           gTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGAt                                                                 <  3:825074/137‑1 (MQ=255)
                                                                                           gTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGAt                                                                 >  4:825074/1‑137 (MQ=255)
                                                                                           gTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCAt                                                     >  4:122747/1‑149 (MQ=255)
                                                                                                 aTCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTcagccagc                                                                                              <  1:856897/102‑1 (MQ=255)
                                                                                                 aTCGGTTACTAATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTcagccagc                                                                                              >  2:856897/1‑102 (MQ=255)
                                                                                                  tCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTcccc                                              >  3:114852/1‑149 (MQ=255)
                                                                                                    ggTTACTGATATATTTTCCTGCATAGATTTCACGGGCGATGTGAGATGGCGGAGAGGTTAACATAATTAGAGATATTTAGGGATAGGCGTTCAGTTAGCGGATGAAAACGTATATATTATCTTGTGATTGAGATGTTAATtctctcgtg                                            >  2:56763/1‑144 (MQ=255)
                                                                                                     gTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGc                                           >  3:761710/1‑149 (MQ=255)
                                                                                                          tGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGtt                                      >  3:151921/1‑149 (MQ=255)
                                                                                                           gATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTc                                     <  3:443374/149‑1 (MQ=255)
                                                                                                           gATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTc                                     <  2:538504/149‑1 (MQ=255)
                                                                                                            atatatTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCAc                                                                                                           >  1:377854/1‑78 (MQ=255)
                                                                                                            atatatTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCg                                    >  2:45594/1‑149 (MQ=255)
                                                                                                               tatTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCAtt                                                    >  3:1054120/1‑130 (MQ=255)
                                                                                                               tatTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCa                                 >  3:884311/1‑149 (MQ=255)
                                                                                                               tatTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCa                                 >  2:410427/1‑149 (MQ=255)
                                                                                                                    tCCTGCATAGATTTCACCGGCCAGCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAg                            <  1:982971/149‑1 (MQ=255)
                                                                                                                        gCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTc                        >  4:712110/1‑149 (MQ=255)
                                                                                                                        gCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTc                        >  1:829564/1‑149 (MQ=255)
                                                                                                                         cATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCg                       >  2:1034362/1‑149 (MQ=255)
                                                                                                                          aTAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGa                             >  2:514279/1‑142 (MQ=255)
                                                                                                                          aTAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGa                             <  1:514279/142‑1 (MQ=255)
                                                                                                                          aTAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGc                      >  4:202667/1‑149 (MQ=255)
                                                                                                                            aGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAg                            >  4:120386/1‑141 (MQ=255)
                                                                                                                            aGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGAGGATGAAAACGCATCTATCCCCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAg                            <  3:120386/141‑1 (MQ=255)
                                                                                                                               tttCACCGGCCATCTGAGAAGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCCGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATaa                 >  2:1234/1‑149 (MQ=255)
                                                                                                                                  tgtcGGCCATCTCAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACt              >  4:463156/4‑149 (MQ=255)
                                                                                                                                   aCCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTc             <  2:1040316/149‑1 (MQ=255)
                                                                                                                                   aCCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTc             <  1:424292/149‑1 (MQ=255)
                                                                                                                                   aCCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTc             <  1:412432/149‑1 (MQ=255)
                                                                                                                                     cGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGa           <  2:62932/149‑1 (MQ=255)
                                                                                                                                     cGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGa           >  4:450873/1‑149 (MQ=255)
                                                                                                                                      ggCCATCTCAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATtc                                                   >  4:818875/1‑108 (MQ=255)
                                                                                                                                      ggCCATCTCAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATtc                                                   <  3:818875/108‑1 (MQ=255)
                                                                                                                                      ggCCATCTCAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGAc          >  3:1007386/1‑149 (MQ=255)
                                                                                                                                      ggCCATCTCAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGAc          >  2:508543/1‑149 (MQ=255)
                                                                                                                                            cTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAATTCGGCACTGAGCTTCGCCATAAACTCGACTaaaaa    <  4:677326/149‑1 (MQ=255)
                                                                                                                                              gagaTGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAAt  <  3:375192/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
CACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAAT  >  NC_000913/2032194‑2032484

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