Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,805,131 +T coding (317/1203 nt) proV → glycine betaine/proline ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,805,1311.T97.2% 272.7 / ‑0.7 72coding (317/1203 nt)proVglycine betaine/proline ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base T (23/47);  total (25/47)
Fisher's exact test for biased strand distribution p-value = 1.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGTTCGGTATGGAATTGGCCGGAATTAATGCCGAAGAACGCCGGGAAAAAGCCCTTGATGCACTGCGTCAGGTCGGGC  >  NC_000913/2804982‑2805274
                                                                                                                                                       |                                                                                                                                               
atTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGaaat                                                                                                                                                 <  2:1166531/151‑2 (MQ=255)
 tttGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGaaat                                                                                                                                                 >  2:385284/1‑149 (MQ=255)
    gTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGaaataaag                                                                                                                                            <  2:1059874/151‑3 (MQ=255)
      catcatGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGa                                                                                                                                            <  2:961060/150‑1 (MQ=255)
      catcatGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGAt                                                                                                                                           <  1:814396/151‑1 (MQ=255)
       atcatGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGAtt                                                                                                                                          <  2:551425/151‑1 (MQ=255)
         catGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGc                                                                                                                                        <  2:952764/151‑1 (MQ=255)
         catGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGc                                                                                                                                        >  2:713435/1‑151 (MQ=255)
          atGGGATTATCCGGCTCGGGTTAATCCACAATGGTATGCCTTCTCAATGGCCTGATTGAACCCACCCGCGGGCGAGTGCTGAGTGATGGAGTGGATATTGCCAAAATATCCGACGCCTAACTATGTGAGGTGCGCAGAAATAAgtattg                                                                                                                                        >  1:383902/1‑144 (MQ=255)
            gggATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGAt                                                                                                                                     <  2:678848/151‑1 (MQ=255)
              gATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATgg                                                                                                                                   <  2:954712/151‑1 (MQ=255)
                    ccGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTc                                                                                                                              <  2:1052451/150‑1 (MQ=255)
                    ccGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTc                                                                                                                              <  1:638484/150‑1 (MQ=255)
                     cGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTAATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCa                                                                                                                            <  2:25713/151‑1 (MQ=255)
                      ggCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTcc                                                                                                                             >  2:603442/1‑149 (MQ=255)
                        cTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTc                                                                                                                         >  1:445482/1‑151 (MQ=255)
                          cGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTc                                                                                                                         <  1:1050656/149‑1 (MQ=255)
                          cGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTc                                                                                                                         <  2:605239/149‑1 (MQ=255)
                            ggTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCttt                                                                                                                     <  1:656836/151‑1 (MQ=255)
                            ggTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCttt                                                                                                                     <  2:111431/151‑1 (MQ=255)
                                aaTCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCt                                                                                                                 <  2:260895/151‑1 (MQ=255)
                                  tCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCtt                                                                                                                >  1:901040/1‑150 (MQ=255)
                                    cacaATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAAt                                                                                                             <  2:146595/151‑1 (MQ=255)
                                      caATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGc                                                                                                           >  1:637688/1‑151 (MQ=255)
                                       aaTGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGcc                                                                                                          <  2:341506/151‑1 (MQ=255)
                                        aTGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCg                                                                                                         <  2:276799/151‑1 (MQ=255)
                                           gTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCat                                                                                                      >  2:36955/1‑151 (MQ=255)
                                               gCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGa                                                                                                  <  1:1140850/151‑1 (MQ=255)
                                                   tctcAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGAcc                                                                                                <  2:1116409/149‑1 (MQ=255)
                                                      cAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGc                                                                                            >  2:316865/1‑150 (MQ=255)
                                                          cGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGAc                                                                                       >  1:766138/1‑151 (MQ=255)
                                                             cTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAAt                                                                                    >  1:63211/1‑151 (MQ=255)
                                                              tGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAAt                                                                                    >  1:557831/1‑150 (MQ=255)
                                                               gATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATa                                                                                   <  2:42177/150‑1 (MQ=255)
                                                               gATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATa                                                                                   <  2:854058/150‑1 (MQ=255)
                                                                  tGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTg                                                                                <  1:10905/150‑1 (MQ=255)
                                                                   gAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCg                                                                              <  2:261885/151‑1 (MQ=255)
                                                                      cccacccGCGGGCAAGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGTTc                                                                           <  1:1053018/151‑1 (MQ=255)
                                                                           ccGCGGGCAAGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGTTCGGTAt                                                                      >  2:565990/1‑151 (MQ=255)
                                                                             gcgGGCAAGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGTTCGGTATgg                                                                    <  1:917255/151‑1 (MQ=255)
                                                                              cgGGCAAGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGTTCGGTATGGa                                                                   <  2:855568/151‑1 (MQ=255)
                                                                               gggCAAGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGTTCGGTATGGaa                                                                  <  1:331404/151‑1 (MQ=255)
                                                                                   aaGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGTTCGGTATGGAATTgg                                                              <  2:628518/151‑1 (MQ=255)
                                                                                       gCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGTTCGGTATGGAATTGGCCg                                                           <  2:21657/150‑1 (MQ=255)
                                                                                           attgatGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGTTCGGTATGGAATTGGCCGGAAtt                                                      <  2:535062/151‑1 (MQ=255)
                                                                                              gatGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGTTCGGTATGGAATTGGCCGGAATTAAt                                                   >  1:1182007/1‑151 (MQ=255)
                                                                                              gatGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGTTCGGTATGGAATTGGCCGGAATTAAt                                                   >  1:582112/1‑151 (MQ=255)
                                                                                              gatGGTGTGGATATTGCCAAAATATCAGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGTTCGGTATGGAATTGGCCGGAATTaa                                                    <  2:1032451/150‑1 (MQ=255)
                                                                                               atGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGTTCGGTATGGAATTGGCCGGAATTAATg                                                  <  2:1184031/151‑1 (MQ=255)
                                                                                                  gtgtGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGTTCGGTATGGAATTGGCCGGAATTAATGCCg                                               <  2:264368/151‑1 (MQ=255)
                                                                                                     tgGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGTTCGGTATGGAATTGGCCGGAATTAATGCCgaa                                             <  1:859437/150‑1 (MQ=255)
                                                                                                      ggATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGTTCGGTATGGAATTGGCCGGAATTAATGCCgaag                                            <  1:152506/150‑1 (MQ=255)
                                                                                                        atatTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGTTCGGTATGGAATTGGCCGGAATTAATGCCGAAGAAc                                         >  2:342453/1‑151 (MQ=255)
                                                                                                          atTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGTTCGGTATGGAATTGGCCGGAATTAATGCCGAAGAACGc                                       >  1:835201/1‑151 (MQ=255)
                                                                                                          atTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGTTCGGTATGGAATTGGCCGGAATTAATGCCGAAGAACGc                                       >  1:766222/1‑151 (MQ=255)
                                                                                                                aaaaTATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGTTCGGTATGGAATTGGCCGGAATTAATGCCGAAGAACGCCGGGaa                                 <  1:849405/151‑1 (MQ=255)
                                                                                                                 aaaTATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGTTCGGTATGGAATTGGCCGGAATTAATGCCGAAGAACGCCGGGaaa                                <  1:1120152/151‑1 (MQ=255)
                                                                                                                  aaTATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGTTCGGTATGGAATTGGCCGGAATTAATGCCGAAGAACGCCGGGaaaa                               <  2:890089/151‑1 (MQ=255)
                                                                                                                  aaTATCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGTTCGGTATGGAATTGGCCGGAATTAATGCCGAAGAACGCCGGGaaa                                >  1:383915/1‑150 (MQ=255)
                                                                                                                   atatCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGTTCGGTATGGAATTGGCCGGAATTAATGCCGAAGAACGCCGGGaaaaa                              <  1:911712/151‑1 (MQ=255)
                                                                                                                     atCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGTTCGGTATGGAATTGGCCGGAATTAATGCCGAAGAACGCCGGGAAAAAGc                            >  1:491066/1‑151 (MQ=255)
                                                                                                                      tCCGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGTTCGGTATGGAATTGGCCGGAATTAATGCCGAAGAACGCCGGGAAAAAGcc                           >  2:862102/1‑151 (MQ=255)
                                                                                                                       ccGACGCCGAACTCCGTGAGGTGCGCAGACATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGTTCGGTATGGAATTGGCCGGAATTAATGCCGAAGAACGCCGGGAAAAAGccc                          >  1:1098749/1‑151 (MQ=255)
                                                                                                                        cGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGTTCGGTATGGAATTGGCCGGAATTAATGCCGAAGAACGCCGGGAAAAAGCCCt                         >  2:465817/1‑151 (MQ=255)
                                                                                                                        cGACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGTTCGGTATGGAATTGGCCGGAATTAATGCCGAAGAACGCCGGGAAAAAGCCCt                         >  2:178497/1‑151 (MQ=255)
                                                                                                                         gACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGTTCGGTATGGAATTGGCCGGAATTAATGCCGAAGAACGCCGGGAAAAAGCCCtt                        >  1:243621/1‑151 (MQ=255)
                                                                                                                         gACGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGTTCGGTATGGAATTGGCCGGAATTAATGCCGAAGAACGCCGGGAAAAAGCCCtt                        <  2:253031/151‑1 (MQ=255)
                                                                                                                           cGCCGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGTTCGGTATGGAATTGGCCGGAATTAATGCCGAAGAACGCCGGGAAAAAGCCCTTGa                      <  2:603926/151‑1 (MQ=255)
                                                                                                                             ccGAACTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGTTCGGTATGGAATTGGCCGGAATTAATGCCGAAGAACGCCGGGAAAAAGCCCTTGATg                    <  1:655774/151‑1 (MQ=255)
                                                                                                                                aaCTCCGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGTTCGGTATGGAATTGGCCGGAATTAATGCCGAAGAACGCCGGGAAAAAGCCCTTGATGCAc                 <  2:186267/151‑1 (MQ=255)
                                                                                                                                    ccGTGAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGTTCGGTATGGAATTGGCCGGAATTAATGCCGAAGAACGCCGGGAAAAAGCCCTTGATGCACTg               <  1:1133422/149‑1 (MQ=255)
                                                                                                                                        gAGGTGCGCAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGTTCGGTATGGAATTGGCCGGAATTAATGCCGAAGAACGCCGGGAAAAAGCCCTTGATGCACTGCGTCa          <  1:603412/150‑1 (MQ=255)
                                                                                                                                             gcgcAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGTTCGGTATGGAATTGGCCGGAATTAATGCCGAAGAACGCCGGGAAAAAGCCCTTGATGCACTGCGTCAGGTCg     <  1:152451/150‑1 (MQ=255)
                                                                                                                                               gcAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGTTCGGTATGGAATTGGCCGGAATTAATGCCGAAGAACGCCGGGAAAAAGCCCTTGATGCACTGCGTCAGGTCGGGc  >  2:697942/1‑151 (MQ=255)
                                                                                                                                               gcAGAAATAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGTTCGGTATGGAATTGGCCGGAATTAATGCCGAAGAACGCCGGGAAAAAGCCCTTGATGCACTGCGTCAGGTCGGGc  <  2:474072/151‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                               
ATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGTTCGGTATGGAATTGGCCGGAATTAATGCCGAAGAACGCCGGGAAAAAGCCCTTGATGCACTGCGTCAGGTCGGGC  >  NC_000913/2804982‑2805274

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