Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,032,336 G→C 5.3% 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,3360GC 5.3% 802.4 / 28.1 228intergenic (‑19/+48)dcm/yedJDNA cytosine methyltransferase/putative HD superfamily phosphohydrolase
Reads supporting (aligned to +/- strand):  ref base G (107/108);  new base C (6/6);  total (113/115)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.63e-03

Only 100 of 233 total aligned reads displayed.

GATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATGGGCGAACTCTTTCTCACTTAAACGGTGCCATGCCGAGTCATTCGCCAGCGCACGTTTTAAAATTGCCGCGCTCCAGTGATTCTCACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAAC  >  NC_000913/2032071‑2032609
                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                 
gatgatGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATGGGCGAACTCTTTCTCACTTAAACGGTGCCATGCCGAGTCATTCGCCAGCGCACGTTTTAAAATTGCCGCGCTCCAGTGATTCTCACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCCTTCCCGATGCTGTATGAATCGTTTACGGATATGTTGCTCTTCATATCTTTCCCCGGCCATGTGAGATGGCCg                                                                                                                                                                                                                                                                          >  4:478215/1‑275 (MQ=255)
gatgatGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATGGGCGAACTCTTTCTCACTTAAACGGTGCCATGCCGAGTCATTCGCCAGCGCACGTTTTAAAATTGCCGCGCTCCAGTGATTCTCACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGa                                                                                                                                                                                                                                                                         <  1:220647/276‑1 (MQ=255)
gatgatGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATGGGCGAACTCTTTCTCACTTAAACGGTGCCATGCCGAGTCATTCGCCAGCGCACGTTTTAAAATTGCCGCGCTCCAGTGATTCTCACCTACACCATTAAGCTACGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGAATTGCCTGTNGGGCATTCCCGGTGCTGTATGCATCGTTTCCGTATATATTTTCCTGCATAAATTTCACAGGATACCTGAAATgggcaaaac                                                                                                                                                                                                                                                                      >  3:220647/1‑272 (MQ=255)
              ggtTTGGGTAATAACGTTTGCAGATGGGCGAACTCTTTCTCACTTAAACGGTGCCATGCCGAGTCATTCGCCAGCGCACGTTTTAAAATTGCCGCGCTCCAGTGATTCTCACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCTGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATaa                                                                                                                                                                                                                                                            >  2:294048/1‑275 (MQ=255)
                 tttggTAATAACGTTTGCCGAGTGGCGAACTCTTTCGCCCTTAAACGGTGCCATGCCGAGTCATTCGCCAGCGCACGTTTTAAAATTGCCGCGCTCCAGTGATTCTCACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAAtt                                                                                                                                                                                                                                                          <  4:534467/271‑1 (MQ=255)
                                          cGAACTCTTTCTCACTTAAACGGTGCCATGCCGAGTCATTCGCCAGCGCACGTTTTAAAATTGCCGCGCTCCAGTGATTCTCACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTca                                                                                                                                                                                                                                 >  2:245148/1‑274 (MQ=255)
                                                     tcACTTAAACGGTGCCATGCCGAGTCATTCGCCAGCGCACGTTTTAAAATTGCCGCGCTCCAGTGATTCTCACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAAc                                                                                                                                                                                                                 <  1:130151/279‑1 (MQ=255)
                                                                                                  aaaaTCTCCGCGCTCCAGTGATTCTATCCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGGGCTGTATGAATCGGTTACTGATATATTTTCCTGAATAGATTGCACCGGCCAGCTGAGATGGACGACAGGTTGAAAAAAGTACAGCCATTTAGAGCCACTCGCTCAGCCAGCGGATGGAAACGCATCTATAACCTTGTGATCGACACCTTCATTCTCCCCCGc                                                                                                                                                                        <  4:245145/275‑1 (MQ=255)
                                                                                                         ccGCGCTCCAGTGATTCTCACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGcc                                                                                                                                                              >  1:17188/1‑277 (MQ=255)
                                                                                                                      aTTCTCACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTga                                                                                                                                                                                                                                                                                  <  8:1002221/149‑1 (MQ=255)
                                                                                                                       ttCTCACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTGACTGATATGTTTTCCTGCATAGATTTCACCGGCCATCTgag                                                                                                                                                                                                                                                                                 >  8:1376233/1‑149 (MQ=255)
                                                                                                                             ccTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGtt                                                                                                                                                                                                                                                                   <  4:133884/157‑1 (MQ=255)
                                                                                                                             ccTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAAc                                                                                                                                                                                                                                                                <  2:463443/160‑1 (MQ=255)
                                                                                                                             ccTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGcc                                                                                                                                                                                                                                                                           <  7:290610/149‑1 (MQ=255)
                                                                                                                              cTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCg                                                                                                                                                                                                                                                                          >  8:768697/1‑149 (MQ=255)
                                                                                                                                acacCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGAc                                                                                                                                                                                                                                                                        >  2:699582/1‑149 (MQ=255)
                                                                                                                                acacCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGAc                                                                                                                                                                                                                                                                        <  8:1675760/149‑1 (MQ=255)
                                                                                                                                    cATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGt                                                                                                                                                                                                                                                                    <  6:399331/149‑1 (MQ=255)
                                                                                                                                    cATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGt                                                                                                                                                                                                                                                                    >  5:574022/1‑149 (MQ=255)
                                                                                                                                    cATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGt                                                                                                                                                                                                                                                                    <  4:699577/149‑1 (MQ=255)
                                                                                                                                     aTTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGtt                                                                                                                                                                                                                                                                   >  5:145980/1‑149 (MQ=255)
                                                                                                                                      ttAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTa                                                                                                                                                                                                                                                                  >  8:966405/1‑149 (MQ=255)
                                                                                                                                            tGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATaa                                                                                                                                                                                                                                                            <  6:1078337/149‑1 (MQ=255)
                                                                                                                                                  accaACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGTCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCg                                                                                                                                                                                                                                                                          >  5:1104556/1‑129 (MQ=255)
                                                                                                                                                  accaACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCg                                                                                                                                                                                                                                                                          <  6:1104556/129‑1 (MQ=255)
                                                                                                                                                    caACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACa                                                                                                                                                                                                                                                                       <  6:865717/130‑1 (MQ=255)
                                                                                                                                                    caACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAg                                                                                                                                                                                                                                                      >  7:87098/1‑147 (MQ=255)
                                                                                                                                                    caACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAg                                                                                                                                                                                                                                                      <  8:87098/147‑1 (MQ=255)
                                                                                                                                                     aaCGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCa                                                                                                                                                                                                                                                   >  7:1168199/1‑149 (MQ=255)
                                                                                                                                                     aaCGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCa                                                                                                                                                                                                                                                   >  7:649745/1‑149 (MQ=255)
                                                                                                                                                      aCGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCAt                                                                                                                                                                                                                                                  >  5:515987/1‑149 (MQ=255)
                                                                                                                                                         ttttAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTACTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTa                                                                                                                                                                                                                                               <  7:1233517/149‑1 (MQ=255)
                                                                                                                                                          tttAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAg                                                                                                                                                                                                                                              <  5:896595/149‑1 (MQ=255)
                                                                                                                                                          tttAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAg                                                                                                                                                                                                                                              <  5:724963/149‑1 (MQ=255)
                                                                                                                                                           ttAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAgg                                                                                                                                                                                                                                             <  8:649745/149‑1 (MQ=255)
                                                                                                                                                                 tcatAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCAcg                                                                                                                                                                                                                                       >  6:1433708/1‑149 (MQ=255)
                                                                                                                                                                  catAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGc                                                                                                    >  1:76682/1‑279 (MQ=255)
                                                                                                                                                                   atAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCGCGGTGCTGTATGGATCGGATACTGATATATTTTCGTGCATAGATTTCACCGGCCATCTGAGATCGCCGACAGGATAACATAATTACAGCCGTTTAGTGCCAcgcg                                                                                                                                                                                                                                     >  7:1042278/1‑149 (MQ=255)
                                                                                                                                                                   atAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACa                                                                                                                                                                                                                                                       <  7:345050/131‑1 (MQ=255)
                                                                                                                                                                    tAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCAttt                                                                                                                                                                                                                                                >  7:1129187/1‑137 (MQ=255)
                                                                                                                                                                    tAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCAttt                                                                                                                                                                                                                                                <  8:1129187/137‑1 (MQ=255)
                                                                                                                                                                           gcagcaGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTcagcca                                                                                                                                                                                                                             >  5:334027/1‑149 (MQ=255)
                                                                                                                                                                             atcagTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTcagccagc                                                                                                                                                                                                                           <  7:1376233/147‑1 (MQ=255)
                                                                                                                                                                                agTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGa                                                                                                                                                                                                                        <  8:413914/149‑1 (MQ=255)
                                                                                                                                                                                agTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGa                                                                                                                                                                                                                        <  8:1439909/149‑1 (MQ=255)
                                                                                                                                                                                agTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGa                                                                                                                                                                                                                        <  5:155786/149‑1 (MQ=255)
                                                                                                                                                                                 gTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGAt                                                                                                                                                                                                                       <  7:1466906/149‑1 (MQ=255)
                                                                                                                                                                                   gttcccggCATGGCCGGTCGGGAATTCCCGGTGCGTTTTGAATCCGTTACTG‑‑GTATTTTCCTGCAAGGATTTCACCGGCCATCGTGGATGGCCTGCAGGTTAACATAATTACGGCCTTTTAGGGCCGCGCGCTCAGCCAGCGGCTGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTTTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAAGGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCa                                                                                     <  4:634500/268‑1 (MQ=255)
                                                                                                                                                                                       ccAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAAc                                                                                                                                                                                                                 <  8:249417/149‑1 (MQ=255)
                                                                                                                                                                                        cAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACg                                                                                                                                                                                                                >  6:1514639/1‑149 (MQ=255)
                                                                                                                                                                                          gCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCAt                                                                                                                                                                                                                                                  >  8:545968/1‑113 (MQ=255)
                                                                                                                                                                                          gCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCAt                                                                                                                                                                                                                                                  <  7:545968/113‑1 (MQ=255)
                                                                                                                                                                                          gCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCa                                                                                                                                                                                                              <  6:946049/149‑1 (MQ=255)
                                                                                                                                                                                               gCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTAt                                                                                                                                                                                                         <  6:334027/149‑1 (MQ=255)
                                                                                                                                                                                                ccTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATc                                                                                                                                                                                                        >  7:716611/1‑149 (MQ=255)
                                                                                                                                                                                                ccTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCACCTATc                                                                                                                                                                                                        <  8:921837/149‑1 (MQ=255)
                                                                                                                                                                                                 cTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCa                                                                                                                                                                                                       >  7:1675692/1‑149 (MQ=255)
                                                                                                                                                                                                 cTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCa                                                                                                                                                                                                       >  7:1193916/1‑149 (MQ=255)
                                                                                                                                                                                                  tgtgCGGCATTCCCGGTGCTGTATCAATCGGTTACTGATCTATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCAc                                                                                                                                                                                                      <  7:1110799/149‑1 (MQ=255)
                                                                                                                                                                                                    tgCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCCCCGGCCATCTGAGATGGCCGACAGGTTTACATAATTACAGCCATTTAGGGCCACGCTGTCAGCCAGCGGATTAAAACGCATCTATCACCTTGTGATCGAAACCTTCATTCTCCCCCGCCAGTTCGGAACTGAGCTTCCCCATAAACTCCACTAAAACAAGCCCATTGTGCTTCGCCAGCTTCTTGCCTCTTGCTGTTTGCTAGGTTTGCGGTAg                                                                        >  4:630450/1‑273 (MQ=255)
                                                                                                                                                                                                     gCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCtt                                                                                                                                                                                                   <  8:716611/149‑1 (MQ=255)
                                                                                                                                                                                                       ggCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTtgt                                                                                                                                                                                                 <  5:1327232/149‑1 (MQ=255)
                                                                                                                                                                                                          aTTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCAt                                                                                                                                                                                                             <  8:1055600/134‑1 (MQ=255)
                                                                                                                                                                                                          aTTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGAt                                                                                                                                                                                              >  6:438384/1‑149 (MQ=255)
                                                                                                                                                                                                                     tGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGc                                                                                                                                                                                                                                           >  5:1393860/1‑93 (MQ=255)
                                                                                                                                                                                                                      gTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCAt                                                                                                                                                                                  <  8:1210484/149‑1 (MQ=255)
                                                                                                                                                                                                                           aaTCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCAcc                                                                                                                                                                                                     <  2:563859/125‑1 (MQ=255)
                                                                                                                                                                                                                            aTCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGt                                                                                          >  2:348510/1‑231 (MQ=255)
                                                                                                                                                                                                                                ggtACTGATATATGTTCCGACATTGATTGCACCGGCCAGCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCGCCTTGTGATCGACACCTTCATCCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGaaa                                                  <  2:636874/265‑1 (MQ=255)
                                                                                                                                                                                                                                gTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAg                                                                                                                                                                                                                                              >  7:515085/1‑79 (MQ=255)
                                                                                                                                                                                                                                gTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAg                                                                                                                                                                                                                                              <  8:515085/79‑1 (MQ=255)
                                                                                                                                                                                                                                     tGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTAt                                                                                                                                                                                                         <  2:532942/111‑1 (MQ=255)
                                                                                                                                                                                                                                     tGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGtt                                                                                                                                                                   >  6:809225/1‑149 (MQ=255)
                                                                                                                                                                                                                                     tGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGtt                                                                                                                                                                   >  6:1339812/1‑149 (MQ=255)
                                                                                                                                                                                                                                       atatatTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGTCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCg                                                                                                                                                                 >  8:1105640/1‑149 (MQ=255)
                                                                                                                                                                                                                                         atatTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCg                                                                                                                                                                 <  6:1027208/147‑1 (MQ=255)
                                                                                                                                                                                                                                         atatTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGCGCTTCGCCATAAACTCGACTAAAACATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGCTTTGCTTGGTTTGCGTTAGCTTCGGCAACTTAGTCTGGAAATGATCCAGCGCGTacc                                  >  4:79921/1‑272 (MQ=255)
                                                                                                                                                                                                                                          tatTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAggg                                                                                                                                                                                                                                            >  7:767555/1‑71 (MQ=255)
                                                                                                                                                                                                                                            ttttCCTGCATGGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGtt                               <  2:79921/274‑1 (MQ=255)
                                                                                                                                                                                                                                             tttCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTg                                                                                                                                                           <  6:14598/149‑1 (MQ=255)
                                                                                                                                                                                                                                              ttCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGa                                                                                                                                                          <  6:515987/149‑1 (MQ=255)
                                                                                                                                                                                                                                               tCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGt                                                                                                                         >  1:5563/1‑181 (MQ=255)
                                                                                                                                                                                                                                                  tGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGtt                                                                                                                                                                                                                                                                   <  3:7464/40‑1 (MQ=255)
                                                                                                                                                                                                                                                    cATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTCGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCCTTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTTCGGTAGCTTCAGCCACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCaa                      >  2:661271/1‑275 (MQ=255)
                                                                                                                                                                                                                                                     aTAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATc                                                                                                                                                                                             <  7:234998/107‑1 (MQ=255)
                                                                                                                                                                                                                                                     aTAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGc                                                                                                                                                   <  5:733206/149‑1 (MQ=255)
                                                                                                                                                                                                                                                      tAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGcc                                                                                                                                                  >  5:801045/1‑149 (MQ=255)
                                                                                                                                                                                                                                                        gATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCAt                                                                                                                                                <  5:834348/149‑1 (MQ=255)
                                                                                                                                                                                                                                                         aTTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATCAACTCGACTAAAAAATGCGCGTTGATCTGCGCCAGCTGCTTTCCTCTTTCCGTTTGCATGGTTTGCTGTATCTTACGCTACTTAGTCTGGAAACGATCCAGCCCTTAGCCTTTaattttcacgcg                  >  4:453970/1‑263 (MQ=255)
                                                                                                                                                                                                                                                          tttCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTca                                                                                                                                                                                                                                 <  6:1591158/66‑1 (MQ=255)
                                                                                                                                                                                                                                                             cACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACt                                                                                                                                           >  5:941267/1‑149 (MQ=255)
                                                                                                                                                                                                                                                             cACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCtt                                                                                                      <  3:194263/186‑1 (MQ=255)
                                                                                                                                                                                                                                                              aCCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTc                                                                                                                                          <  8:302906/149‑1 (MQ=255)
                                                                                                                                                                                                                                                               ccGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCAt                                                                                                                                                                                  >  5:360315/1‑108 (MQ=255)
                                                                                                                                                                                                                                                                cGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGa                                                                                                                                        <  7:979829/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                 ggCCATCTCAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGAc                                                                                                                                       >  7:618352/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                   ccATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACt                                                                                                                                                            >  1:182138/1‑126 (MQ=255)
                                                                                                                                                                                                                                                                   ccATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCg                                                                                            <  1:405/190‑1 (MQ=255)
                                                                                                                                                                                                                                                                     aTCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTaaa                                                                                                                                   >  5:961764/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                      tCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTaaaa                                                                                                                                  >  5:1208093/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                 
GATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATGGGCGAACTCTTTCTCACTTAAACGGTGCCATGCCGAGTCATTCGCCAGCGCACGTTTTAAAATTGCCGCGCTCCAGTGATTCTCACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAAC  >  NC_000913/2032071‑2032609

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