Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 964,922 T→A 5.7% G201G (GGT→GGA ycaI → ComEC family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913964,9220TA 5.7% 299.0 / 5.6 86G201G (GGT→GGAycaIComEC family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (37/44);  new base A (2/3);  total (39/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

ATGCCATTGCCCAGCATCAGCCGCTCACCGGCCGTTTTCTGCAGGCAAGTGTCATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGCTTTTGCGGCGTTGCTGGCTGCCGGACTCATTCGCAGTGGACAAATTTTTCT  >  NC_000913/964774‑965068
                                                                                                                                                    |                                                                                                                                                  
aTGCCATTGCCCAGCATCAGCCGCTCACCGGCCGTTTTCTGCAGGCAAGTGTCATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGt                                                                                                                                                    >  2:603465/1‑149 (MQ=255)
 tGCCATTGCCCAGCATCAGCCGCTCACCGGCCGTTTTCTGCAGGCAAGTGTCATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTa                                                                                                                                                   >  2:11479/1‑149 (MQ=255)
  gCCATTGCCCAGCATCAGCCGCTCACCGGCCGTTTTCTGCAGGCAAGTGTCATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTAt                                                                                                                                                  >  3:794797/1‑149 (MQ=255)
    cATTGCCCAGCATCAGCCGCTCACCGGCCGTTTTCTGCAGGCAAGTGTCATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATgg                                                                                                                                                <  1:767947/149‑1 (MQ=255)
    cATTGCCCAGCATCAGCCGCTCACCGGCCGTTTTCTGCAGGCAAGTGTCATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATgg                                                                                                                                                <  4:23186/149‑1 (MQ=255)
     aTTGCCCAGCATCAGCCGCTCACCGGCCGTTTTCTGCAGGCAAGTGTCATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATggg                                                                                                                                               <  4:215185/149‑1 (MQ=255)
     aTTGCCCAGCATCAGCCGCTCACCGGCCGTTTTCTGCAGGCAAGTGTCATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATggg                                                                                                                                               <  4:100955/149‑1 (MQ=255)
      ttGCCCAGCATCAGCCGCTCACCGGCCGTTTTCTGCAGGCAAGTGTCATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATgggg                                                                                                                                              <  4:740244/149‑1 (MQ=255)
      ttGCCCAGCATCAGCCGCTCACCGGCCGTTTTCTGCAGGCAAGTGTCATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATgggg                                                                                                                                              <  2:232620/149‑1 (MQ=255)
         cccAGCATCAGCCGCTCACCGGCCGTTTTCTGCAGGCAAGTGTCATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGaa                                                                                                                                           <  4:107428/149‑1 (MQ=255)
          ccAGCATCAGCCGCTCACCGGCCGTTTTCTGCAGGCAAGTGTCATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAAc                                                                                                                                          >  1:238247/1‑149 (MQ=255)
            aGCATCAGCCGCTCACCGGCCGTTTTCTGCAGGCAAGTGTCATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACgg                                                                                                                                        <  4:523625/149‑1 (MQ=255)
             gCATCAGCCGCTCACCGGCCGTTTTCTGCAGGCAAGTGTCATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGt                                                                                                                                       <  4:794797/149‑1 (MQ=255)
                tCAGCCGCTCACCGGCCGTTTTCTGCAGGCAAGTGTCATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTAt                                                                                                                                    <  1:289318/149‑1 (MQ=255)
                  aGCCGCTCACCGGCCGTTTTCTGCAGGCAAGTGTCATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATcc                                                                                                                                  <  1:786897/149‑1 (MQ=255)
                     cGCTCACCGGCCGTTTTCAGCAGGCAAGTGTCATTAAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAATAACGCTGCAACCCTATCCGTGGAATGCGGTTATCCTTGGTTTAGGTATGGGGGAACGGTTATCCGTc                                                                                                                               <  4:572029/149‑1 (MQ=255)
                      gCTCACCGGCCGTTTTCTGCAGGCAAGTGTCATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTcc                                                                                                                              >  1:440911/1‑149 (MQ=255)
                           ccGGCCGTTTTCTGCAGGCAAGTGTCATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCaaa                                                                                                                         <  1:381267/149‑1 (MQ=255)
                                cGTTTTCTGCAGGCAAGTGTCATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAAt                                                                                                                    <  4:96622/149‑1 (MQ=255)
                                      cTGCAGGCAAGTGTCATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGGATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCaaaaa                                                                                                              <  1:438744/149‑1 (MQ=255)
                                      cTGCAGGCAAGTGTCATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCaaaaa                                                                                                              <  2:440911/149‑1 (MQ=255)
                                       tGCAGGCAAGTGTCATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGAATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAat                                                                                                             >  3:293252/1‑149 (MQ=255)
                                       tGCAGGCAAGTGTCATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGAATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAat                                                                                                             >  1:425187/1‑149 (MQ=255)
                                          aGGCAAGTGTCATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATc                                                                                                          <  4:293896/149‑1 (MQ=255)
                                           ggCAAGTGTCATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCaaa                                                                                                                <  3:2169/142‑1 (MQ=255)
                                           ggCAAGTGTCATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCaaa                                                                                                                >  4:2169/1‑142 (MQ=255)
                                                 tgtCATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGcgcg                                                                                                   >  4:582244/1‑149 (MQ=255)
                                                   tCATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGAt                                                                                                 >  4:159782/1‑149 (MQ=255)
                                                   tCATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAAGCAAAAATATCATGCGCGAt                                                                                                 >  2:584764/1‑149 (MQ=255)
                                                     aTTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATAc                                                                                               >  2:159228/1‑149 (MQ=255)
                                                      ttGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACt                                                                                              >  2:622407/1‑149 (MQ=255)
                                                              aaTTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATc                                                                                                          >  3:554015/1‑129 (MQ=255)
                                                              aaTTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGCAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATc                                                                                                          <  4:554015/129‑1 (MQ=255)
                                                                ttGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGgcgc                                                                                    <  3:159782/149‑1 (MQ=255)
                                                                      ccTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTaa                                                                              <  1:622407/149‑1 (MQ=255)
                                                                       cTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAAt                                                                             >  3:435782/1‑149 (MQ=255)
                                                                         gCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATgg                                                                           <  1:334859/149‑1 (MQ=255)
                                                                         gCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATgg                                                                           <  1:811062/149‑1 (MQ=255)
                                                                          cGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGc                                                                          <  2:812055/149‑1 (MQ=255)
                                                                              cacaGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGat                                                                       >  4:802564/1‑148 (MQ=255)
                                                                              cacaGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGat                                                                       <  3:802564/148‑1 (MQ=255)
                                                                              cacaGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGat                                                                       >  4:179379/1‑148 (MQ=255)
                                                                              cacaGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGat                                                                       <  3:179379/148‑1 (MQ=255)
                                                                                  gTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGGGCGATACTGGAACGGCGCACTTAACGGCGTTATCgg                                                                  <  4:398460/149‑1 (MQ=255)
                                                                                   tATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCggg                                                                 >  3:240829/1‑149 (MQ=255)
                                                                                    aTCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGAATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGa                                                                <  4:293252/149‑1 (MQ=255)
                                                                                    aTCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGAATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGa                                                                <  2:425187/149‑1 (MQ=255)
                                                                                       tGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTg                                                             >  2:127659/1‑149 (MQ=255)
                                                                                        ggCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGc                                                            <  1:11479/149‑1 (MQ=255)
                                                                                        ggCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGc                                                            <  1:603465/149‑1 (MQ=255)
                                                                                         gCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGca                                                           >  1:517958/1‑149 (MQ=255)
                                                                                          cGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGcac                                                          >  2:766650/1‑149 (MQ=255)
                                                                                              aCTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCg                                                      >  3:93580/1‑149 (MQ=255)
                                                                                                tACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGCt                                                    >  4:85356/1‑149 (MQ=255)
                                                                                                tACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGCt                                                    >  2:717008/1‑149 (MQ=255)
                                                                                                tACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGCt                                                    >  2:434676/1‑149 (MQ=255)
                                                                                                   aatcaacGCTGCAACCCTGTCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCg                                                                   >  4:535580/4‑131 (MQ=255)
                                                                                                   aaacaaCGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCg                                                                   <  3:535580/131‑1 (MQ=255)
                                                                                                       aacGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCggg                                                                 >  1:671492/1‑129 (MQ=255)
                                                                                                       aacGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCggg                                                                 <  2:671492/129‑1 (MQ=255)
                                                                                                           cTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGCTTTTGCGGCGtt                                         <  1:566595/149‑1 (MQ=255)
                                                                                                           cTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGAATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATCGCACATCGCTTTTGCGGCGtt                                         <  1:584764/149‑1 (MQ=255)
                                                                                                            tGCAACCCTATCCGGGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGCTTTTGCGGCGTTg                                        >  3:55442/1‑149 (MQ=255)
                                                                                                                aCCCTATCCGTGGCATGCGGGTATTCTTGGTCTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGCTTTTGCGGCGTTgctgg                                    >  2:90502/1‑149 (MQ=255)
                                                                                                                      tCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCAtt                                                                                 >  1:192314/1‑98 (MQ=255)
                                                                                                                      tCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCAtt                                                                                 <  2:192314/98‑1 (MQ=255)
                                                                                                                        cGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGCTTTTGCGGCGTTGCTGGCTGCCGGa                            <  3:582244/149‑1 (MQ=255)
                                                                                                                               tGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGCTTTTGCGGCGTTGCTGGCTGcc                               <  3:136453/139‑1 (MQ=255)
                                                                                                                               tGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGCTTTTGCGGCGTTGCTGGCTGcc                               >  4:136453/1‑139 (MQ=255)
                                                                                                                               tGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGCTTTTGCGGCGTTGCTGGCTGcc                               <  3:157812/139‑1 (MQ=255)
                                                                                                                               tGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGCTTTTGCGGCGTTGCTGGCTGcc                               >  4:157812/1‑139 (MQ=255)
                                                                                                                                   gTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGCTTTTGCGGCGTTgct                                      >  3:614528/1‑128 (MQ=255)
                                                                                                                                   gTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGCTTTTGCGGCGTTgct                                      <  4:614528/128‑1 (MQ=255)
                                                                                                                                   gTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGCTTTTGCGGCGTTGCTGGCTGCCGGACTCATTCGCAg                 >  1:341830/1‑149 (MQ=255)
                                                                                                                                   gTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGCTTTTGCGGCGTTGCTGGCTGCCGGACTCATTCGCAg                 <  1:717008/149‑1 (MQ=255)
                                                                                                                                   gTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGCTTTTGCGGCGTTGCTGGCTGCCGGACTCATTCGCAg                 <  1:434676/149‑1 (MQ=255)
                                                                                                                                    ttATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGCTTTTGCGGCGTTGCTGGCTGCCGGACTCATTCGCAGt                <  1:159228/149‑1 (MQ=255)
                                                                                                                                    ttATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGCTTTTGCGGCGTTGCTGGCTGCCGGACTCATTCGCAGt                <  4:93580/149‑1 (MQ=255)
                                                                                                                                     tATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGCTTTTGCGGCGTTGCTGGCTGCCGGACTCATTCGCAGTg               <  3:15855/149‑1 (MQ=255)
                                                                                                                                       ttCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGCTTTTGCGGCGTTGCTGGCTGCCGGACTCATTCGCAGTGGa             >  2:608120/1‑149 (MQ=255)
                                                                                                                                        tCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGCTTTTGCGGCGTTGCTGGCTGCCGGACTCATTCGCAGTGGAc            >  4:301250/1‑149 (MQ=255)
                                                                                                                                          ttGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGc                                                     <  2:543414/106‑1 (MQ=255)
                                                                                                                                          ttGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGc                                                     >  1:543414/1‑106 (MQ=255)
                                                                                                                                          ttGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGCTTTTGCGGCGTTGCTGGCTGCCGGACTCATTCGCAGTGGACaa          >  4:156411/1‑149 (MQ=255)
                                                                                                                                              tttAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGCTTTTGCGGCGTTGCTGGCTGCCGGACTCATTCGCAGTGGACAAAttt      >  2:316521/1‑149 (MQ=255)
                                                                                                                                              tttAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGCTTTTGCGGCGTTGCTGGCTGCCGGACTCATTCGCAGTGGACAAAttt      >  3:371112/1‑149 (MQ=255)
                                                                                                                                                  ggTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGCTTTTGCGGCGTTGCTGGCTGCCGGACTCATTCGCAGTGGACAAATTTTTCt  <  1:316521/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
ATGCCATTGCCCAGCATCAGCCGCTCACCGGCCGTTTTCTGCAGGCAAGTGTCATTGAACCGAATTGTAGCCTGCGTGCACAGTATCTGGCGTCACTACAAACAACGCTGCAACCCTATCCGTGGAATGCGGTTATTCTTGGTTTAGGTATGGGGGAACGGTTATCCGTCCCCAAAGAAATCAAAAATATCATGCGCGATACTGGAACGGCGCATTTAATGGCGATATCGGGATTGCACATCGCTTTTGCGGCGTTGCTGGCTGCCGGACTCATTCGCAGTGGACAAATTTTTCT  >  NC_000913/964774‑965068

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