Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,014,958 A→G 5.4% K27E (AAG→GAG) ‡ fliG → flagellar motor switching and energizing component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,014,9580AG 5.4% 335.0 / 6.9 111K27E (AAG→GAG) ‡fliGflagellar motor switching and energizing component
Reads supporting (aligned to +/- strand):  ref base A (48/57);  new base G (3/3);  total (51/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

Only 100 of 111 total aligned reads displayed.

TCCGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCAAGAAGCTGAACAGTTTGCCGCACTGAATATCAACGCCAACGATTA  >  NC_000913/2014811‑2015106
                                                                                                                                                   |                                                                                                                                                    
 ccGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGc                                                                                                                                                   <  2:1572955/150‑1 (MQ=255)
 ccGTGAAATGTCTGATAACGATCCGCGCGTGGGGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGc                                                                                                                                                   <  1:172795/150‑1 (MQ=255)
   gTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGGCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCAc                                                                                                                                                 <  2:923871/150‑1 (MQ=255)
    tGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCAcc                                                                                                                                                <  1:418470/150‑1 (MQ=255)
    tGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCAcc                                                                                                                                                <  2:1380238/150‑1 (MQ=255)
         tGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTcc                                                                                                                                           >  2:422850/1‑150 (MQ=255)
         tGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGCTTGAACACctctgc                                                                                                                                           <  1:1452014/150‑3 (MQ=255)
          gTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTccc                                                                                                                                          <  1:924257/150‑1 (MQ=255)
            cTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGCTTGAACACCTCTgccgc                                                                                                                                        <  1:1256502/150‑3 (MQ=255)
              gATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCg                                                                                                                                      <  1:1427337/150‑1 (MQ=255)
                tAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTg                                                                                                                                    <  2:1236152/150‑1 (MQ=255)
                 aaCGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCACAGCGTGa                                                                                                                                   >  2:1285924/1‑150 (MQ=255)
                   cGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGCCCGGGCGGCAGAGGTGTTAAAGCACCTCTCCCAGCGTGAAg                                                                                                                                 >  2:1050685/1‑150 (MQ=255)
                   cGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAg                                                                                                                                 <  1:1446333/150‑1 (MQ=255)
                   cGATCCGCGCGGGGTGGCGCTGGGCATTCGCCAGTGGATAAATAACGATCATGAGTACCCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAg                                                                                                                                 <  2:347928/150‑1 (MQ=255)
                   cGATCAGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCTTGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCATGCACCTCTCCCAGCGTGAAg                                                                                                                                 <  2:491985/150‑1 (MQ=255)
                     aTCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTa                                                                                                                               <  1:700303/150‑1 (MQ=255)
                     aTCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGAAGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTa                                                                                                                               <  1:1474643/150‑1 (MQ=255)
                      tCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTAc                                                                                                                              <  2:499273/150‑1 (MQ=255)
                        cgcgcgTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCGGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACaa                                                                                                                            <  1:266487/150‑1 (MQ=255)
                           gcgTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAAcc                                                                                                                         <  1:1217459/150‑1 (MQ=255)
                           gcgTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAAcc                                                                                                                         <  2:1542038/150‑1 (MQ=255)
                             gtggtgGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCt                                                                                                                       <  1:812146/150‑1 (MQ=255)
                              tgggggCGCTGGTCAATCGCCAGTGGATTAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTg                                                                                                                      <  1:124481/150‑1 (MQ=255)
                               ggtggCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTgc                                                                                                                     >  2:1201501/1‑149 (MQ=255)
                                    cgcTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCt                                                                                                                >  2:1656306/1‑150 (MQ=255)
                                          tCATTCGCCAGTGGATAAATAACGATCATGAGTAAACTGACACGCACCGATATAAGCGGCATCCTGCTGATGAGCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATg                                                                                                          <  1:1254116/150‑1 (MQ=255)
                                          gcATTCGCCCGGGGCTAACTCACGCTCATGAGTAACCTGCCAGGCACCGATAAAAGCGTCCTCCTGCTGCTGACCATTGGCGAAGACCGGGCGGCAGAGGGGTTCAAGCACCTCTCCCAGCGTGAAGTCCAAACCCTGAGCGCTGCAATg                                                                                                          <  2:135859/149‑1 (MQ=255)
                                           cATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGCGCGCTGCAATgg                                                                                                         >  2:1069388/1‑150 (MQ=255)
                                           cATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATgg                                                                                                         >  2:121761/1‑150 (MQ=255)
                                           cATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATgg                                                                                                         <  2:799106/150‑1 (MQ=255)
                                           cATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATgg                                                                                                         <  1:422850/150‑1 (MQ=255)
                                              tCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGa                                                                                                      >  2:1602631/1‑150 (MQ=255)
                                                   aGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTc                                                                                                 <  1:648239/150‑1 (MQ=255)
                                                    gTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCa                                                                                                <  2:1186303/150‑1 (MQ=255)
                                                      ggATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACg                                                                                              >  1:176597/1‑150 (MQ=255)
                                                       gATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGc                                                                                             >  1:350174/1‑150 (MQ=255)
                                                       gATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTAAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGc                                                                                             >  1:1087538/1‑150 (MQ=255)
                                                        ataaataaCGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCa                                                                                            >  1:1425547/1‑150 (MQ=255)
                                                         taaataaCGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAg                                                                                           >  1:1462211/1‑150 (MQ=255)
                                                          aaataaCGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGa                                                                                          <  1:527267/150‑1 (MQ=255)
                                                              aaCGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGAtctc                                                                                      <  1:1635942/150‑1 (MQ=255)
                                                              aaCGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGAtctc                                                                                      >  2:221284/1‑150 (MQ=255)
                                                               aCGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTcc                                                                                     <  2:70165/150‑1 (MQ=255)
                                                                cGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCca                                                                                    >  1:1669293/1‑150 (MQ=255)
                                                                cGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCca                                                                                    <  1:1271408/150‑1 (MQ=255)
                                                                 gATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCcaa                                                                                   >  2:348281/1‑150 (MQ=255)
                                                                  aTCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAACACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGCGCGCTGCAATGGCGAACGTCCCGCAGATCGCcaac                                                                                  >  2:1025185/1‑150 (MQ=255)
                                                                   tCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCcaaca                                                                                 <  1:44537/150‑1 (MQ=255)
                                                                      tGAGTAACCTGACAGGCACCGATAAAAGCGTAATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTGCCAGCGTGAAGTACAAATCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAagc                                                                              <  2:279353/150‑1 (MQ=255)
                                                                       gAGTAACCTGACAGGCCCCGATAAACGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGTGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAagca                                                                             <  1:62950/150‑1 (MQ=255)
                                                                         gTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAagcagc                                                                           >  2:729939/1‑150 (MQ=255)
                                                                          tacccggACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCt                                                                          <  2:1560739/144‑1 (MQ=255)
                                                                            aCCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTaa                                                                        >  1:72183/1‑150 (MQ=255)
                                                                            aCCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGACCGCTGAAATGGCGAACGTCACGCAGATCTCCAACAAGCATCTaa                                                                        >  1:1148273/1‑150 (MQ=255)
                                                                                  cAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGAtg                                                                  <  1:1128259/150‑1 (MQ=255)
                                                                                  cAGGCAACGATAAAATCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCACGCACCTCTCCCAGCGGGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGAtg                                                                  <  1:1119011/150‑1 (MQ=255)
                                                                                   aGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGAtgt                                                                 <  2:778655/150‑1 (MQ=255)
                                                                                   aGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGAtgt                                                                 <  1:700175/150‑1 (MQ=255)
                                                                                      cACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCt                                                              <  1:976972/150‑1 (MQ=255)
                                                                                             aaaaGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACACACCCTGAGCGCGGAAATTGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGCTGTGCTGgaggag                                                       >  2:77055/1‑150 (MQ=255)
                                                                                              aaaGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGt                                                      >  2:1183607/1‑150 (MQ=255)
                                                                                              aaaGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCAGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGt                                                      >  2:1049334/1‑150 (MQ=255)
                                                                                                    tCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTAGCGGAGTTTGAGc                                                <  2:1402457/150‑1 (MQ=255)
                                                                                                     cATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAACCCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCa                                               >  2:782634/1‑150 (MQ=255)
                                                                                                     cATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCa                                               >  1:1085140/1‑150 (MQ=255)
                                                                                                     cATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCa                                               >  2:156300/1‑150 (MQ=255)
                                                                                                         ctgctgATGACCATTGGCGAAGCCCGGGCGGCAGAGGTGTTCAAGCCCCTCTCCCAGCGGGAAGTCCAACCCCTGAGCGCTGCACGGGCGAACGTCACGGCGATCTCCCACAAGCCGCTAACCGATGTGCGGGAGGAGTTTGAGCaagac                                           >  1:224903/1‑149 (MQ=255)
                                                                                                         ctgctgATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCaagaa                                           >  1:590775/1‑150 (MQ=255)
                                                                                                          tgctgATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCaagaag                                          >  1:528970/1‑150 (MQ=255)
                                                                                                           gctgATGACCATTGGCGACGACCGGGCGGCAGAGGTGCTCAAGCACCTCTCCCAGCGTGAAGAACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCAAGAAGc                                         <  1:246703/150‑1 (MQ=255)
                                                                                                             tgatgaCCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCAAGAAGCTg                                       <  2:129103/150‑1 (MQ=255)
                                                                                                             tgatgaCCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCAAGAAGCTg                                       <  1:1575753/150‑1 (MQ=255)
                                                                                                               atgaCCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCAAGAAGCTGaa                                     <  2:762015/150‑1 (MQ=255)
                                                                                                                tgaCCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCAAGAAGCTGAAc                                    >  2:1028394/1‑150 (MQ=255)
                                                                                                                   ccATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCAAGAAGCTGAACAGt                                 >  2:1377092/1‑150 (MQ=255)
                                                                                                                      ttGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCAAGAAGCTGAACAGTTTg                              <  1:1656306/150‑1 (MQ=255)
                                                                                                                        ggCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCAAGAAGCTGAACAGTTTGcc                            <  1:1606026/150‑1 (MQ=255)
                                                                                                                        ggCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCAAGAAGCTGAACAGTTTGcc                            <  1:159252/150‑1 (MQ=255)
                                                                                                                         gCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGATCGTCACGCAGATCTCCATCAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCAAGAAGCTGAACAGTTTGCCg                           >  1:1282833/1‑150 (MQ=255)
                                                                                                                         gCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAGCGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCAAGAAGCTGAACAGTTTGCCg                           >  2:675514/1‑150 (MQ=255)
                                                                                                                         gCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCAAGAAGCTGAACAGTTTGCCg                           <  1:552084/150‑1 (MQ=255)
                                                                                                                          cGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTTCAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCAAGAAGCTGAACAGTTTGCCGc                          >  1:1072791/1‑150 (MQ=255)
                                                                                                                          cGAAGACAGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACACTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGCTGTGCTGGCGGAGTTTGAGCAAGAAGCTGAACAGTTTGCCGc                          >  1:108416/1‑150 (MQ=255)
                                                                                                                             aGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCAAGAAGCTGAACAGTTTGCCGCACt                       >  1:483320/1‑150 (MQ=255)
                                                                                                                              gACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCAAGAAGCTGAACAGTTTGCCGCACTg                      >  2:1366801/1‑150 (MQ=255)
                                                                                                                                ccGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCAAGAAGCTGAACAGTTTGCCGCACTGaa                    <  1:1602631/150‑1 (MQ=255)
                                                                                                                                    ctgggagagGTGCTTGAACACCTCTGCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCAAGAAGCTGAACAGTTTGCCGCACTg                      >  1:1649673/6‑144 (MQ=255)
                                                                                                                                    ctgggagagGTGCTTGAACACCTCTGCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCAAGAAGCTGAACAGTTTGCCGCACTg                      <  2:1649673/139‑1 (MQ=255)
                                                                                                                                       gcAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCAAGAAGCTGAACAGTTTGCCGCACTGAATATCAAc             <  1:1075143/150‑1 (MQ=255)
                                                                                                                                        gagagGTGCTTGAACACCTCTGCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCAAGAAGCTGAACAGTTTGCCGCACTGAATATCAACg            >  2:1149793/2‑150 (MQ=255)
                                                                                                                                            ggTGTTCAAGCACCTCTCCCAGTGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCAAGAAGCTGAACAGTTTCCCGCACTGAATATCAACGCCaa        >  2:491482/1‑150 (MQ=255)
                                                                                                                                            ggTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCAAGAAGCTGAACAGTTTGCCGCACTGAATATCAACGCCaa        >  1:479299/1‑150 (MQ=255)
                                                                                                                                            ggTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCAAGAAGCTGAACAGTTTGCCGCACTGAATATCAACGCCaa        >  1:1538446/1‑150 (MQ=255)
                                                                                                                                              tgtTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTGTGAGCAAGAAGCTGAACAGTTTGCCGCACTGAATATCAACGCCAACg      >  2:1033789/1‑150 (MQ=255)
                                                                                                                                               gtTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCAAGAAGCTGAACAGTTTGCCGCACTGAATATCAACGCCAACGa     >  1:113726/1‑150 (MQ=255)
                                                                                                                                               gtTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCAAGAAGCTGAACAGTTTGCCGCACTGAATATCAACGCCAACGa     >  1:990591/1‑150 (MQ=255)
                                                                                                                                                ttCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCAAGAAGCTGAACAGTTTGCCGCACTGAATATCAACGCCAACGAt    <  1:1050685/150‑1 (MQ=255)
                                                                                                                                                  cAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCAAGAAGCTGAACAGTTTGCCGCACTGAATATCAACGCCAACGATTa  >  2:574501/1‑150 (MQ=255)
                                                                                                                                                  cAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCAAGAAGCTGAACAGTTTGCCGCACTGAATATCAACGCCAACGATTa  >  2:699769/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TCCGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCAAGAAGCTGAACAGTTTGCCGCACTGAATATCAACGCCAACGATTA  >  NC_000913/2014811‑2015106

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