Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,279,711 A→T 9.9% L17Q (CTG→CAG)  bcr ← multidrug efflux pump Bcr

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,279,7110AT 9.9% 202.6 / 11.4 71L17Q (CTG→CAG) bcrmultidrug efflux pump Bcr
Reads supporting (aligned to +/- strand):  ref base A (35/29);  new base T (4/3);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.47e-01

AGATTAACTGCCCCAACGCAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGACGATATTCACCGGGGGCTAAATCAGCATCCAGCGTA  >  NC_000913/2279568‑2279860
                                                                                                                                               |                                                                                                                                                     
aGATTAACTGCCCCAACGCAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCaa                                                                                                                                                 >  1:950321/1‑150 (MQ=255)
 gATTAACTGCCCCAACGCAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAg                                                                                                                                                <  2:1369466/150‑1 (MQ=255)
  aTTAACTGCCCCAACGCAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCTGCTGGCCatgt                                                                                                                                               <  1:855338/150‑4 (MQ=255)
  aTTAACTGCCCCAACGCAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCTGCTGGCCatgt                                                                                                                                               <  1:849996/150‑4 (MQ=255)
     aaCTGCCCCAACGCAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATa                                                                                                                                            >  1:439707/1‑150 (MQ=255)
     aaCTGCCCCAACGCAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGAGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCGAAGGATa                                                                                                                                            <  2:967492/150‑1 (MQ=255)
         gCCCCAACGCAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATaaaaa                                                                                                                                        >  1:1302370/1‑150 (MQ=255)
             cAACGCAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAAt                                                                                                                                    <  1:1002596/150‑1 (MQ=255)
             cAACGCAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAAt                                                                                                                                    >  2:1325041/1‑150 (MQ=255)
             cAACGCAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAAt                                                                                                                                    <  1:553649/150‑1 (MQ=255)
              aaCGCAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATa                                                                                                                                   <  1:1325041/150‑1 (MQ=255)
              aaCGCAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATa                                                                                                                                   >  1:194305/1‑150 (MQ=255)
                cGCAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGc                                                                                                                                 >  2:912329/1‑150 (MQ=255)
                  cAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGTGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGAt                                                                                                                                             <  2:1059969/136‑1 (MQ=255)
                  cAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGTGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGAt                                                                                                                                             >  1:1059969/1‑136 (MQ=255)
                  cAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCaa                                                                                                                               >  2:289416/1‑150 (MQ=255)
                    aaGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAac                                                                                                                             >  1:988929/1‑150 (MQ=255)
                      gCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGTGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAacga                                                                                                                           >  1:815912/1‑150 (MQ=255)
                      gCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAacga                                                                                                                           >  1:415381/1‑150 (MQ=255)
                      gCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAacga                                                                                                                           >  1:1104162/1‑150 (MQ=255)
                        ccAGAATATAAGTACTGAGGGTCATCTGCGTACCGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAacgacg                                                                                                                         >  2:1019819/1‑150 (MQ=255)
                                taAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGt                                                                                                                 >  2:1396971/1‑150 (MQ=255)
                                taAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGt                                                                                                                 >  2:374349/1‑150 (MQ=255)
                                    tACTGAGTGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCggg                                                                                                             <  1:233302/150‑1 (MQ=255)
                                       tGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGggtgg                                                                                                          >  2:523449/1‑150 (MQ=255)
                                         aGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGTCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTc                                                                                                        <  2:1104162/150‑1 (MQ=255)
                                         aGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTc                                                                                                        >  1:160684/1‑150 (MQ=255)
                                                   gCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCt                                                                                              <  2:1302370/150‑1 (MQ=255)
                                                    cGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTc                                                                                             >  1:1344645/1‑150 (MQ=255)
                                                           cccGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGGTACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAg                                                                                      <  2:1377863/150‑1 (MQ=255)
                                                           cccGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAg                                                                                      <  1:271866/150‑1 (MQ=255)
                                                            ccgccgGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGt                                                                                     >  1:1364911/1‑150 (MQ=255)
                                                              gccgGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCa                                                                                   <  2:73516/150‑1 (MQ=255)
                                                                 ggTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTg                                                                                <  2:299884/150‑1 (MQ=255)
                                                                     cGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGaaaa                                                                            >  1:490280/1‑150 (MQ=255)
                                                                     cGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGaaaa                                                                            >  1:339658/1‑150 (MQ=255)
                                                                      gCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAg                                                                           <  2:232907/150‑1 (MQ=255)
                                                                      gCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAg                                                                           >  1:432798/1‑150 (MQ=255)
                                                                         aaaCTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTa                                                                        <  1:1019819/150‑1 (MQ=255)
                                                                           aCTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAga                                                                      <  1:289416/150‑1 (MQ=255)
                                                                             tGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAgacg                                                                    <  1:989077/150‑1 (MQ=255)
                                                                                 cTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGc                                                                >  2:9476/1‑150 (MQ=255)
                                                                                   gAAATTACCGGTAGCGCGGGCAGATACATATCAATCAACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTg                                                              <  2:586187/150‑1 (MQ=255)
                                                                                      aTTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCa                                                           >  2:228618/1‑150 (MQ=255)
                                                                                         aCCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAAtt                                                        <  1:122176/150‑1 (MQ=255)
                                                                                            ggTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATttct                                                     >  1:495054/1‑150 (MQ=255)
                                                                                            ggTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATttct                                                     <  2:746029/150‑1 (MQ=255)
                                                                                                 cgcgGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATttcttcttc                                                <  1:64923/150‑1 (MQ=255)
                                                                                                        aGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAAc                                         <  2:815912/150‑1 (MQ=255)
                                                                                                        aGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAAc                                         <  1:767843/150‑1 (MQ=255)
                                                                                                         gATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACg                                        >  1:653883/1‑150 (MQ=255)
                                                                                                             cATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAAc                                         <  1:298874/145‑1 (MQ=255)
                                                                                                             cATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAAc                                         >  2:298874/1‑145 (MQ=255)
                                                                                                             cATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGACg                                    <  2:640235/150‑1 (MQ=255)
                                                                                                               tatCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGACGat                                  <  2:1344645/150‑1 (MQ=255)
                                                                                                                    aTCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGACGATATTCa                             <  1:100428/150‑1 (MQ=255)
                                                                                                                      cGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGACGATATTCAcc                           >  1:1237888/1‑150 (MQ=255)
                                                                                                                           gCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGACGATATTCACCggggg                      >  2:1123612/1‑150 (MQ=255)
                                                                                                                           gCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCGGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGACGATATTCACCggggg                      <  2:194305/150‑1 (MQ=255)
                                                                                                                                aTCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGACGATATTCACCGGGGGCTaaa                 >  2:962253/1‑150 (MQ=255)
                                                                                                                                 tCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGACGATATTCACCGGGGGCTAAAt                >  1:519554/1‑150 (MQ=255)
                                                                                                                                  atccttGGCCTGCTGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGACGATATTCACCGGGGGCTAAATc               >  1:639148/6‑150 (MQ=255)
                                                                                                                                  atccttGGCCTGCTGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGACGATATTCACCGGGGGCTAAATc               >  1:514447/6‑150 (MQ=255)
                                                                                                                                   tccttGGCCTGCTGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGACGATATTCACCGGGGGCTAAATCa              >  1:1301877/5‑150 (MQ=255)
                                                                                                                                    ccttGGCCTGCTGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGACGATATTCACCGGGGGCTAAATCAg             <  2:1301877/147‑1 (MQ=255)
                                                                                                                                    aCATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGACGATATTCACCGGGGGCTAAATCAg             <  2:439707/150‑1 (MQ=255)
                                                                                                                                     cttGGCCTGCTGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGACGATATTCACCGGGGGCTAAATCAGc            >  1:994850/3‑150 (MQ=255)
                                                                                                                                        ggCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGACGATATTCACCGGGGGCTAAATCAGCATc         >  2:189549/1‑150 (MQ=255)
                                                                                                                                         gCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGACGATATTCACCGGGGGCTAAATCAGCATcc        >  1:886012/1‑150 (MQ=255)
                                                                                                                                         gCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGGCACAACGGGCGCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCGTCTGCAGGTAACGGACGAGATTCACCGGGGGCTAAATCAGCATcc        >  1:880856/1‑150 (MQ=255)
                                                                                                                                          ccAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGACGATATTCACCGGGGGCTAAATCAGCATCCa       >  1:227590/1‑150 (MQ=255)
                                                                                                                                           cagcagGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGACGATATTCACCGGGGGCTAAATCAGCATCCAg      <  1:9476/150‑1 (MQ=255)
                                                                                                                                               agGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGACGATATTCACCGGGGGCTAAATCAGCATCCAGCGTa  <  1:1123612/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
AGATTAACTGCCCCAACGCAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGACGATATTCACCGGGGGCTAAATCAGCATCCAGCGTA  >  NC_000913/2279568‑2279860

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