Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,166,156 G→A 5.3% P521P (CCC→CCT)  ygiS ← putative ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,166,1560GA 5.3% 359.2 / 4.7 112P521P (CCC→CCT) ygiSputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (60/46);  new base A (4/2);  total (64/48)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.64e-01

Only 100 of 119 total aligned reads displayed.

CAAAGGGGCAATATTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCTTTGTCGCATCAGTGATCTGCGT  >  NC_000913/3166007‑3166305
                                                                                                                                                     |                                                                                                                                                      
cAAAGGGGCAATATTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTgg                                                                                                                                                        >  2:65435/1‑150 (MQ=255)
 aaaGGGGCAATATTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTggg                                                                                                                                                       <  1:1671072/150‑1 (MQ=255)
  aaGGGGCAATATTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTgggg                                                                                                                                                      <  2:17419/150‑1 (MQ=255)
  aaGGGGCAATATTGCCGGAGATACTACATTATCTGCCGCGTCAAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTgggg                                                                                                                                                      <  2:1578688/150‑1 (MQ=255)
     gggCAATATTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGAtt                                                                                                                                                   <  2:643895/150‑1 (MQ=255)
      ggCAATATTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTa                                                                                                                                                  >  2:1426915/1‑150 (MQ=255)
         aaTATTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCATTGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGc                                                                                                                                               <  2:1303778/150‑1 (MQ=255)
         aaTATTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGc                                                                                                                                               <  1:706278/150‑1 (MQ=255)
           tatTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAg                                                                                                                                             <  2:886758/150‑1 (MQ=255)
           tatTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCCCAGGATTATgtct                                                                                                                                             <  1:1703602/150‑4 (MQ=255)
           tatTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCCCAGGATTATgtct                                                                                                                                             <  1:1710288/150‑4 (MQ=255)
            atTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGc                                                                                                                                            <  1:457593/150‑1 (MQ=255)
                ccGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGaa                                                                                                                                        <  1:220016/150‑1 (MQ=255)
                  ggAGATACTACATTATCTGCCGCGTCGAGTCACAGGATATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACTGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGaaaa                                                                                                                                      <  2:1130635/150‑1 (MQ=255)
                    agaTACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATATGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAAcc                                                                                                                                    >  1:1218686/1‑150 (MQ=255)
                       tactacATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGcc                                                                                                                                 >  1:1601674/1‑150 (MQ=255)
                             aTTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTa                                                                                                                           <  2:752474/150‑1 (MQ=255)
                             aTTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTa                                                                                                                           >  1:426991/1‑150 (MQ=255)
                             aTTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTa                                                                                                                           >  1:616693/1‑150 (MQ=255)
                              ttATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGCCATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑g                                                                                                                         >  2:62343/1‑150 (MQ=255)
                                aTCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑ggg                                                                                                                       <  1:554648/150‑1 (MQ=255)
                                 tCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACACCTCTTTGCTGTAGACATAATCCTGGGGATTATGAAGCGGAAAACCGCCAACGTA‑GGGt                                                                                                                      >  1:1106921/1‑150 (MQ=255)
                                  cTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGtt                                                                                                                     >  1:1483931/1‑150 (MQ=255)
                                   tGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGttt                                                                                                                    <  2:89156/150‑1 (MQ=255)
                                    tccGCTTCGAGTCACAGGATATCAACGTGAATAGCGCCGCCGGACTTGCTCCGACGGCGATAAAGGCATAAATTTGCCTTGATATACAACTATTTGCTGTAGAAATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTa                                                                                                                   <  1:1047353/149‑1 (MQ=255)
                                       gcgTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGc                                                                                                                <  1:342501/150‑1 (MQ=255)
                                       gcgTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCGCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGc                                                                                                                <  1:354894/150‑1 (MQ=255)
                                           cGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCCGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGtt                                                                                                            >  2:1085788/1‑150 (MQ=255)
                                            gAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGttt                                                                                                           <  2:1062838/150‑1 (MQ=255)
                                            gAGGCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGttt                                                                                                           <  2:1489489/150‑1 (MQ=255)
                                             aGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTg                                                                                                          <  1:1613612/150‑1 (MQ=255)
                                              gTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGa                                                                                                         >  1:1449604/1‑150 (MQ=255)
                                              gTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGa                                                                                                         >  1:1165833/1‑150 (MQ=255)
                                               tCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGAt                                                                                                        <  1:1170119/150‑1 (MQ=255)
                                                cacaGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGAtt                                                                                                       >  2:362456/1‑150 (MQ=255)
                                                   aGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAAc                                                                                                    <  1:1154447/150‑1 (MQ=255)
                                                       gaTCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCt                                                                                                >  1:477788/1‑150 (MQ=255)
                                                       gaTCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCt                                                                                                <  1:655653/150‑1 (MQ=255)
                                                        ctCAACGTGTATTCCGCCTCCGGAAGTTCTCCTACGTAGATAAAGGCATAAATGTGCCTTGATATCCAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTg                                                                                               <  2:1140021/149‑1 (MQ=255)
                                                          cAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGAt                                                                                             >  2:843376/1‑150 (MQ=255)
                                                           aaCGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCCGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGAta                                                                                            >  2:375624/1‑150 (MQ=255)
                                                           aaCGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGAta                                                                                            >  1:955488/1‑150 (MQ=255)
                                                            aCGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGAtag                                                                                           >  1:1143838/1‑150 (MQ=255)
                                                               tGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGAtagtag                                                                                        <  2:1725274/150‑1 (MQ=255)
                                                               tGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGGA‑GGGTTTAAGCAGTTTGATTAACGGCTGAtagtag                                                                                        >  1:671470/1‑150 (MQ=255)
                                                                gAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGa                                                                                       >  2:1115026/1‑150 (MQ=255)
                                                                    aGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAgg                                                                                   >  2:275248/1‑150 (MQ=255)
                                                                      cgccgcCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACAAAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGaa                                                                                 <  1:445251/150‑1 (MQ=255)
                                                                       gccgccGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAAt                                                                                >  2:126634/1‑150 (MQ=255)
                                                                          gccgGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAg                                                                             >  1:818365/1‑150 (MQ=255)
                                                                           ccgGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGc                                                                            <  1:65435/150‑1 (MQ=255)
                                                                            cgGACGCGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGAGAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCg                                                                           >  2:871429/1‑150 (MQ=255)
                                                                                cGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGc                                                                       >  2:195742/1‑150 (MQ=255)
                                                                                 gTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGcc                                                                      >  2:922140/1‑150 (MQ=255)
                                                                                 gTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGcc                                                                      <  1:843376/150‑1 (MQ=255)
                                                                                  tGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCt                                                                     >  1:602944/1‑150 (MQ=255)
                                                                                   gCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTg                                                                    >  1:1535242/1‑150 (MQ=255)
                                                                                    cTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGAAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTgt                                                                   >  1:943280/1‑150 (MQ=255)
                                                                                     tCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTgtt                                                                  >  2:164859/1‑150 (MQ=255)
                                                                                      ccGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTAGGGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTgtt                                                                  >  1:1494290/1‑150 (MQ=255)
                                                                                      ccGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTgttg                                                                 <  1:1641375/150‑1 (MQ=255)
                                                                                        gACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTgttggt                                                               >  2:1638483/1‑150 (MQ=255)
                                                                                          cGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTgttggttg                                                             >  1:1209378/1‑150 (MQ=255)
                                                                                            gCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTtgat                                                           <  1:86667/150‑1 (MQ=255)
                                                                                            gCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTtgat                                                           <  2:798851/150‑1 (MQ=255)
                                                                                            gCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGAAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTtgat                                                           <  2:75662/150‑1 (MQ=255)
                                                                                             cGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCCAGTTA‑GGGTTTAAGCAGTTTGATTAACGTCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTtgatg                                                          >  1:127160/1‑150 (MQ=255)
                                                                                             cGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTtgatg                                                          >  1:874461/1‑150 (MQ=255)
                                                                                             cGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTtgatg                                                          <  2:489251/150‑1 (MQ=255)
                                                                                             cGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGACAACGTA‑GGGTTTAAGCAGTTTGATAAACGGCTGATAGTTGATAGTAATCAGCGTTGCCTGTTGGTtgatg                                                          >  2:1640639/1‑150 (MQ=255)
                                                                                                tAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATc                                                       <  2:448870/150‑1 (MQ=255)
                                                                                                 aaaGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCa                                                      <  1:1365587/150‑1 (MQ=255)
                                                                                                      cATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAATCCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTTATAGTAGATAGGAATCAGCGGTTCCTGTTGGTTGATGTTCACTTct                                                 >  2:78531/1‑150 (MQ=255)
                                                                                                         cAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTctgcc                                              >  1:560595/1‑150 (MQ=255)
                                                                                                         cAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTctgcc                                              <  2:813743/150‑1 (MQ=255)
                                                                                                           aTGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGAGTTTGCAGCGGAAAACCGCCAACGTA‑GGGTGTAAGGAGTTTGAGTACCGGCTGATTGTAGATAGGAATCAGAGGTGGATGTTGGGTGATGATCCCTTctgcctg                                            >  2:272860/1‑150 (MQ=255)
                                                                                                                ccTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGAt                                       >  1:738577/1‑150 (MQ=255)
                                                                                                                 cTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATa                                      >  1:1441503/1‑150 (MQ=255)
                                                                                                                    gATATACAACTCTTTGCTGTAGACATAATCCTGGGGAGTATGCAGCGGAAAACCGCAAACGTA‑GGGTTTAAGCAGTTTGATTAAGGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACaa                                   >  2:669680/1‑150 (MQ=255)
                                                                                                                      tataCAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACg                                 >  1:338389/1‑150 (MQ=255)
                                                                                                                      tataCAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACg                                 >  2:654031/1‑150 (MQ=255)
                                                                                                                          cAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTACCGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCAGATACAACGCAtt                             >  2:299667/1‑150 (MQ=255)
                                                                                                                          cAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCAtt                             >  2:738384/1‑150 (MQ=255)
                                                                                                                              tctTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGGTTAAGCCGTTTGAGTAACGGCTGAGTGTAGAGAGGACTCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGc                         >  2:172056/1‑150 (MQ=255)
                                                                                                                               ctTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCt                        >  2:635099/1‑150 (MQ=255)
                                                                                                                                tttGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCtt                       <  2:204125/150‑1 (MQ=255)
                                                                                                                                tttGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCtt                       <  2:1449604/150‑1 (MQ=255)
                                                                                                                                 ttGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCttt                      <  2:426991/150‑1 (MQ=255)
                                                                                                                                    cTGTAGACATAATCCTGGTGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCTTTGTc                   <  1:458476/150‑1 (MQ=255)
                                                                                                                                      gTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGACAACGTA‑GGCTTTAAGCAGTTTGAATAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAAAGCATTACGCTTTGTCGc                 >  2:1214736/1‑150 (MQ=255)
                                                                                                                                      cgctgcATAATCCCCAGGATTATGTAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCTTTGTCGc                 >  1:1187529/6‑150 (MQ=255)
                                                                                                                                       gctgcATAATCCCCAGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCTTTGTCGCa                >  1:1128885/5‑150 (MQ=255)
                                                                                                                                       gctgcATAATCCCCAGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCTTTGTCGCa                >  1:1521216/5‑150 (MQ=255)
                                                                                                                                          aCATAATCCTGGGGATTATGCAGCGGAAAACCTCCAACGTA‑GGGTTTAATCAGTTTGATTAACGGCTGATAGTAGATAGGTATCAGCGGTGCCTGTTTGTTGATGATCACTTCTGCCGGCTGATTCAACGCATTACGCTTTGTCGCATCa             >  2:210005/1‑150 (MQ=255)
                                                                                                                                          aCATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCCCTTCTGCCTGCTGATACAACGCATTACGCTTTGTCGCATCa             >  1:1355284/1‑150 (MQ=255)
                                                                                                                                          aCATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCTTTGTCGCATCa             >  1:202770/1‑150 (MQ=255)
                                                                                                                                             tAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCTTTGTCGCATCAGTg          >  2:57532/1‑150 (MQ=255)
                                                                                                                                               aTCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTTCCTGCTGATACAACGCATTACGCTTTGTCGCATCAGTGAt        >  1:113191/1‑150 (MQ=255)
                                                                                                                                                    ggggATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACGACGCATTACGCTTTGTCGCATCAGTGATCTGCg   >  1:755477/1‑150 (MQ=255)
                                                                                                                                                     gggATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGACTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCTTTGTCGCATCAGTGATCTGCGt  >  1:284581/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CAAAGGGGCAATATTGCCGGAGATACTACATTATCTGCCGCGTCGAGTCACAGGAGATCAACGTGAATAGCGCCGCCGGACGTGCTCCGACGGCGATAAAGGCATCAATGTGCCTTGATATACAACTCTTTGCTGTAGACATAATCCTGGGGATTATGCAGCGGAAAACCGCCAACGTA‑GGGTTTAAGCAGTTTGATTAACGGCTGATAGTAGATAGGAATCAGCGGTGCCTGTTGGTTGATGATCACTTCTGCCTGCTGATACAACGCATTACGCTTTGTCGCATCAGTGATCTGCGT  >  NC_000913/3166007‑3166305

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑