Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 604,753 A→C 7.4% intergenic (‑90/+18) ybdG ← / ← nfsB mechanosensitive channel protein, miniconductance/dihydropteridine reductase, NAD(P)H‑dependent, oxygen‑insensitive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913604,7530AC 7.4% 317.5 / 14.8 95intergenic (‑90/+18)ybdG/nfsBmechanosensitive channel protein, miniconductance/dihydropteridine reductase, NAD(P)H‑dependent, oxygen‑insensitive
Reads supporting (aligned to +/- strand):  ref base A (44/44);  new base C (3/4);  total (47/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

GTGGTGTGATCGATCTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCTTTCAGACCAAATT  >  NC_000913/604605‑604901
                                                                                                                                                    |                                                                                                                                                    
gtggtgTGATCGATCTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGCAATAGCCGGGCAGa                                                                                                                                                      <  1:799187/149‑1 (MQ=255)
 tggtgTGATCGATCTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGAt                                                                                                                                                     <  1:121420/149‑1 (MQ=255)
        atcgatcTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCATCTGCCCGGc                                                                                                                                              <  3:670144/149‑1 (MQ=255)
         tcgatcTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCATCTGCCCGGct                                                                                                                                             <  1:273758/149‑2 (MQ=255)
            atcTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAaga                                                                                                                                          <  2:667417/149‑1 (MQ=255)
            atcTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAaga                                                                                                                                          >  2:113755/1‑149 (MQ=255)
              ctcAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAagaga                                                                                                                                        <  2:263017/149‑1 (MQ=255)
               tcAATACCCGCTAAATCTTCAACCTGGGATATTATATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAagagag                                                                                                                                       >  1:227187/1‑149 (MQ=255)
                 aaTACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGaa                                                                                                                                     <  4:769042/149‑1 (MQ=255)
                 aaTACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGaa                                                                                                                                     >  1:728401/1‑149 (MQ=255)
                      ccGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCATCTGCCCGGctatt                                                                                                                                          >  2:491883/1‑135 (MQ=255)
                      ccGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCATCTGCCCGGctatt                                                                                                                                          <  1:491883/139‑5 (MQ=255)
                      ccGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAaga                                                                                                                                          <  1:174243/139‑1 (MQ=255)
                      ccGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAaga                                                                                                                                          >  2:174243/1‑139 (MQ=255)
                      ccGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTaca                                                                                                                                <  3:108314/149‑1 (MQ=255)
                      ccGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTaca                                                                                                                                >  2:308186/1‑149 (MQ=255)
                          tAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTc                                                                                                                            >  2:976917/1‑149 (MQ=255)
                           aaaTCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCg                                                                                                                           <  2:943560/149‑1 (MQ=255)
                              tCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGtt                                                                                                                        <  3:846465/149‑1 (MQ=255)
                                ttCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTaa                                                                                                                      >  4:292236/1‑149 (MQ=255)
                                 tCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAg                                                                                                                     <  2:482734/149‑1 (MQ=255)
                                 tCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAg                                                                                                                     <  3:983539/149‑1 (MQ=255)
                                   aaCCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGt                                                                                                                   <  1:748550/149‑1 (MQ=255)
                                   aaCCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGt                                                                                                                   <  2:514427/149‑1 (MQ=255)
                                   aaCCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGt                                                                                                                   <  3:123471/149‑1 (MQ=255)
                                     ccTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGa                                                                                                                 <  2:283172/149‑1 (MQ=255)
                                      cTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGAt                                                                                                                >  1:554414/1‑149 (MQ=255)
                                       tGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATg                                                                                                               >  4:544768/1‑149 (MQ=255)
                                       tGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATg                                                                                                               >  2:624865/1‑149 (MQ=255)
                                         ggATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGtt                                                                                                             >  4:76209/1‑149 (MQ=255)
                                          gATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGttt                                                                                                            >  2:412229/1‑149 (MQ=255)
                                           atatTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGtttt                                                                                                           >  2:727756/1‑149 (MQ=255)
                                             atTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGc                                                                                                         <  1:981027/149‑1 (MQ=255)
                                              ttAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCg                                                                                                        >  3:505340/1‑149 (MQ=255)
                                                 aaTCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCa                                                                                                     <  2:728401/149‑1 (MQ=255)
                                                   tCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGa                                                                                                   >  4:388608/1‑149 (MQ=255)
                                                         aTAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGAt                                                                                             <  1:412229/149‑1 (MQ=255)
                                                           aaCGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGAttt                                                                                           >  1:249951/1‑149 (MQ=255)
                                                                ctcCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCa                                                                                      <  2:830485/149‑1 (MQ=255)
                                                                 tcCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCag                                                                                     >  3:765023/1‑149 (MQ=255)
                                                                 tcCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCag                                                                                     >  4:430783/1‑149 (MQ=255)
                                                                   cTTTATACAACAGCAGCCTATGATGACGACTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACGTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCagcg                                                                                   >  1:400072/1‑149 (MQ=255)
                                                                    tttAAACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGATGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCagcgt                                                                                  >  4:753752/1‑149 (MQ=255)
                                                                         acaacaGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTaa                                                                                                                      >  2:729808/1‑108 (MQ=255)
                                                                         acaacaGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTaa                                                                                                                      <  1:729808/108‑1 (MQ=255)
                                                                         acaacaGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCagcgtagcgt                                                                             >  1:607887/1‑149 (MQ=255)
                                                                           aacaGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGtt                                                                                                             >  1:971469/1‑115 (MQ=255)
                                                                           aacaGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGtt                                                                                                             <  2:971469/115‑1 (MQ=255)
                                                                           aacaGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTa                                                                           >  1:930085/1‑149 (MQ=255)
                                                                           aacaGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTa                                                                           >  4:564061/1‑149 (MQ=255)
                                                                            acaGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTaa                                                                          <  3:321557/149‑1 (MQ=255)
                                                                              agcagcCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTaaaa                                                                        <  3:544768/149‑1 (MQ=255)
                                                                                 agcCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCt                                                                     <  3:292236/149‑1 (MQ=255)
                                                                                 agcCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCt                                                                     >  4:988841/1‑149 (MQ=255)
                                                                                  gcCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACgaga                                                                                              <  3:331745/123‑1 (MQ=255)
                                                                                  gcCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACgaga                                                                                              >  4:331745/1‑123 (MQ=255)
                                                                                    cTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGCGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGCATTACACTTCGGTTAAAGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCa                                                                  >  1:247485/1‑149 (MQ=255)
                                                                                     tATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCaa                                                                 <  4:217450/149‑1 (MQ=255)
                                                                                      atgatgACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAAc                                                                >  4:554652/1‑149 (MQ=255)
                                                                                           aaCGGCTGAAACAGGGTTATGCAAATCAAGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCtgt                                                           <  2:339922/148‑1 (MQ=255)
                                                                                             cGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGtga                                                         <  3:988841/149‑1 (MQ=255)
                                                                                              ggCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGtgat                                                        <  4:1013664/149‑1 (MQ=255)
                                                                                               gCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGtgatg                                                       <  3:684918/149‑1 (MQ=255)
                                                                                                     aCAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTAc                                                 <  1:985036/149‑1 (MQ=255)
                                                                                                      cAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTAcc                                                >  2:319032/1‑149 (MQ=255)
                                                                                                         ggTTATGCAAATCAGGAGAATGTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGa                                             >  4:715552/1‑149 (MQ=255)
                                                                                                           ttATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGaac                                           >  2:279778/1‑149 (MQ=255)
                                                                                                              tGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGaacaac                                        <  1:754583/149‑1 (MQ=255)
                                                                                                                 aaaTCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAaccac                                     <  3:388608/149‑1 (MQ=255)
                                                                                                                   aTCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAAccacca                                   <  4:420742/149‑1 (MQ=255)
                                                                                                                   aTCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGATAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAAccacca                                   >  3:590710/1‑149 (MQ=255)
                                                                                                                    tCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAg                                  <  4:389305/149‑1 (MQ=255)
                                                                                                                    tCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAg                                  >  2:574489/1‑149 (MQ=255)
                                                                                                                      aGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGAc                                >  1:635968/1‑149 (MQ=255)
                                                                                                                      aGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGAc                                >  4:753272/1‑149 (MQ=255)
                                                                                                                        gagaATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTg                              >  1:588437/1‑149 (MQ=255)
                                                                                                                        gagaATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCCGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAaccac                                     <  1:388056/142‑1 (MQ=255)
                                                                                                                        gagaATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCCGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAaccac                                     >  2:388056/1‑142 (MQ=255)
                                                                                                                          gaATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGgt                            <  4:662825/149‑1 (MQ=255)
                                                                                                                           aaTCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGgtg                           <  3:715552/149‑1 (MQ=255)
                                                                                                                            aTCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGgtgt                          >  2:366505/1‑149 (MQ=255)
                                                                                                                             tCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTa                         <  3:45779/149‑1 (MQ=255)
                                                                                                                               tGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGtgat                                                        <  3:969680/116‑1 (MQ=255)
                                                                                                                               tGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGtgat                                                        >  4:969680/1‑116 (MQ=255)
                                                                                                                               tGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGc                       <  3:430783/149‑1 (MQ=255)
                                                                                                                                   agGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCttt                   >  1:670852/1‑149 (MQ=255)
                                                                                                                                     gAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTtct                 >  1:144610/1‑149 (MQ=255)
                                                                                                                                       tcttgCCGGGCATCTGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATc                                                                      >  4:1004681/5‑94 (MQ=255)
                                                                                                                                       tcttgCCGGGCATCTGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATc                                                                      <  3:1004681/90‑1 (MQ=255)
                                                                                                                                       aaTAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCtgt                                                           >  4:224207/1‑105 (MQ=255)
                                                                                                                                       aaTAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCtgt                                                           <  3:224207/105‑1 (MQ=255)
                                                                                                                                        aTAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCTGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCtt              >  3:1030971/1‑149 (MQ=255)
                                                                                                                                        aTAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCtt              <  3:846828/149‑1 (MQ=255)
                                                                                                                                         ttgCCGGGCATCTGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCttt             >  2:819673/3‑149 (MQ=255)
                                                                                                                                           gCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCTTTCa           <  4:488607/149‑1 (MQ=255)
                                                                                                                                             cGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCTTTCAGa         <  2:141813/149‑1 (MQ=255)
                                                                                                                                                    aTGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCTTTCAGACCAAAtt  <  2:1032106/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTGGTGTGATCGATCTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCTTTCAGACCAAATT  >  NC_000913/604605‑604901

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