Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 484,337 (C)6→7 coding (67/3150 nt) acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913484,3311.C98.6% 237.7 / ‑2.9 70coding (73/3150 nt)acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (31/38);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

GTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTGTTCAGGCTGAGCACCGCTTGCGGCTTGCTGGTTATT  >  NC_000913/484182‑484476
                                                                                                                                                       |                                                                                                                                                 
gTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAg                                                                                                                                                    >  2:1196466/1‑150 (MQ=255)
    tCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGccccc                                                                                                                                              >  1:1181609/1‑151 (MQ=255)
      gATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGccccccc                                                                                                                                            >  2:623968/1‑150 (MQ=255)
       aTACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGccccccc                                                                                                                                            >  1:674754/1‑149 (MQ=255)
       aTACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGGTTGAGGATCGCCAGccccccc                                                                                                                                            >  1:652438/1‑149 (MQ=255)
            gTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCaa                                                                                                                                       <  1:648393/151‑1 (MQ=255)
             ttCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCaa                                                                                                                                       >  1:60320/1‑150 (MQ=255)
               cATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACat                                                                                                                                    <  2:355015/151‑1 (MQ=255)
                atatTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACatg                                                                                                                                   >  2:413059/1‑151 (MQ=255)
                  atTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACatgat                                                                                                                                 <  2:1123613/151‑1 (MQ=255)
                        ttCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATaa                                                                                                                            >  2:108429/1‑150 (MQ=255)
                         tCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATg                                                                                                                          >  1:1106633/1‑151 (MQ=255)
                          cGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATg                                                                                                                          <  1:102888/150‑1 (MQ=255)
                          cGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATg                                                                                                                          <  2:100990/150‑1 (MQ=255)
                          cGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAACGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGTTGCAATCGTTGGATATTGCGCCACCGGCAGTTTGCGGGTAGCCAGCCCCCCCTGCCAACAGGTTGGTAATgg                                                                                                                         >  2:881611/1‑151 (MQ=255)
                           gATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGc                                                                                                                        <  1:1035991/151‑1 (MQ=255)
                           gATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGGGCAATCGTAGGATATTGCGCCACCGGCAGGTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATgg                                                                                                                         >  2:270246/1‑150 (MQ=255)
                              aaCCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGa                                                                                                                      <  1:1016534/150‑1 (MQ=255)
                              aaCCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGAt                                                                                                                     <  1:325084/151‑1 (MQ=255)
                               aCCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGGCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGGGCAATCCTAGGGTATTTCGCCCCCGGCAGGTTTGGGATCGCCAGCCCCCCCCGCCCACAAGATGAAAAAgg                                                                                                                         >  1:1018939/1‑146 (MQ=255)
                                     gtCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCa                                                                                                               <  1:256624/150‑1 (MQ=255)
                                        aCCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGc                                                                                                            >  1:1142784/1‑150 (MQ=255)
                                          cGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCaa                                                                                                          >  1:253307/1‑150 (MQ=255)
                                          cGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGGTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCaa                                                                                                          >  1:106972/1‑150 (MQ=255)
                                           gtgtCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCaaaa                                                                                                        <  1:1092881/151‑1 (MQ=255)
                                            tgtCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCaaaa                                                                                                        <  1:623988/150‑1 (MQ=255)
                                            tgtCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCaaaa                                                                                                        >  2:805975/1‑150 (MQ=255)
                                               ccTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATc                                                                                                     <  1:1135638/150‑1 (MQ=255)
                                                cTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCgg                                                                                                   >  2:585631/1‑151 (MQ=255)
                                                    aCTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGGTCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCggg                                                                                                  >  2:1057252/1‑148 (MQ=255)
                                                    aCTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGcga                                                                                               <  2:165774/151‑1 (MQ=255)
                                                      tGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGcgatc                                                                                             >  2:283710/1‑151 (MQ=255)
                                                       gTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGcgatcg                                                                                            <  1:586587/151‑1 (MQ=255)
                                                       gTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGcgatc                                                                                             <  1:505216/150‑1 (MQ=255)
                                                       gTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGcgat                                                                                              <  2:575870/149‑1 (MQ=255)
                                                          ttCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGcgatcgat                                                                                          <  2:932948/150‑1 (MQ=255)
                                                          ttCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGcgatcga                                                                                           <  2:917187/149‑1 (MQ=255)
                                                           tCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATa                                                                                         >  1:769805/1‑150 (MQ=255)
                                                           tCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATa                                                                                         >  1:421205/1‑150 (MQ=255)
                                                           tCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATa                                                                                         <  1:1196452/150‑1 (MQ=255)
                                                             gCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAg                                                                                      <  1:146865/151‑1 (MQ=255)
                                                              cctCAGCGCCGGAGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCTCCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATaa                                                                                        <  2:820580/146‑1 (MQ=255)
                                                              cATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGa                                                                                     <  2:480521/151‑1 (MQ=255)
                                                               aTCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGaa                                                                                    <  1:1215405/151‑1 (MQ=255)
                                                                tCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGaa                                                                                    >  1:241696/1‑150 (MQ=255)
                                                                 cAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAAt                                                                                  <  1:1155723/151‑1 (MQ=255)
                                                                    cgcCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGCTATTGCGCCCCCGGCAGCTTGAGGAACGCCAGCCCCCCCTGCCAACACGATGACAACGGCGGTCACCCGCGCCAAAAACGGGCGAACGATAAAGCAAtt                                                                                 >  2:314151/1‑149 (MQ=255)
                                                                       cGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCa                                                                            <  2:1166155/151‑1 (MQ=255)
                                                                       cGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCAAGCAAAAATCGGGCGATCGATAAAGAAATTAgg                                                                              >  2:78637/1‑149 (MQ=255)
                                                                                 ggCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACg                                                                  >  2:1062735/1‑151 (MQ=255)
                                                                                   cGGAGATCGTTACTGCCGGCGGTGCAATTGTAGGATATTGCGCCACCCGCCGTTTGAGGATCGCCAGCCCCCCCTGCCAACATTAAGATAATGGCGCTCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGc                                                                >  2:705751/1‑151 (MQ=255)
                                                                                       gaTCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCt                                                            <  2:1243675/151‑1 (MQ=255)
                                                                                          cGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGtt                                                         <  1:129491/151‑1 (MQ=255)
                                                                                           gTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGttt                                                        <  2:179906/151‑1 (MQ=255)
                                                                                             tACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTtta                                                       <  2:1055456/150‑1 (MQ=255)
                                                                                              aCTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTtta                                                       >  1:807964/1‑149 (MQ=255)
                                                                                                 gCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTttaagtt                                                   <  1:901447/150‑1 (MQ=255)
                                                                                                        gTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGAGATAACCCAAGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTg                                           >  2:99559/1‑151 (MQ=255)
                                                                                                            aaTCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACt                                       <  2:601369/151‑1 (MQ=255)
                                                                                                             aTCGTAGGATATTGTGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTg                                      >  1:1220935/1‑151 (MQ=255)
                                                                                                                 tAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTGTTc                                   <  2:115702/150‑1 (MQ=255)
                                                                                                                  aGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTGTTc                                   >  1:982774/1‑149 (MQ=255)
                                                                                                                  aGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTGTTCa                                  <  2:620963/150‑1 (MQ=255)
                                                                                                                  aGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTGTTCa                                  <  1:916473/150‑1 (MQ=255)
                                                                                                                        ttGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTGTTCAGGCTg                             >  1:719680/1‑149 (MQ=255)
                                                                                                                          gcgcCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTGTTCAGGCTGAGc                          >  1:938165/1‑150 (MQ=255)
                                                                                                                           cgcCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTGTTCAGGCTGAGCAc                        >  1:641878/1‑151 (MQ=255)
                                                                                                                               aCCGGCAGTTTGAGGATCGCCAGGCCCCCCTGCCAACATGATGATAATGGAGATCACCCACGCAAAAATCGGGGGAACGATAAAGCAGTTAGGCCTGTCTTAACGGCTACTGTTTAAGTTAAGGCTTGGACTGTTCAGGCTGATCAc                        >  2:87517/1‑147 (MQ=255)
                                                                                                                                      gTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTGTTCAGGCTGAGCACCGCTTGCGGCt             >  2:172805/1‑151 (MQ=255)
                                                                                                                                      gTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTGTTCAGGCTGAGCACCGCTTGCGGCt             <  2:251941/151‑1 (MQ=255)
                                                                                                                                         tGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTGTTCAGGCTGAGCACCGCTTGCGGCTTg           <  2:167705/150‑1 (MQ=255)
                                                                                                                                             gATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTGTTCAGGCTGAGCACCGCTTGGGGCTTGCTGGt      <  1:109571/151‑1 (MQ=255)
                                                                                                                                             gATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTGTTCAGGCTGAGCACCGCTTGCGGCTTGCTGGt      >  1:23598/1‑151 (MQ=255)
                                                                                                                                                 gCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTGTTCAGGCTGAGCACCGCTTGCGGCTTGCTGGttatt  <  1:413071/151‑1 (MQ=255)
                                                                                                                                                 gCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTGTTCAGGCTGAGCACCGCTTGCGGCTTGCTGGttatt  >  1:187307/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                 
GTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTGTTCAGGCTGAGCACCGCTTGCGGCTTGCTGGTTATT  >  NC_000913/484182‑484476

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