Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,014,960 G→A 5.3% K27K (AAG→AAA) ‡ 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,9600GA 5.3% 364.8 / 5.5 112K27K (AAG→AAA) ‡fliGflagellar motor switching and energizing component
Reads supporting (aligned to +/- strand):  ref base G (48/58);  new base A (3/3);  total (51/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 113 total aligned reads displayed.

TCCGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCAAGAAGCTGAACAGTTTGCCGCACTGAATATCAACGCCAACGATTATC  >  NC_000913/2014811‑2015108
                                                                                                                                                     |                                                                                                                                                    
tCCGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAg                                                                                                                                                      <  2:927348/150‑1 (MQ=255)
 ccGTGAAATGTCTGATAACGATCCGCGCGTGGGGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGc                                                                                                                                                     <  1:172795/150‑1 (MQ=255)
   gTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGGCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCAc                                                                                                                                                   <  2:923871/150‑1 (MQ=255)
    tGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCAcc                                                                                                                                                  <  2:1380238/150‑1 (MQ=255)
    tGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCAcc                                                                                                                                                  <  1:418470/150‑1 (MQ=255)
         tGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGCTTGAACACctctgc                                                                                                                                             <  1:1452014/150‑3 (MQ=255)
         tGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGAAAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTcc                                                                                                                                             <  1:145173/150‑1 (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)
                      cccGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTACCCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTAc                                                                                                                                <  2:277808/149‑1 (MQ=255)
                        cgcgcgTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACaa                                                                                                                              >  1:1186303/1‑150 (MQ=255)
                        cgcgcgTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCGGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACaa                                                                                                                              <  1:266487/150‑1 (MQ=255)
                          cgcgTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAAc                                                                                                                            <  1:1645830/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)
                              gggtggCGCTGGGCATTCGCCAGTGGATAAATAACGCGCATGACTAAGCTGACAGGCACCGATAACACCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTg                                                                                                                        <  2:1390417/149‑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:799106/150‑1 (MQ=255)
                                           cATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATgg                                                                                                           >  2:121761/1‑150 (MQ=255)
                                            aTTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGc                                                                                                          >  1:416052/1‑150 (MQ=255)
                                              tCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGa                                                                                                        >  2:1602631/1‑150 (MQ=255)
                                                   aGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTc                                                                                                   <  1:648239/150‑1 (MQ=255)
                                                    gTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCa                                                                                                  <  2:1186303/150‑1 (MQ=255)
                                                     tGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTccc                                                                                                 >  2:689632/1‑148 (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                                                                                        >  2:221284/1‑150 (MQ=255)
                                                              aaCGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGAtctc                                                                                        <  1:1635942/150‑1 (MQ=255)
                                                               aCGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTcc                                                                                       <  2:70165/150‑1 (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)
                                                                          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)
                                                                                              aaaGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCAGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGt                                                        >  2:1049334/1‑150 (MQ=255)
                                                                                                    tCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGc                                                  <  2:1198832/150‑1 (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)
                                                                                                               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)
                                                                                                                         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)
                                                                                                                                ccGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCAAGAAGCTGAACAGTTTGCCGCACTGaa                      <  1:1602631/150‑1 (MQ=255)
                                                                                                                                    ctgggagagGTGCTTGAACACCTCTGCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCAAGAAGCTGAACAGTTTGCCGCACTg                        <  2:1649673/139‑1 (MQ=255)
                                                                                                                                    ctgggagagGTGCTTGAACACCTCTGCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCAAGAAGCTGAACAGTTTGCCGCACTg                        >  1:1649673/6‑144 (MQ=255)
                                                                                                                                     tgggagagGTGCTTGAACACCTCTGCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCAAGAAGCTGAACAGTTTGCCGCACTGAATATCa                 >  1:996913/5‑150 (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:1538446/1‑150 (MQ=255)
                                                                                                                                            ggTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCAAGAAGCTGAACAGTTTGCCGCACTGAATATCAACGCCaa          >  1:479299/1‑150 (MQ=255)
                                                                                                                                              tgtTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTGTGAGCAAGAAGCTGAACAGTTTGCCGCACTGAATATCAACGCCAACg        >  2:1033789/1‑150 (MQ=255)
                                                                                                                                              ggtTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCAAGAAGCTGAACAGTTTGCCGCACTGAATATCAACGCCAACg        <  1:600419/149‑1 (MQ=255)
                                                                                                                                                ttCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCAAGAAGCTGAACAGTTTGCCGCACTGAATATCAACGCCAACGAt      <  1:1050685/150‑1 (MQ=255)
                                                                                                                                                  cAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCAAGAAGCTGAACAGTTTGCCGCACTGAATATCAACGCCAACGATTa    >  2:574501/1‑150 (MQ=255)
                                                                                                                                                  cAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCAAGAAGCTGAACAGTTTGCCGCACTGAATATCAACGCCAACGATTa    >  2:699769/1‑150 (MQ=255)
                                                                                                                                                    aGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCAAGAAGCTGAACAGTTTGCCGCACTGAATATCAACGCCAACGATTATc  <  1:1183607/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCCGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACCTCTCCCAGCGTGAAGTACAAACCCTGAGCGCTGCAATGGCGAACGTCACGCAGATCTCCAACAAGCAGCTAACCGATGTGCTGGCGGAGTTTGAGCAAGAAGCTGAACAGTTTGCCGCACTGAATATCAACGCCAACGATTATC  >  NC_000913/2014811‑2015108

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: