Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,854,965 C→G G348G (GGC→GGG emrD → multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,854,9650CG96.4% 295.1 / 1.4 83G348G (GGC→GGGemrDmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (3/0);  new base G (41/39);  total (44/39)
Fisher's exact test for biased strand distribution p-value = 2.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

CTGGACGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAGGGGCAGCCCGTTTAAGC  >  NC_000913/3854818‑3855108
                                                                                                                                                   |                                                                                                                                               
cTGGACGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTg                                                                                                                                              >  1:550185/1‑151 (MQ=255)
  ggACGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGt                                                                                                                                               <  1:670991/148‑1 (MQ=255)
  ggACGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTg                                                                                                                                              <  1:290753/149‑1 (MQ=255)
  ggACGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGc                                                                                                                                             <  2:237599/150‑1 (MQ=255)
  ggACGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGc                                                                                                                                             <  1:965275/150‑1 (MQ=255)
     cGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGGCTGCAAAACATTGGTTCCGGGGt                                                                                                                                               >  2:740405/1‑145 (MQ=255)
      gctgctCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTgg                                                                                                                                          >  1:1074932/1‑149 (MQ=255)
      gctgctCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGc                                                                                                                                         >  2:1142033/1‑150 (MQ=255)
         gctCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGGCGGCGGTCTGCAAAACATTTGGTCCGGGGGGCTGGCGt                                                                                                                                       >  2:691831/1‑149 (MQ=255)
          ctCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTc                                                                                                                                      <  1:1039650/149‑1 (MQ=255)
           tCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTc                                                                                                                                      <  2:836570/148‑1 (MQ=255)
           tCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGct                                                                                                                                   >  2:770715/1‑151 (MQ=255)
               tCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGctctc                                                                                                                                <  2:545659/150‑1 (MQ=255)
               tCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACCGCTGGGGCGCTGGTCGGCGGGCTGCACAACAATGGTTCCCGGGTGCTGGCGTCGctctc                                                                                                                                >  2:350164/1‑150 (MQ=255)
                     cgcgCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGGGCTGGCGTCGCTCTCTGCgct                                                                                                                          >  1:255204/1‑148 (MQ=255)
                      gcgcTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGGTCCGGGGGGCTGGCGTCGCTCTCTGCGAt                                                                                                                          >  1:576212/1‑149 (MQ=255)
                      gcgcTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGGCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATg                                                                                                                         >  1:1269237/1‑150 (MQ=255)
                       cgcTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGGCGCTCTCTGCGATGtt                                                                                                                       >  1:354918/1‑151 (MQ=255)
                          tGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGGGCTGGCGTCGCTCTCTGCGATGTTGc                                                                                                                     >  1:315695/1‑150 (MQ=255)
                           gTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCg                                                                                                                   <  1:1076783/151‑1 (MQ=255)
                                ttttCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGGGCTGGCGTCGCTCTTTGCGATGTTGCCGCaaa                                                                                                               >  1:807088/1‑150 (MQ=255)
                                   tCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACCTTGGTTCCGGGGGGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCgg                                                                                                           >  1:276976/1‑151 (MQ=255)
                                          cGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAggg                                                                                                     >  1:81676/1‑150 (MQ=255)
                                           gggATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGGCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGc                                                                                                    >  1:542662/1‑150 (MQ=255)
                                            ggATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGGGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAg                                                                                                  <  2:634967/151‑1 (MQ=255)
                                            ggATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCa                                                                                                   <  2:436850/150‑1 (MQ=255)
                                            ggATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAg                                                                                                  <  2:775147/151‑1 (MQ=255)
                                               tGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGcc                                                                                                <  1:1261839/150‑1 (MQ=255)
                                                  tGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCCGGGGCGCTGGTCGGCCGGCTGCAAAACATTGGTTCCGGGGTGCTGGCGCCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGGCGGCTggg                                                                                            >  1:802736/1‑151 (MQ=255)
                                                   gTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTgg                                                                                             <  1:784900/149‑1 (MQ=255)
                                                      tCCGCTGGCGCCCCGCGGCGCGCTGGCGCCGTTCCCCTTCCTGGCGGGCCCTTCTTGCGCGCTGGTCGGCGGTCTGCAACACCTTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGggtg                                                                                        >  1:655034/1‑151 (MQ=255)
                                                       ccGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGgttgt                                                                                       <  2:401785/151‑1 (MQ=255)
                                                       ccGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGgttg                                                                                        <  1:1296980/150‑1 (MQ=255)
                                                         gCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGttgttg                                                                                     <  2:129781/151‑1 (MQ=255)
                                                           tGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGGGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCCGGGCAGCCTGGGGgtgttg                                                                                     >  1:669873/1‑149 (MQ=255)
                                                           tGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGGCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGttgttggt                                                                                   >  1:225229/1‑149 (MQ=255)
                                                            ggCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTtga                                                                                    <  1:601839/149‑1 (MQ=255)
                                                             gCGACCAGCGGCGCGATGGAGCCGTTCCCCGTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTtgatga                                                                                 <  2:709029/151‑1 (MQ=255)
                                                                 ccAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCtt                                                                             <  1:772307/151‑1 (MQ=255)
                                                                  cAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCtt                                                                             >  1:177112/1‑150 (MQ=255)
                                                                  cAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCtt                                                                             >  2:210583/1‑150 (MQ=255)
                                                                   aGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTaa                                                                           <  2:1104366/151‑1 (MQ=255)
                                                                       gcgcGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATggg                                                                       >  2:1048109/1‑151 (MQ=255)
                                                                          cgATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGAtt                                                                    >  2:1086109/1‑151 (MQ=255)
                                                                            aTGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGAttgt                                                                  <  1:887288/151‑1 (MQ=255)
                                                                            aTGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCCGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGAttg                                                                   >  2:54484/1‑150 (MQ=255)
                                                                                 gCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGAt                                                              <  1:169620/150‑1 (MQ=255)
                                                                                  ccGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGTTCAGGGCAGCCTTCGGTTGTTGATGCCCTTAAGGGGTTTGTTGa                                                               <  1:1199859/148‑1 (MQ=255)
                                                                                  ccGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCg                                                            <  2:1247275/151‑1 (MQ=255)
                                                                                   cGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCg                                                            >  1:795433/1‑150 (MQ=255)
                                                                                     ttCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCgtg                                                          >  1:431042/1‑150 (MQ=255)
                                                                                        cccTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCgtgctgt                                                      >  2:1267073/1‑151 (MQ=255)
                                                                                          cTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGGTGTTGATGACCTTAATGGGATTGTTGATCgtgctgtgc                                                    >  2:691562/1‑151 (MQ=255)
                                                                                           ttCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGGGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGGTGTTGATGACCTTAATGGGATTGTTGATCgtgctgtg                                                     >  1:366373/1‑149 (MQ=255)
                                                                                            tCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCgtgctgtgc                                                    <  1:507696/149‑1 (MQ=255)
                                                                                             ccTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCgtgctgtgcta                                                  <  2:612499/150‑2 (MQ=255)
                                                                                             ccTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGtgctgtgctg                                                  <  2:107616/150‑1 (MQ=255)
                                                                                             ccTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGtgctgtgctg                                                  <  2:84324/150‑1 (MQ=255)
                                                                                             ccTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTgctgg                                                 <  1:33221/151‑1 (MQ=255)
                                                                                               tGGCGGGCACGGCTGGCGCGCTGGTCGGCGGGCTGCAAAACATTGGTTCCGGGGGGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCCGCCTGGGGGTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTgctggct                                               >  2:277050/1‑151 (MQ=255)
                                                                                                    ggCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGc                                          <  2:763315/151‑1 (MQ=255)
                                                                                                     gCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCt                                         <  2:1249925/151‑1 (MQ=255)
                                                                                                        cGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGc                                      <  2:1259034/151‑1 (MQ=255)
                                                                                                            tGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGAcg                                  >  1:1259414/1‑151 (MQ=255)
                                                                                                             ggCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGAcgc                                 >  2:808268/1‑151 (MQ=255)
                                                                                                             ggCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGAcg                                  <  1:923991/150‑1 (MQ=255)
                                                                                                                  gcTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGAt                             <  1:356630/150‑1 (MQ=255)
                                                                                                                    tGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTc                          >  1:663843/1‑151 (MQ=255)
                                                                                                                     ggTCGGCGTTCTGCAAAACATTGGTTCCGGGGTGCGGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGGTGGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTc                          <  1:263429/150‑1 (MQ=255)
                                                                                                                     ggTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCg                         <  1:740378/151‑1 (MQ=255)
                                                                                                                       tCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGc                        >  1:346396/1‑150 (MQ=255)
                                                                                                                         ggcggTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATc                     >  1:401351/1‑151 (MQ=255)
                                                                                                                          gcggTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCa                    <  1:231113/151‑1 (MQ=255)
                                                                                                                                 gCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAGGGGCAg             <  2:1116552/151‑1 (MQ=255)
                                                                                                                                  cAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAGGGGCAGc            >  1:1089078/1‑151 (MQ=255)
                                                                                                                                   aaaaCATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAGGGGCAGcc           >  2:900337/1‑151 (MQ=255)
                                                                                                                                    aaaCATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAGGGGCAGccc          >  1:320087/1‑151 (MQ=255)
                                                                                                                                        aTTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAGGGGCAGCCCGttt      >  2:320116/1‑151 (MQ=255)
                                                                                                                                        aTTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAGGGGCAGCCCGttt      >  1:404191/1‑151 (MQ=255)
                                                                                                                                        aTTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAGGGGCAGCCCGttt      <  1:1253237/151‑1 (MQ=255)
                                                                                                                                        aTTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAGGGGCAGCCCGttt      >  1:1136731/1‑151 (MQ=255)
                                                                                                                                            gTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAGGGGCAGCCCGTTTAAgc  >  2:616475/1‑151 (MQ=255)
                                                                                                                                            gTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAGGGGCAGCCCGTTTAAgc  >  2:237545/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
CTGGACGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAGGGGCAGCCCGTTTAAGC  >  NC_000913/3854818‑3855108

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