Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 604,753 A→C 7.1% 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.1% 330.1 / 16.5 99intergenic (‑90/+18)ybdG/nfsBmechanosensitive channel protein, miniconductance/dihydropteridine reductase, NAD(P)H‑dependent, oxygen‑insensitive
Reads supporting (aligned to +/- strand):  ref base A (52/40);  new base C (3/4);  total (55/44)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.57e-01

TGGTGTGATCGATCTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCTTTCAGACCAAA  >  NC_000913/604606‑604899
                                                                                                                                                   |                                                                                                                                                  
tggtgTGATCGATCTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGAt                                                                                                                                                   <  4:596654/149‑1 (MQ=255)
   tgtgATCGATCTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGcc                                                                                                                                                <  2:205843/149‑1 (MQ=255)
     tgATCGATCTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCg                                                                                                                                              <  2:890647/149‑1 (MQ=255)
      gatcgatcTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCgg                                                                                                                                             >  4:336179/1‑149 (MQ=255)
       atcgatcTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGc                                                                                                                                            >  1:521035/1‑149 (MQ=255)
       atcgatcTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGc                                                                                                                                            <  3:611724/149‑1 (MQ=255)
         cgatcTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGc                                                                                                                                                 >  3:250911/1‑142 (MQ=255)
         cgatcTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGc                                                                                                                                                 <  4:250911/142‑1 (MQ=255)
          gatcTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGTAGAATCTGAGAGGAAATAGCCGGGCATCTGCCCGGctat                                                                                                                                         <  4:841160/149‑4 (MQ=255)
          gatcTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCATCTGCCCGGctat                                                                                                                                         <  4:748360/149‑4 (MQ=255)
          gatcTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAag                                                                                                                                         >  3:471251/1‑149 (MQ=255)
          gatcTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAag                                                                                                                                         >  4:547391/1‑149 (MQ=255)
           atcTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAaga                                                                                                                                        >  1:193664/1‑149 (MQ=255)
             ctcAATACCCGCTAAATCTTCCACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAagaga                                                                                                                                      >  4:374614/1‑149 (MQ=255)
             ctcAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAagaga                                                                                                                                      <  4:418080/149‑1 (MQ=255)
              tcAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAagagag                                                                                                                                     >  2:14608/1‑149 (MQ=255)
                aaTACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGaa                                                                                                                                   >  2:58736/1‑149 (MQ=255)
                    cccGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTac                                                                                                                               <  1:447894/149‑1 (MQ=255)
                    cccGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTac                                                                                                                               <  1:785055/149‑1 (MQ=255)
                     ccGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTaca                                                                                                                              >  2:33306/1‑149 (MQ=255)
                     ccGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTaca                                                                                                                              <  3:336179/149‑1 (MQ=255)
                      cGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTacac                                                                                                                             >  1:598719/1‑149 (MQ=255)
                       gCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACt                                                                                                                            >  1:591918/1‑149 (MQ=255)
                        cTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACtt                                                                                                                           >  1:817534/1‑149 (MQ=255)
                           aaTCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCgg                                                                                                                        <  3:562432/149‑1 (MQ=255)
                             tCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGtt                                                                                                                      >  4:814618/1‑149 (MQ=255)
                              cTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATTACGGCTGTAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTa                                                                                                                     >  2:219204/1‑149 (MQ=255)
                              cTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTa                                                                                                                     >  3:125112/1‑149 (MQ=255)
                               ttCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCTGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTaa                                                                                                                    >  1:103986/1‑149 (MQ=255)
                               ttCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTaa                                                                                                                    >  1:566610/1‑149 (MQ=255)
                                tCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAg                                                                                                                   <  1:454715/149‑1 (MQ=255)
                                 cAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAgg                                                                                                                  >  3:921221/1‑149 (MQ=255)
                                            atTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGc                                                                                                       <  4:495107/149‑1 (MQ=255)
                                             ttAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCATCTGCCCGGc                                                                                                                                            >  3:700237/1‑111 (MQ=255)
                                             ttAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCATCTGCCCGGc                                                                                                                                            <  4:700237/111‑1 (MQ=255)
                                             ttAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCg                                                                                                      >  2:806508/1‑149 (MQ=255)
                                             ttAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCg                                                                                                      <  2:677120/149‑1 (MQ=255)
                                              tAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCgg                                                                                                     <  4:921221/149‑1 (MQ=255)
                                                aaTCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCa                                                                                                   <  2:193664/149‑1 (MQ=255)
                                                  tCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGa                                                                                                 >  3:531089/1‑149 (MQ=255)
                                                   ccTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGAc                                                                                                <  1:33306/149‑1 (MQ=255)
                                                     tGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACga                                                                                              >  2:743101/1‑149 (MQ=255)
                                                           aCGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGAc                                                                                                <  1:27703/141‑1 (MQ=255)
                                                           aCGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGAc                                                                                                >  2:27703/1‑141 (MQ=255)
                                                           aCGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTc                                                                                        <  2:472189/149‑1 (MQ=255)
                                                           aCGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTc                                                                                        <  2:817534/149‑1 (MQ=255)
                                                             gTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCgg                                                                                      >  1:35554/1‑149 (MQ=255)
                                                                  cTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCagcg                                                                                 >  1:207599/1‑149 (MQ=255)
                                                                  cTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCagcg                                                                                 >  2:87810/1‑149 (MQ=255)
                                                                  cTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCagcg                                                                                 >  4:839753/1‑149 (MQ=255)
                                                                   tttATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTaca                                                                                                                              <  2:505534/103‑1 (MQ=255)
                                                                   tttATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTaca                                                                                                                              >  1:505534/1‑103 (MQ=255)
                                                                   tttATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCagcgt                                                                                <  3:114684/149‑1 (MQ=255)
                                                                   tttATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCagcgt                                                                                >  2:362683/1‑149 (MQ=255)
                                                                          aacaGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTa                                                                         >  4:336402/1‑149 (MQ=255)
                                                                          aacaGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTa                                                                         <  3:336402/149‑1 (MQ=255)
                                                                          aacaGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTa                                                                         >  3:580289/1‑149 (MQ=255)
                                                                          aacaGCAGCCAATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTa                                                                         >  4:124563/1‑149 (MQ=255)
                                                                           acaGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTaa                                                                        >  1:511152/1‑149 (MQ=255)
                                                                                        atgaCGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCt                                                           >  4:543651/1‑149 (MQ=255)
                                                                                        atgaCGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCt                                                           <  3:814618/149‑1 (MQ=255)
                                                                                         tgaCGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCtg                                                          <  4:580289/149‑1 (MQ=255)
                                                                                             ggCTGAAACAGGGTTATGCGAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGtgat                                                      >  1:334462/1‑149 (MQ=255)
                                                                                             ggCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGtgat                                                      >  2:147871/1‑149 (MQ=255)
                                                                                              gCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGtgatg                                                     >  2:324573/1‑149 (MQ=255)
                                                                                                tGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGAc                                                   <  2:167090/149‑1 (MQ=255)
                                                                                                    aCAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTAc                                               >  1:515634/1‑149 (MQ=255)
                                                                                                                aaaTCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAaccac                                   >  3:790954/1‑149 (MQ=255)
                                                                                                                 aaTCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAaccacc                                  >  4:444803/1‑149 (MQ=255)
                                                                                                                   tCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAg                                >  3:436995/1‑149 (MQ=255)
                                                                                                                   tCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAg                                <  4:306899/149‑1 (MQ=255)
                                                                                                                     aGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGAc                              >  4:207791/1‑149 (MQ=255)
                                                                                                                      ggAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACt                             >  2:279724/1‑149 (MQ=255)
                                                                                                                       gagaATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTg                            >  4:342535/1‑149 (MQ=255)
                                                                                                                       gagaATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTg                            >  1:428124/1‑149 (MQ=255)
                                                                                                                       gagaATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTg                            >  4:419693/1‑149 (MQ=255)
                                                                                                                        agaATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTgg                           <  3:288713/149‑1 (MQ=255)
                                                                                                                        agaATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTgg                           >  2:284655/1‑149 (MQ=255)
                                                                                                                         gaATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGgt                          <  1:284655/149‑1 (MQ=255)
                                                                                                                         gaATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGgt                          <  4:436995/149‑1 (MQ=255)
                                                                                                                             cTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAg                      <  1:70907/149‑1 (MQ=255)
                                                                                                                                agagGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCt                   >  4:925906/1‑149 (MQ=255)
                                                                                                                                 gagGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCtt                  >  2:713094/1‑149 (MQ=255)
                                                                                                                                 gagGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCtt                  >  4:297023/1‑149 (MQ=255)
                                                                                                                                  agGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCttt                 >  1:359123/1‑149 (MQ=255)
                                                                                                                                   ggAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTtc                <  1:147871/149‑1 (MQ=255)
                                                                                                                                   ggAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTtc                <  4:471251/149‑1 (MQ=255)
                                                                                                                                   ggAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTtc                <  3:547391/149‑1 (MQ=255)
                                                                                                                                    tctcttgCCGGGCATCTGCCCGGCTAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAaccac                                   <  2:743061/123‑1 (MQ=255)
                                                                                                                                    tctcttgCCGGGCATCTGCCCGGCTAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAaccac                                   >  1:743061/7‑129 (MQ=255)
                                                                                                                                    gAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTtct               >  3:745609/1‑149 (MQ=255)
                                                                                                                                       aTAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCtt            <  2:521035/149‑1 (MQ=255)
                                                                                                                                        ttgCCGGGCATCTGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCttt           >  4:473219/3‑149 (MQ=255)
                                                                                                                                        tAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCttt           <  3:124563/149‑1 (MQ=255)
                                                                                                                                          gCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCTTTCa         <  2:207599/149‑1 (MQ=255)
                                                                                                                                          gCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCTTTCa         <  3:917586/149‑1 (MQ=255)
                                                                                                                                            cGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCTTTCAGa       <  3:207791/149‑1 (MQ=255)
                                                                                                                                            cGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCTTTCAGa       <  1:713094/149‑1 (MQ=255)
                                                                                                                                                 aGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCTTTCAGACCaaa  <  2:511152/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TGGTGTGATCGATCTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCTTTCAGACCAAA  >  NC_000913/604606‑604899

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