Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 314,255 T→C 14.1% intergenic (‑450/‑102) ykgP ← / → eaeH pseudogene, oxidoreductase family/pseudogene, attaching and effacing protein homology;factor; Not classified

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913314,2550TC14.1% 249.8 / 24.8 93intergenic (‑450/‑102)ykgP/eaeHpseudogene, oxidoreductase family/pseudogene, attaching and effacing protein homology;factor; Not classified
Reads supporting (aligned to +/- strand):  ref base T (40/39);  new base C (7/6);  total (48/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCTATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATTC  >  NC_000913/314111‑314403
                                                                                                                                                |                                                                                                                                                    
aTCTATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCCGcag                                                                                                                                                  >  4:6546/1‑147 (MQ=255)
    aTTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGaaaa                                                                                                                                              <  3:257853/149‑1 (MQ=255)
    aTTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGaaa                                                                                                                                               <  2:805871/148‑1 (MQ=255)
    aTTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGaaa                                                                                                                                               <  2:16155/148‑1 (MQ=255)
    aTTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGaaa                                                                                                                                               <  1:471712/148‑1 (MQ=255)
    aTTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCCGCAGaaa                                                                                                                                               <  1:448143/148‑1 (MQ=255)
     tttAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGaaaa                                                                                                                                              <  1:343886/148‑1 (MQ=255)
       tAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTTCATAAAGGTATTTCCTTTTCCGCAGAAAAGGa                                                                                                                                           <  3:592366/149‑1 (MQ=255)
       tAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCCGCAGAAAAGGa                                                                                                                                           <  3:6546/149‑1 (MQ=255)
        aaTGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCCGCAGAAAAGGaa                                                                                                                                          <  3:731095/149‑1 (MQ=255)
              ttGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCgg                                                                                                                                     >  3:225329/1‑148 (MQ=255)
              ttGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCgg                                                                                                                                     <  4:225329/148‑1 (MQ=255)
               tGCACAAAAAGCTAAACTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTTCGGAAAAGGAAATCggg                                                                                                                                    <  2:691695/148‑1 (MQ=255)
                   caAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACCTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATc                                                                                                                               >  4:286195/1‑149 (MQ=255)
                      aaaGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATcccc                                                                                                                            >  3:708690/1‑149 (MQ=255)
                       aaGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATcccc                                                                                                                            >  1:876359/1‑148 (MQ=255)
                       aaGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCg                                                                                                                           >  4:805325/1‑149 (MQ=255)
                       aaGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCg                                                                                                                           <  4:341593/149‑1 (MQ=255)
                           tAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACTTAAAGGTATTTCCTTTTCCGCAGAAAAGGaa                                                                                                                                          <  1:447244/130‑1 (MQ=255)
                           tAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCCGCAGAAAAGGaa                                                                                                                                          >  2:447244/1‑130 (MQ=255)
                           tAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCCGCAGAAAAGGAAAtaccttt                                                                                                                                  >  3:420357/1‑132 (MQ=255)
                           tAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCCGCAGAAAAGGAAAtaccttt                                                                                                                                  <  4:420357/138‑7 (MQ=255)
                            aaaTTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGttt                                                                                                                       >  1:303926/1‑148 (MQ=255)
                              aTTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGttttt                                                                                                                     <  2:303926/148‑1 (MQ=255)
                                    ccccATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACaa                                                                                                              >  4:121575/1‑149 (MQ=255)
                                      ccATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGc                                                                                                            <  3:678681/149‑1 (MQ=255)
                                      ccATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTAACAAGc                                                                                                            >  4:644206/1‑149 (MQ=255)
                                         taaataaaAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGa                                                                                                         >  3:869052/1‑149 (MQ=255)
                                            ataaaAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACaa                                                                                                              <  3:346349/141‑1 (MQ=255)
                                            ataaaAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACaa                                                                                                              >  4:346349/1‑141 (MQ=255)
                                                aaaTATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGAc                                                                                                                >  4:260163/1‑135 (MQ=255)
                                                aaaTATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGAc                                                                                                                <  3:260163/135‑1 (MQ=255)
                                                aaaTATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCatt                                                                                                  >  4:189348/1‑149 (MQ=255)
                                                aaaTATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCatt                                                                                                  >  3:130897/1‑149 (MQ=255)
                                                  atatAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCattat                                                                                                <  3:121575/149‑1 (MQ=255)
                                                  atatAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCattat                                                                                                >  3:890577/1‑149 (MQ=255)
                                                   tataATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCattatt                                                                                               >  3:190815/1‑149 (MQ=255)
                                                      aaTCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTtgt                                                                                            <  4:890577/149‑1 (MQ=255)
                                                       aTCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCgg                                                                                                                                     <  3:901242/107‑1 (MQ=255)
                                                       aTCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCgg                                                                                                                                     >  4:901242/1‑107 (MQ=255)
                                                            gcgcCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGc                                                                                      <  4:869052/149‑1 (MQ=255)
                                                                      accTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGt                                                                            >  3:460582/1‑149 (MQ=255)
                                                                       ccTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATTAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGtt                                                                           <  3:644206/149‑1 (MQ=255)
                                                                          gatgaGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTATTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTaa                                                                         >  1:267394/1‑148 (MQ=255)
                                                                            tgaGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACCTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTc                                                                      <  3:286195/149‑1 (MQ=255)
                                                                               gTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGtt                                                                   <  4:460582/149‑1 (MQ=255)
                                                                                tGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTCCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTtc                                                                  <  3:288656/149‑1 (MQ=255)
                                                                                   cTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCatta                                                                                                 >  3:690754/1‑115 (MQ=255)
                                                                                   cTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCatta                                                                                                 <  4:690754/115‑1 (MQ=255)
                                                                                    tataGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCAGCGGAAAAGGAAATCGGGTAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTtcatca                                                              >  4:419129/1‑149 (MQ=255)
                                                                                     ataGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAg                                                             >  3:588228/1‑149 (MQ=255)
                                                                                          cACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAAt                                                        >  4:830000/1‑149 (MQ=255)
                                                                                            cTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAAtt                                                       <  2:185666/148‑1 (MQ=255)
                                                                                              ggATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTCCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTAtc                                                    >  3:156296/1‑149 (MQ=255)
                                                                                               gATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTAtct                                                   >  3:169192/1‑149 (MQ=255)
                                                                                               gATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTAtct                                                   >  3:724977/1‑149 (MQ=255)
                                                                                                 tatatTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCa                                                 >  3:86031/1‑149 (MQ=255)
                                                                                                  atatTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCaa                                                >  4:257915/1‑149 (MQ=255)
                                                                                                     ttAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGt                                             >  3:629585/1‑149 (MQ=255)
                                                                                                      tAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGcc                                                                                     <  1:806885/108‑1 (MQ=255)
                                                                                                      tAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGcc                                                                                     >  2:806885/1‑108 (MQ=255)
                                                                                                      tAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGt                                             >  2:286320/1‑148 (MQ=255)
                                                                                                        ggTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCgg                                                                                                                          <  3:635731/69‑1 (MQ=255)
                                                                                                        ggTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCgg                                                                                                                          >  4:635731/1‑69 (MQ=255)
                                                                                                        ggTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTAGGTTCTTTAATTCGTTCATCAGGAAATTATCTCAATGTcac                                          >  4:587730/1‑149 (MQ=255)
                                                                                                          tGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTcaca                                         >  1:510257/1‑148 (MQ=255)
                                                                                                           ggCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACAtt                                       <  3:189348/149‑1 (MQ=255)
                                                                                                              ggTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATtata                                    <  3:875968/149‑1 (MQ=255)
                                                                                                               gTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATaa                                   >  4:5272/1‑149 (MQ=255)
                                                                                                                 gCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATaaaa                                 <  4:130897/149‑1 (MQ=255)
                                                                                                                  cACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAAc                                >  3:849467/1‑149 (MQ=255)
                                                                                                                    cTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAg                              <  4:724977/149‑1 (MQ=255)
                                                                                                                     tttCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAgg                             <  3:257915/149‑1 (MQ=255)
                                                                                                                      ttCTTACATAAAGGTATTCCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGt                            <  4:156296/149‑1 (MQ=255)
                                                                                                                           aCATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAAc                                <  2:831680/140‑1 (MQ=255)
                                                                                                                           aCATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAAc                                >  1:831680/1‑140 (MQ=255)
                                                                                                                            cATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATaaa                      <  4:849467/149‑1 (MQ=255)
                                                                                                                               aaaGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAaacaa                   <  3:419129/149‑1 (MQ=255)
                                                                                                                                aaGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTCTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAaacaa                   >  2:275368/1‑148 (MQ=255)
                                                                                                                                  ccgaTTTCCTTTTCCGCAGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCa                >  3:134457/4‑149 (MQ=255)
                                                                                                                                     aTTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTCTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACg              <  1:275368/148‑1 (MQ=255)
                                                                                                                                        tCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGAttt          >  3:8310/1‑149 (MQ=255)
                                                                                                                                        tCCTTTTCCGCAGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGAttt          >  4:258485/1‑149 (MQ=255)
                                                                                                                                        tCCTTTTCCGCAGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGAtt           >  2:321871/1‑148 (MQ=255)
                                                                                                                                        tCCTTTTCCGCAGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGAtt           >  2:883303/1‑148 (MQ=255)
                                                                                                                                          cTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGt                                                                    >  4:392810/1‑89 (MQ=255)
                                                                                                                                          cTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGt                                                                    <  3:392810/89‑1 (MQ=255)
                                                                                                                                          cTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAAc                                <  4:708995/125‑1 (MQ=255)
                                                                                                                                          cTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAAc                                >  3:708995/1‑125 (MQ=255)
                                                                                                                                           ttttCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAg                              <  3:885949/126‑1 (MQ=255)
                                                                                                                                           ttttCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAg                              >  4:885949/1‑126 (MQ=255)
                                                                                                                                           ttttCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCg        <  2:267394/148‑1 (MQ=255)
                                                                                                                                            tttCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGtt      <  3:805325/149‑1 (MQ=255)
                                                                                                                                                tGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATTc  >  4:730775/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
ATCTATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATTC  >  NC_000913/314111‑314403

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: