Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,279,708 A→T 9.9% L18Q (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,7080AT 9.9% 204.1 / 11.4 71L18Q (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 = 4.10e-01

GGCCCGTAGATTAACTGCCCCAACGCAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGACGATATTCACCGGGGGCTAAATCAGCATCCAG  >  NC_000913/2279561‑2279856
                                                                                                                                                   |                                                                                                                                                    
ggCCCGTAGATTAACTGCCCCAACGCAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCcagc                                                                                                                                                    >  1:1394794/1‑150 (MQ=255)
       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:374349/1‑150 (MQ=255)
                                       taAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGt                                                                                                             >  2:1396971/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                                                 <  2:746029/150‑1 (MQ=255)
                                                                                                   ggTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATttct                                                 >  1:495054/1‑150 (MQ=255)
                                                                                                        cgcgGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATttcttcttc                                            <  1:64923/150‑1 (MQ=255)
                                                                                                               aGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAAc                                     <  1:767843/150‑1 (MQ=255)
                                                                                                               aGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAAc                                     <  2:815912/150‑1 (MQ=255)
                                                                                                                gATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACg                                    >  1:653883/1‑150 (MQ=255)
                                                                                                                    cATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAAc                                     >  2:298874/1‑145 (MQ=255)
                                                                                                                    cATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAAc                                     <  1:298874/145‑1 (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)
                                                                                                                                                   |                                                                                                                                                    
GGCCCGTAGATTAACTGCCCCAACGCAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGACGATATTCACCGGGGGCTAAATCAGCATCCAG  >  NC_000913/2279561‑2279856

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: ‑