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.3% 289.3 / ‑0.1 81G348G (GGC→GGGemrDmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base G (42/36);  total (43/38)
Fisher's exact test for biased strand distribution p-value = 5.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

GGACGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAGGGGCAGCCCGTTTAAGCGC  >  NC_000913/3854820‑3855110
                                                                                                                                                 |                                                                                                                                                 
ggACGCTGCTCTTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGt                                                                                                                                                 <  1:995084/148‑1 (MQ=255)
ggACGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGt                                                                                                                                                 <  1:593317/148‑1 (MQ=255)
   cGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGc                                                                                                                                           <  2:1199818/151‑1 (MQ=255)
    gttccgcGTGCCCGCCGCGCTGTTCTTTTTCTGTGCCGGGATGCTGATTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCGTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTTCTGGGTGCTGGc                                                                                                                                           <  2:738415/144‑1 (MQ=255)
    gctgctCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGc                                                                                                                                           <  1:886420/150‑1 (MQ=255)
       gctCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTc                                                                                                                                        <  1:41376/150‑1 (MQ=255)
        ctCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTgg                                                                                                                                            <  2:1044338/145‑1 (MQ=255)
         tCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTc                                                                                                                                        <  2:872821/148‑1 (MQ=255)
          cGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGt                                                                                                                                         <  2:160996/146‑1 (MQ=255)
           gTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGGGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGct                                                                                                                                     <  1:615254/149‑1 (MQ=255)
           gTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGctc                                                                                                                                    <  1:709249/150‑1 (MQ=255)
           gTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGGCGGCGGTCTGCAAAACATTGGTTCCGGGGGGCTGGCGTCGc                                                                                                                                      >  2:831666/1‑148 (MQ=255)
            ttCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGctctc                                                                                                                                  >  1:810634/1‑151 (MQ=255)
            ttCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGct                                                                                                                                     <  2:315690/148‑1 (MQ=255)
               ccgccgCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCCAAACATTGGTTCCGGGGGGCTGGCGTCGCTCTCTg                                                                                                                                >  2:499419/1‑150 (MQ=255)
               ccgccgCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGCCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCCGGCGTCGctctc                                                                                                                                  >  1:117350/1‑148 (MQ=255)
                cgccgcGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCg                                                                                                                              >  1:589708/1‑151 (MQ=255)
                 gccgcGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGGTCCGGGGGGCTGGCGTCGCTCTCTg                                                                                                                                >  2:159687/1‑148 (MQ=255)
                   cgcgCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATg                                                                                                                           <  1:464490/151‑1 (MQ=255)
                      gcTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGGCGGCGGGCTGCAAAAAATTGGTTCCGGGGTGCTGGCGGCGCTCTCTGCGATGTTg                                                                                                                        >  2:58442/1‑151 (MQ=255)
                        tGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGcc                                                                                                                      >  1:664648/1‑151 (MQ=255)
                        tGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGcc                                                                                                                      >  1:356013/1‑151 (MQ=255)
                        tGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGGGCTGGCGTCGCCCTCTGCGATGTTGcc                                                                                                                      >  2:1051602/1‑151 (MQ=255)
                        tGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGCCTGCAAAACATTGGTTCCGGGGGGCTGGCGTCGCTCTCTGCGATGtt                                                                                                                         >  2:380075/1‑148 (MQ=255)
                          ttCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGCCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGc                                                                                                                    >  1:397098/1‑151 (MQ=255)
                          ttCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGGCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGcc                                                                                                                      >  2:140436/1‑149 (MQ=255)
                           tCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGc                                                                                                                    <  1:69053/150‑1 (MQ=255)
                               tttCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAAcc                                                                                                               >  1:14471/1‑151 (MQ=255)
                                 tCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCg                                                                                                              <  2:588735/150‑1 (MQ=255)
                                   ggTGCCGGGATGCTGTTTCCGCTGGCGACCATCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTc                                                                                                           >  1:210223/1‑151 (MQ=255)
                                   ggTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTc                                                                                                           <  1:1165558/151‑1 (MQ=255)
                                      gCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAg                                                                                                         <  1:302395/150‑1 (MQ=255)
                                         gggATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCa                                                                                                     <  2:246624/151‑1 (MQ=255)
                                          ggATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGc                                                                                                      >  1:1099536/1‑149 (MQ=255)
                                              gCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCt                                                                                                 <  2:79702/150‑1 (MQ=255)
                                              gCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCt                                                                                                 <  1:201435/150‑1 (MQ=255)
                                                tGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGGGCTGGCGTCGCTCTCTGCGATGTTGGCGCAAACCGGTCCAGGCAGCCTggg                                                                                              >  2:703753/1‑151 (MQ=255)
                                                 gTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCACAAACCGGTCAGGGCAGCCTggg                                                                                              <  1:319231/150‑1 (MQ=255)
                                                     ccGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGgttgt                                                                                         <  1:897827/151‑1 (MQ=255)
                                                      cGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGgttgtt                                                                                        <  2:1064784/151‑1 (MQ=255)
                                                            cGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGAc                                                                                  <  2:774630/151‑1 (MQ=255)
                                                                cAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGGTGTTGATGACCtt                                                                               >  2:1209459/1‑150 (MQ=255)
                                                                 aGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGGCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGGTGATGACCt                                                                                >  2:733760/1‑148 (MQ=255)
                                                                  gcggcgCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAAt                                                                            <  1:1099725/151‑1 (MQ=255)
                                                                       gcgATGGAGCCGTTCCCCGTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAGCCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGa                                                                        <  2:988788/150‑1 (MQ=255)
                                                                        cgATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGAtt                                                                      >  1:245724/1‑151 (MQ=255)
                                                                          aTGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGGTGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGAttgt                                                                    >  2:767468/1‑151 (MQ=255)
                                                                          aTGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGAttgt                                                                    >  1:1026989/1‑151 (MQ=255)
                                                                               gCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATc                                                               <  1:842424/151‑1 (MQ=255)
                                                                                       ccTTCCTGGCGGGCACGGCTGGCGCGCCGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCgtgctgtg                                                       >  2:865245/1‑151 (MQ=255)
                                                                                       ccTTCCTGGCGGGCACGGCCGGCGCGCTGGGCGGCGGGCTGCAAAACAATGGGTTCGGGGTGCTGGCGTAGACCTCCGCGATGTTGCCGCAAACCGGTAAGGGGAGCCGGGGGTTGTTGATGACCTTAATGGGCGAGGTGATCgtgctgt                                                        >  2:745097/1‑150 (MQ=255)
                                                                                           ccTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTgctgg                                                   >  1:901683/1‑151 (MQ=255)
                                                                                             tGGCGGGCACGGTTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCCCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGGTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTgctggc                                                  >  2:625420/1‑150 (MQ=255)
                                                                                             tGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTgctgg                                                   >  1:418724/1‑149 (MQ=255)
                                                                                              ggCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGGTCCGGGGTGCTGGCGCCGCCCCCCGCGATGTTGCCGCAAACCGGTAAGGGCAACCTGGGGGTGTTGATGACCTTGATGGGATTGTTTTTCgtgcctcggtggct                                                 >  2:90957/1‑140 (MQ=255)
                                                                                                 gggCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCg                                             >  1:236994/1‑151 (MQ=255)
                                                                                                 gggCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCg                                             <  1:431264/151‑1 (MQ=255)
                                                                                                          tGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGAc                                     >  1:49101/1‑150 (MQ=255)
                                                                                                              gcgcTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGGTGACCTTAATGGGATTGTTGATCGGGCTGTGCTGGCTGCCGCTGGCGACGCgg                                 >  2:692585/1‑150 (MQ=255)
                                                                                                               cgcTGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGAt                               <  2:1104925/151‑1 (MQ=255)
                                                                                                               cgcTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGa                                <  2:700424/150‑1 (MQ=255)
                                                                                                                  tGGTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTc                            >  2:561008/1‑151 (MQ=255)
                                                                                                                   ggTCGGCGGTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTc                            <  2:182871/150‑1 (MQ=255)
                                                                                                                         cggTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATc                       >  2:1057908/1‑149 (MQ=255)
                                                                                                                          ggTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAgg                    >  2:1042469/1‑151 (MQ=255)
                                                                                                                           gTCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAggg                   >  1:1110836/1‑151 (MQ=255)
                                                                                                                            tCTGCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAgggg                  >  1:44998/1‑151 (MQ=255)
                                                                                                                               gCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAgggg                  <  2:396841/148‑1 (MQ=255)
                                                                                                                               gCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAGGGGc                 >  2:372767/1‑149 (MQ=255)
                                                                                                                               gCAAAACATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAGGGGc                 <  2:554783/149‑1 (MQ=255)
                                                                                                                                 aaaaCATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAGGGGCAGcc             >  1:890528/1‑151 (MQ=255)
                                                                                                                                  aaaCATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAGGGGCAGccc            >  2:915927/1‑151 (MQ=255)
                                                                                                                                    aCATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGAGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAGGGGCAGCCCGt          >  2:793614/1‑151 (MQ=255)
                                                                                                                                     cATTGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAGGGGCAGCCCGt          >  2:1055813/1‑150 (MQ=255)
                                                                                                                                       ttGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAGGGGCAGCCCGttt        >  1:879784/1‑150 (MQ=255)
                                                                                                                                       ttGGTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAGGGGCAGCCCGttt        <  2:845801/150‑1 (MQ=255)
                                                                                                                                         ggTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACGGGTCATGGCAGGCTGGGGTTGTTGGTGACCTTAATGGGATTGTTGGTCGTGCTGTGCTGGCGGCCGCTGGCGACGCGGATGTCACATCAGGGGCAGACCGTTTaa      <  2:1135037/150‑1 (MQ=255)
                                                                                                                                         ggTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAGGGGCAGCCCGTTTAAg     <  1:1043792/151‑1 (MQ=255)
                                                                                                                                          gTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAGGGGCAGCCCGTTTAAgc    <  1:1186718/151‑1 (MQ=255)
                                                                                                                                          gTTCCGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCCGCCGCTGGCGACGCGGATGTCGCATCAGGGGCAGCCCGTTTAAgc    >  1:1039096/1‑151 (MQ=255)
                                                                                                                                             ccGGGGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAGGGGCAGCCCGTTTAAgcgc  <  2:648231/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GGACGCTGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCATCAGGGGCAGCCCGTTTAAGCGC  >  NC_000913/3854820‑3855110

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