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

GTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑G‑CCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTGTTCAGGCTGAGCACCGCTTGCGGCTTGCTGGTTATT  >  NC_000913/484182‑484476
                                                                                                                                                       |                                                                                                                                                 
gTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑g                                                                                                                                                    >  2:1196466/1‑150 (MQ=255)
    tCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑G‑ccccc                                                                                                                                              >  1:1181609/1‑151 (MQ=255)
      gATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑G‑ccccccc                                                                                                                                            >  2:623968/1‑150 (MQ=255)
       aTACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑G‑ccccccc                                                                                                                                            >  1:674754/1‑149 (MQ=255)
       aTACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGGTTGAGGATCGCCA‑G‑ccccccc                                                                                                                                            >  1:652438/1‑149 (MQ=255)
            gTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCaa                                                                                                                                       <  1:648393/151‑1 (MQ=255)
             ttCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCaa                                                                                                                                       >  1:60320/1‑150 (MQ=255)
               cATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACat                                                                                                                                    <  2:355015/151‑1 (MQ=255)
                atatTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACatg                                                                                                                                   >  2:413059/1‑151 (MQ=255)
                  atTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACatgat                                                                                                                                 <  2:1123613/151‑1 (MQ=255)
                        ttCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATaa                                                                                                                            >  2:108429/1‑150 (MQ=255)
                         tCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATg                                                                                                                          >  1:1106633/1‑151 (MQ=255)
                          cGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATg                                                                                                                          <  1:102888/150‑1 (MQ=255)
                          cGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATg                                                                                                                          <  2:100990/150‑1 (MQ=255)
                          cGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAACGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGTTGCAATCGTTGGATATTGCGCCACCGGCAGTTTGCGGGTAGCCA‑GCCCCCCCTGCCAACAGGTTGGTAATgg                                                                                                                         >  2:881611/1‑151 (MQ=255)
                           gATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGc                                                                                                                        <  1:1035991/151‑1 (MQ=255)
                           gATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGGGCAATCGTAGGATATTGCGCCACCGGCAGGTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATgg                                                                                                                         >  2:270246/1‑150 (MQ=255)
                              aaCCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGa                                                                                                                      <  1:1016534/150‑1 (MQ=255)
                              aaCCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGAt                                                                                                                     <  1:325084/151‑1 (MQ=255)
                               aCCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGGCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGGGCAATCCTAGGGTATTTCGCCCCCGGCAGGTTTGGGATCGCCA‑GCCCCCCCCGCCCACAAGATGAAAAAgg                                                                                                                         >  1:1018939/1‑146 (MQ=255)
                                     gtCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCa                                                                                                               <  1:256624/150‑1 (MQ=255)
                                        aCCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGc                                                                                                            >  1:1142784/1‑150 (MQ=255)
                                          cGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCaa                                                                                                          >  1:253307/1‑150 (MQ=255)
                                          cGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGGTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCaa                                                                                                          >  1:106972/1‑150 (MQ=255)
                                           gtgtCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCaaaa                                                                                                        <  1:1092881/151‑1 (MQ=255)
                                            tgtCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCaaaa                                                                                                        <  1:623988/150‑1 (MQ=255)
                                            tgtCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCaaaa                                                                                                        >  2:805975/1‑150 (MQ=255)
                                               ccTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATc                                                                                                     <  1:1135638/150‑1 (MQ=255)
                                                cTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCgg                                                                                                   >  2:585631/1‑151 (MQ=255)
                                                    aCTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGGTCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCggg                                                                                                  >  2:1057252/1‑148 (MQ=255)
                                                    aCTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGcga                                                                                               <  2:165774/151‑1 (MQ=255)
                                                      tGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGcgatc                                                                                             >  2:283710/1‑151 (MQ=255)
                                                       gTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGcgatcg                                                                                            <  1:586587/151‑1 (MQ=255)
                                                       gTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGcgatc                                                                                             <  1:505216/150‑1 (MQ=255)
                                                       gTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGcgat                                                                                              <  2:575870/149‑1 (MQ=255)
                                                          ttCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGcgatcgat                                                                                          <  2:932948/150‑1 (MQ=255)
                                                          ttCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGcgatcga                                                                                           <  2:917187/149‑1 (MQ=255)
                                                           tCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATa                                                                                         >  1:769805/1‑150 (MQ=255)
                                                           tCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATa                                                                                         >  1:421205/1‑150 (MQ=255)
                                                           tCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATa                                                                                         <  1:1196452/150‑1 (MQ=255)
                                                             gCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAg                                                                                      <  1:146865/151‑1 (MQ=255)
                                                              cctCAGCGCCGGAGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCTCCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATaa                                                                                        <  2:820580/146‑1 (MQ=255)
                                                              cATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGa                                                                                     <  2:480521/151‑1 (MQ=255)
                                                               aTCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGaa                                                                                    <  1:1215405/151‑1 (MQ=255)
                                                                tCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGaa                                                                                    >  1:241696/1‑150 (MQ=255)
                                                                 cAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAAt                                                                                  <  1:1155723/151‑1 (MQ=255)
                                                                    cgcCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGCTATTGCGCCCCCGGCAGCTTGAGGAACGCCA‑GCCCCCCCTGCCAACACGATGACAACGGCGGTCACCCGCGCCAAAAACGGGCGAACGATAAAGCAAtt                                                                                 >  2:314151/1‑149 (MQ=255)
                                                                       cGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCa                                                                            <  2:1166155/151‑1 (MQ=255)
                                                                       cGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCAAGCAAAAATCGGGCGATCGATAAAGAAATTAgg                                                                              >  2:78637/1‑149 (MQ=255)
                                                                                 ggCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACg                                                                  >  2:1062735/1‑151 (MQ=255)
                                                                                   cGGAGATCGTTACTGCCGGCGGTGCAATTGTAGGATATTGCGCCACCCGCCGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATTAAGATAATGGCGCTCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGc                                                                >  2:705751/1‑151 (MQ=255)
                                                                                       gaTCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCt                                                            <  2:1243675/151‑1 (MQ=255)
                                                                                          cGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGtt                                                         <  1:129491/151‑1 (MQ=255)
                                                                                           gTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGttt                                                        <  2:179906/151‑1 (MQ=255)
                                                                                             tACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTtta                                                       <  2:1055456/150‑1 (MQ=255)
                                                                                              aCTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTtta                                                       >  1:807964/1‑149 (MQ=255)
                                                                                                 gCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTttaagtt                                                   <  1:901447/150‑1 (MQ=255)
                                                                                                        gTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGAGATAACCCAAGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTg                                           >  2:99559/1‑151 (MQ=255)
                                                                                                            aaTCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACt                                       <  2:601369/151‑1 (MQ=255)
                                                                                                             aTCGTAGGATATTGTGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTg                                      >  1:1220935/1‑151 (MQ=255)
                                                                                                                 tAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTGTTc                                   <  2:115702/150‑1 (MQ=255)
                                                                                                                  aGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTGTTc                                   >  1:982774/1‑149 (MQ=255)
                                                                                                                  aGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTGTTCa                                  <  2:620963/150‑1 (MQ=255)
                                                                                                                  aGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTGTTCa                                  <  1:916473/150‑1 (MQ=255)
                                                                                                                        ttGCGCCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTGTTCAGGCTg                             >  1:719680/1‑149 (MQ=255)
                                                                                                                          gcgcCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTGTTCAGGCTGAGc                          >  1:938165/1‑150 (MQ=255)
                                                                                                                           cgcCACCGGCAGTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTGTTCAGGCTGAGCAc                        >  1:641878/1‑151 (MQ=255)
                                                                                                                               aCCGGCAGTTTGAGGATCGCCAGG‑CCCCCCTGCCAACATGATGATAATGGAGATCACCCACGCAAAAATCGGGGGAACGATAAAGCAGTTAGGCCTGTCTTAACGGCTACTGTTTAAGTTAAGGCTTGGACTGTTCAGGCTGATCAc                        >  2:87517/1‑147 (MQ=255)
                                                                                                                                      gTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTGTTCAGGCTGAGCACCGCTTGCGGCt             >  2:172805/1‑151 (MQ=255)
                                                                                                                                      gTTTGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTGTTCAGGCTGAGCACCGCTTGCGGCt             <  2:251941/151‑1 (MQ=255)
                                                                                                                                         tGAGGATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTGTTCAGGCTGAGCACCGCTTGCGGCTTg           <  2:167705/150‑1 (MQ=255)
                                                                                                                                             gATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTGTTCAGGCTGAGCACCGCTTGGGGCTTGCTGGt      <  1:109571/151‑1 (MQ=255)
                                                                                                                                             gATCGCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTGTTCAGGCTGAGCACCGCTTGCGGCTTGCTGGt      >  1:23598/1‑151 (MQ=255)
                                                                                                                                                 gCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTGTTCAGGCTGAGCACCGCTTGCGGCTTGCTGGttatt  <  1:413071/151‑1 (MQ=255)
                                                                                                                                                 gCCA‑GCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTGTTCAGGCTGAGCACCGCTTGCGGCTTGCTGGttatt  >  1:187307/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                 
GTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCA‑G‑CCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTGTTCAGGCTGAGCACCGCTTGCGGCTTGCTGGTTATT  >  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: ‑