Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,000,046 G→A R114C (CGC→TGC)  fliY ← cystine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,000,0460GA98.7% 293.9 / NA 77R114C (CGC→TGC) fliYcystine transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (42/34);  total (42/35)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTTCATAGTTGGTGCCCAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCAAATTCCACTTCAAAAC  >  NC_000913/1999901‑2000192
                                                                                                                                                 |                                                                                                                                                  
tctTCATAGTTGGTGCCCAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACtcat                                                                                                                                               >  2:884623/1‑151 (MQ=255)
    cATAGTTGGTGCCCAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAg                                                                                                                                            <  1:165920/150‑1 (MQ=255)
       aGTTGGTGCCCAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAAt                                                                                                                                        <  2:380387/151‑1 (MQ=255)
        gTTGGTGCCCAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATg                                                                                                                                       >  2:354402/1‑151 (MQ=255)
          tGGTGCCCAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATgg                                                                                                                                      <  2:73059/150‑1 (MQ=255)
          tGGTGCCCAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATgg                                                                                                                                      <  2:43071/150‑1 (MQ=255)
           ggTGCCCAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGt                                                                                                                                     <  2:1244754/150‑1 (MQ=255)
              gCCCAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGGGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGAcc                                                                                                                                 >  1:1266637/1‑151 (MQ=255)
                ccAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCt                                                                                                                                >  1:661441/1‑150 (MQ=255)
                ccAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTg                                                                                                                               <  2:366915/151‑1 (MQ=255)
                ccAGACCGACCCCCACTTTTTTGACTTTCAGATCATCGGCTGCTTTAAAGGTGCCTTCGTTACCTTTTTTTACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCt                                                                                                                                >  1:1009687/1‑150 (MQ=255)
                 cAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTg                                                                                                                               >  2:1151045/1‑150 (MQ=255)
                 cAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTg                                                                                                                               >  2:284859/1‑150 (MQ=255)
                  aGACCGACCCCCACTTTTTTGCCTTTCAGATCGTCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGAt                                                                                                                             <  2:393188/151‑1 (MQ=255)
                     ccgaccCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTaa                                                                                                                          <  2:948201/151‑1 (MQ=255)
                      cgaccCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAAt                                                                                                                         <  1:54345/151‑1 (MQ=255)
                       gaccCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATc                                                                                                                        >  2:1286296/1‑151 (MQ=255)
                         cccccACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATcac                                                                                                                      >  2:662622/1‑151 (MQ=255)
                         cccccACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATcac                                                                                                                      >  2:224779/1‑151 (MQ=255)
                            ccACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATcaccac                                                                                                                   <  1:772825/151‑1 (MQ=255)
                             cACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACcaca                                                                                                                  <  1:929241/151‑1 (MQ=255)
                              aCTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACAt                                                                                                                 <  2:276983/151‑1 (MQ=255)
                               cTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATc                                                                                                                >  2:138480/1‑151 (MQ=255)
                                 tttttGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATc                                                                                                                <  1:349413/149‑1 (MQ=255)
                                      gCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGt                                                                                                         <  2:981034/151‑1 (MQ=255)
                                       ccTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGtt                                                                                                        <  2:670941/151‑1 (MQ=255)
                                           tCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGttt                                                                                                       >  2:560102/1‑148 (MQ=255)
                                               atcatcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTcc                                                                                                >  2:753911/1‑151 (MQ=255)
                                                tcatcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGGGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTcc                                                                                                >  1:551479/1‑150 (MQ=255)
                                                 catcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAg                                                                                              <  2:1160169/151‑1 (MQ=255)
                                                  atcGGCTGGTTTAAAGGAGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGGGTACGGGGTTGAGGAATCGTATTTTTTCTTGCACTCAACAGAAATGGTGACCTGTTTAATCACCACATCAATACGTTTAGAGTCCAg                                                                                              >  1:752102/1‑150 (MQ=255)
                                                       cTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGc                                                                                        >  1:183737/1‑151 (MQ=255)
                                                          ttttAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGcc                                                                                       >  2:1093827/1‑149 (MQ=255)
                                                           attAATGGAGCCTTCGTTACCTTTTTTAACCAGCGCCTGAATACCCGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAg                                                                                     >  1:875889/2‑150 (MQ=255)
                                                            ttAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGc                                                                                    >  2:1276169/1‑150 (MQ=255)
                                                            ttAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACAAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCa                                                                                      <  1:893843/148‑1 (MQ=255)
                                                             tAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCa                                                                                   <  1:1151011/150‑1 (MQ=255)
                                                                 ggTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCg                                                                              >  2:414326/1‑151 (MQ=255)
                                                                  gTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGt                                                                             >  2:50401/1‑151 (MQ=255)
                                                                     ccTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTccc                                                                          >  1:241240/1‑151 (MQ=255)
                                                                      cTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTccc                                                                          >  2:900418/1‑150 (MQ=255)
                                                                        tCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCAtt                                                                       >  2:104953/1‑151 (MQ=255)
                                                                         cGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCAttt                                                                      <  1:495707/151‑1 (MQ=255)
                                                                          gTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTg                                                                     >  2:785888/1‑151 (MQ=255)
                                                                          gTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTg                                                                     <  2:327149/151‑1 (MQ=255)
                                                                            tACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGt                                                                   >  1:1205875/1‑151 (MQ=255)
                                                                              ccTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCg                                                                 >  2:947645/1‑151 (MQ=255)
                                                                               cTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCgg                                                                >  1:104073/1‑151 (MQ=255)
                                                                               cTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCgg                                                                >  2:747576/1‑151 (MQ=255)
                                                                                    ttCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGtttg                                                            <  2:554016/150‑2 (MQ=255)
                                                                                          aGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACg                                                     >  1:184504/1‑151 (MQ=255)
                                                                                              ccTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTc                                                 >  1:1186126/1‑151 (MQ=255)
                                                                                                 gAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCaa                                               <  1:297611/150‑1 (MQ=255)
                                                                                                  aaTACCCGAAATGGTGTACGGGGTTGAGCAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATAAGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGACTCAACg                                             >  1:751249/1‑151 (MQ=255)
                                                                                                     aCCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGcc                                           <  2:282292/150‑1 (MQ=255)
                                                                                                        aGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCa                                          >  1:5550/1‑148 (MQ=255)
                                                                                                        aGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCa                                          >  2:926642/1‑148 (MQ=255)
                                                                                                            aTGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGt                                    >  2:1269485/1‑150 (MQ=255)
                                                                                                              ggTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGtttt                                 <  2:830925/151‑1 (MQ=255)
                                                                                                                tgtACGGGGTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGc                               >  1:994022/1‑151 (MQ=255)
                                                                                                                   aCGGGGTTGATAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAg                            <  1:472608/151‑1 (MQ=255)
                                                                                                                      gggTTGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTg                         <  2:889319/151‑1 (MQ=255)
                                                                                                                         ttGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTg                      >  1:1111907/1‑151 (MQ=255)
                                                                                                                          tGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTg                      <  2:346983/150‑1 (MQ=255)
                                                                                                                          tGAGAAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTg                      <  2:793686/150‑1 (MQ=255)
                                                                                                                           gagaAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTggg                    <  2:912734/151‑1 (MQ=255)
                                                                                                                             gaAATCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGc                   >  2:843414/1‑150 (MQ=255)
                                                                                                                              aaaTCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCaa                 <  2:210804/151‑1 (MQ=255)
                                                                                                                              aaaTCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCaa                 <  2:259359/151‑1 (MQ=255)
                                                                                                                                aTCGTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCAAAt               >  2:287969/1‑151 (MQ=255)
                                                                                                                                   gTATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCAAATTcc            >  2:163148/1‑151 (MQ=255)
                                                                                                                                    tATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCACGATGTTTTGCCAGCTGTTGGGCAAATTcc            >  1:55747/1‑150 (MQ=255)
                                                                                                                                    tATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCAAATTcc            >  2:40003/1‑150 (MQ=255)
                                                                                                                                    tATTTTTTCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCAAATTcc            <  1:54702/150‑1 (MQ=255)
                                                                                                                                       tttttCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCAAATTCCACt         <  1:997565/150‑1 (MQ=255)
                                                                                                                                       tttttCTTGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCAAATTCCACt         <  2:420293/150‑1 (MQ=255)
                                                                                                                                              tGCACTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCAAATTCCACTTCAAAAc  <  1:409189/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TCTTCATAGTTGGTGCCCAGACCGACCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACGCCAGCATACCGTCCCATTTGGTCGGTTTTAGTGACGCCTCAACGCCAAGATGTTTTGCCAGCTGTTGGGCAAATTCCACTTCAAAAC  >  NC_000913/1999901‑2000192

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