Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,160,028 C→T P548L (CCG→CTG)  mdtC → multidrug efflux system, subunit C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,160,0280CT98.6% 278.0 / NA 72P548L (CCG→CTG) mdtCmultidrug efflux system, subunit C
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (34/37);  minor base G (1/0);  total (35/37)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGACTGCGTGGTTTTGGTCGCATGTTGGTAGCCCTGCAACAAGGCTACGGCAAGTCACTAAAATGGGTGCTCAATCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCCGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACGATCCGGCAGTGGATAATGTCACCGG  >  NC_000913/2159881‑2160175
                                                                                                                                                   |                                                                                                                                                   
aGCGACTGCGTGGTTTTGGTCGCATGTTGGTAGCCCTGCAACAAGGCTACGGCAAGTCACTAAAATGGGTGCTCAATCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGaa                                                                                                                                                  <  2:604954/151‑1 (MQ=255)
       gCGTGGTTTTGGTCGCATGTTGGTAGCCCTGCAACAAGGCTACGGCAAGTCACTAAAATGGGTGCTCAATCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACct                                                                                                                                             <  2:518396/149‑1 (MQ=255)
        cGTGGTTTTGGTCGCATGTTGGTAGCCCTGCAACAAGGCTACGGCAAGTCACTAAAATGGGTGCTCAATCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACcttct                                                                                                                                          <  2:1185876/151‑1 (MQ=255)
        cGTGGTTTTGGTCGCATGTTGGTAGCCCTGCAACAAGGCTACGGCAAGTCACTAAAATGGGTGCTCAATCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACcttct                                                                                                                                          >  1:528809/1‑151 (MQ=255)
              tttGGTCGCATGTTGGTAGCCCTGCAACAAGGCTACGGCAAGTCACTAAAATGGGTGCTCAATCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCg                                                                                                                                     <  1:7036/150‑1 (MQ=255)
               ttGGTCGCATGTTGGTAGCCCTGCAACAAGGCTACGGCAAGTCACTAAAATGGGTGCTCAATCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCgg                                                                                                                                    >  2:901692/1‑150 (MQ=255)
               ttGGTCGCATGTTGGTAGCCCTGCAACAAGGCTACGGCAAGTCACTAAAATGGGTGCTCAATCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCgg                                                                                                                                    <  1:926412/150‑1 (MQ=255)
                 ggTCGCATGTTGGTAGCCCTGCAACAAGGCTACGGCAAGTCACTAAAATGGGTGCTCAATCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAg                                                                                                                                  <  1:1116824/150‑1 (MQ=255)
                  gTCGCATGTTGGTAGCCCTGCAACAAGGCTACGGCAAGTCACTAAAATGGGTGCTCAATCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCa                                                                                                                                <  2:206494/151‑1 (MQ=255)
                   tCGCATGTTGGTAGCCCTGCAACAAGGCTACGGCAAGTCACTAAAATGGGTGCTCAATCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGc                                                                                                                                 >  1:766630/1‑149 (MQ=255)
                    cGCATGTTGGTAGCCCTGCAACAAGGCTACGGCAAGTCACTAAAATGGGTGCTCAATCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAg                                                                                                                               <  1:751707/150‑1 (MQ=255)
                     gCATGTTGGTAGCCCTGCAACAAGGCTACGGCAAGTCACTAAAATGGGTGCTCAATCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAgg                                                                                                                              >  1:273617/1‑150 (MQ=255)
                      cATGTTGGTAGCCCTGCAACAAGGCTACGGCAAGTCACTAAAATGGGTGCTCAATCATACCCGTCTGGTTGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGac                                                                                                                            <  2:265104/151‑1 (MQ=255)
                        tGTTGGTAGCCCTGCAACAAGGCTACGGCAAGTCACTAAAATGGGTGCTCAATCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGacac                                                                                                                          >  1:130225/1‑151 (MQ=255)
                           tGGTAGCCCTGCAACAAGGCTACGGCAAGTCACTAAAATGGGTGCTCAATCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTgg                                                                                                                       >  1:43736/1‑151 (MQ=255)
                            ggTAGCCCTGCAACAAGGCTACGGCAAGTCACTAAAATGGGTGCTCAATCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGc                                                                                                                      <  2:202903/151‑1 (MQ=255)
                                  ccTGCAACAAGGCTACGGCAAGTCACTAAAATGGGTGCTCAATCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTg                                                                                                                <  1:141978/151‑1 (MQ=255)
                                        acaaGGCTACGGCAAGTCACTAAAATGGGTGCTCAATCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGc                                                                                                          <  1:316460/151‑1 (MQ=255)
                                         caaGGCTACGGCAAGTCACTAAAATGGGTGCTCAATCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCg                                                                                                         <  1:963887/151‑1 (MQ=255)
                                                aCGGCAAGTCACTAAAATGGGTGCTCAATCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTc                                                                                                   >  1:48784/1‑150 (MQ=255)
                                                  ggCAAGTCACTAAAATGGGTGCTCAATCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAgg                                                                                                <  2:477788/151‑1 (MQ=255)
                                                    cAAGTCACTAAAATGGGTGCTCAATCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCg                                                                                              >  2:402536/1‑151 (MQ=255)
                                                     aaGTCACTAAAATGGGTGCTCAATCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCgg                                                                                             >  1:208066/1‑151 (MQ=255)
                                                     aaGTCACTAAAATGAGTGCTCAATCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATAGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGTTGGGCGGGCTTCa                                                                                                  >  2:852265/1‑146 (MQ=255)
                                                      aGTCACTAAAATGGGTGCTCAATCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGa                                                                                            <  2:1160079/151‑1 (MQ=255)
                                                         cACTAAAATGGGTGCTCAATCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATc                                                                                          <  1:539195/150‑1 (MQ=255)
                                                           cTAAAATGGGTGCTCAATCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATc                                                                                          >  1:893193/1‑148 (MQ=255)
                                                           cTAAAATGGGTGCTCAATCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCag                                                                                        <  2:640550/150‑1 (MQ=255)
                                                            tAAAATGGGTGCTCAATCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCag                                                                                        >  1:305702/1‑149 (MQ=255)
                                                              aaaTGGGTGCTCAATCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCagag                                                                                      >  2:206992/1‑149 (MQ=255)
                                                                aTGGGTGCTCAATCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTAtt                                                                                  <  1:532455/151‑1 (MQ=255)
                                                                  gggTGCTCAATCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTAttt                                                                                 <  1:903914/150‑1 (MQ=255)
                                                                  gggTGCTCAATCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTATTTc                                                                                <  2:754911/151‑1 (MQ=255)
                                                                   ggTGCTCAATCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTATTTCg                                                                               <  1:497779/151‑1 (MQ=255)
                                                                         cAATCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTATTTCGTTTc                                                                           >  2:837295/1‑149 (MQ=255)
                                                                         cAATCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTATTTCGTTTCAg                                                                         <  1:916310/151‑1 (MQ=255)
                                                                         cAATCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTATTTCGTTTCAg                                                                         <  2:48784/151‑1 (MQ=255)
                                                                          aaTCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTATTTCGTTTCAgg                                                                        <  2:1152772/151‑1 (MQ=255)
                                                                              aTACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTATTTCGTTTCAGGCGAt                                                                    <  2:1239198/151‑1 (MQ=255)
                                                                              aTACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTATTTCGTTTCAGGCGAt                                                                    <  1:896162/151‑1 (MQ=255)
                                                                                aCCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTATTTCGTTTCAGGCGATgc                                                                  >  2:39106/1‑151 (MQ=255)
                                                                                  ccGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTATTTCGTTTCAGGCGATgcg                                                                 >  2:511304/1‑150 (MQ=255)
                                                                                    gTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTATTTCGTTTCAGGCGATGCGCgg                                                              >  1:792523/1‑151 (MQ=255)
                                                                                       tggtggGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTATTTCGTTTCAGGCGATGCGCGGTa                                                            >  1:725840/1‑150 (MQ=255)
                                                                                           gggCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTATTTCGTTTCAGGCGATGCGCGGTAAGtt                                                        <  1:948606/150‑1 (MQ=255)
                                                                                                    gctgctTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGAttt                                               <  1:259502/150‑1 (MQ=255)
                                                                                                    gctgctTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTc                                              >  2:61405/1‑151 (MQ=255)
                                                                                                         ttGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATg                                           >  1:90508/1‑149 (MQ=255)
                                                                                                         ttGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGa                                          >  1:746840/1‑150 (MQ=255)
                                                                                                          tGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGaa                                         >  2:899876/1‑150 (MQ=255)
                                                                                                            gCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAAt                                      <  1:819969/151‑1 (MQ=255)
                                                                                                                cATTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATc                                  <  1:1108358/151‑1 (MQ=255)
                                                                                                                 aTTGCGCTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATcc                                 >  1:846975/1‑151 (MQ=255)
                                                                                                                    gcgcTGAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGt                               <  1:813161/150‑1 (MQ=255)
                                                                                                                     ggctgggtaaggggCTGTTTAGCGGGAACCGGATAACCCTCTTCGCGGATCAGGGGGCTGGCGTGTTGATTAAAAAAAATAAGGCGGATCAGAGTATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTg                              >  2:517414/13‑150 (MQ=255)
                                                                                                                         gAATATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACGAt                         >  2:788088/1‑151 (MQ=255)
                                                                                                                          aaTATCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACGATc                        >  2:92995/1‑151 (MQ=255)
                                                                                                                           atatCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACGATcc                       >  2:406696/1‑151 (MQ=255)
                                                                                                                            tatCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGTACACTGGCGTGTTGATGGGCGGGATTCACGCGGGTCAGAGTATTTCGTTTCAGGCGGTGCGCGGTAAGTTGCAGGCTTTCATGAAGATTATCCGTGACGATCCg                      >  2:1076857/1‑151 (MQ=255)
                                                                                                                            tatCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACGATCCg                      >  2:454120/1‑151 (MQ=255)
                                                                                                                             atCTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACGATCCgg                     >  1:16353/1‑151 (MQ=255)
                                                                                                                               cTGGCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACGATCCGGc                    >  2:585328/1‑150 (MQ=255)
                                                                                                                                 ggCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACGATCCgg                     <  1:1174585/147‑1 (MQ=255)
                                                                                                                                 ggCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACGATCCGGCAGt                 <  2:177761/151‑1 (MQ=255)
                                                                                                                                  gCTGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACGATCCGGCAGTg                <  1:380085/151‑1 (MQ=255)
                                                                                                                                    tGTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACGATCCGGCAGTgg               >  1:15273/1‑150 (MQ=255)
                                                                                                                                     gTATATCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACGATCCGGCAGTGGAt             <  1:885542/151‑1 (MQ=255)
                                                                                                                                        tatCTCGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACGATCCGGCAGTGGATAAt          >  2:675744/1‑151 (MQ=255)
                                                                                                                                            tcGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACGATCCGGCAGTGGATAATGTc       >  2:1108635/1‑150 (MQ=255)
                                                                                                                                             cGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACGATCCGGCAGTGGATAATg         <  1:351367/147‑1 (MQ=255)
                                                                                                                                             cGATCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACGATCCGGCAGTGGATAATGTCAc     <  2:208080/151‑1 (MQ=255)
                                                                                                                                                tCCTGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACGATCCGGCAGTGGATACTGTCACCgg  >  1:1175909/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGCGACTGCGTGGTTTTGGTCGCATGTTGGTAGCCCTGCAACAAGGCTACGGCAAGTCACTAAAATGGGTGCTCAATCATACCCGTCTGGTGGGCGTGGTGCTGCTTGGCACCATTGCGCTGAATATCTGGCTGTATATCTCGATCCCGAAAACCTTCTTCCCGGAGCAGGACACTGGCGTGTTGATGGGCGGGATTCAGGCGGATCAGAGTATTTCGTTTCAGGCGATGCGCGGTAAGTTGCAGGATTTCATGAAAATTATCCGTGACGATCCGGCAGTGGATAATGTCACCGG  >  NC_000913/2159881‑2160175

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