Predicted mutation
evidence seq id position mutation annotation gene description
RA AC_000091 483,561 (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
*AC_000091483,5551.C98.2% 402.4 / 1.1 109coding (73/3150 nt)acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (40/67);  total (41/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 135 total aligned reads displayed.

CGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAG‑‑CCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCT  >  AC_000091/483466‑483639
                                                                                           |                                                                                       
cGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAg                                                                                          >  2:763601/1‑90 (MQ=255)
 gCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAG‑‑c                                                                                       >  2:1563621/1‑90 (MQ=255)
   aTCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAG‑‑ccc                                                                                    >  1:2629457/1‑90 (MQ=255)
   aTCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAG‑‑ccc                                                                                    >  1:2763629/1‑90 (MQ=255)
    tCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAG‑‑cccc                                                                                   >  1:1157150/1‑90 (MQ=255)
    tCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAG‑‑cccc                                                                                   >  1:1931379/1‑90 (MQ=255)
     cAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAG‑‑ccccc                                                                                 >  1:2203975/1‑90 (MQ=255)
     cAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAG‑‑ccccc                                                                                 <  2:1210128/90‑1 (MQ=255)
     cAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAG‑‑ccccc                                                                                 <  1:1820078/90‑1 (MQ=255)
       gcgcCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAG‑‑ccccccc                                                                               <  1:1244278/90‑2 (MQ=255)
        cgcCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAG‑‑ccccccct                                                                              <  1:2822757/90‑3 (MQ=255)
        cgcCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAG‑‑ccccccct                                                                              <  1:2400501/90‑3 (MQ=255)
        cgcCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAG‑‑ccccccct                                                                              <  2:2546283/90‑3 (MQ=255)
        cgcCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAG‑‑ccccccct                                                                              <  1:3056795/90‑3 (MQ=255)
        cgcCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAG‑‑ccccccct                                                                              <  2:1000167/90‑3 (MQ=255)
         gcCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAG‑‑ccccccctg                                                                             >  2:2781476/1‑87 (MQ=255)
         gcCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAG‑‑ccccccctg                                                                             >  2:1884989/1‑87 (MQ=255)
         gcCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAG‑‑ccccccctg                                                                             <  1:3263885/90‑4 (MQ=255)
          ccGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAG‑‑ccccccctgc                                                                            <  1:1034688/90‑5 (MQ=255)
          ccGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAG‑‑ccccccctgc                                                                            >  2:2861329/1‑86 (MQ=255)
            ggggTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCa                                                                           <  2:2234262/90‑1 (MQ=255)
            ggggTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCa                                                                           >  1:688879/1‑90 (MQ=255)
            ggggTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCa                                                                           >  2:2431064/1‑90 (MQ=255)
              ggTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAAc                                                                         <  2:171053/90‑1 (MQ=255)
              ggTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAAc                                                                         <  1:608953/90‑1 (MQ=255)
              ggTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAAc                                                                         >  1:2018044/1‑90 (MQ=255)
                tAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACat                                                                       >  2:1148329/1‑90 (MQ=255)
                tAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACat                                                                       >  2:3304955/1‑90 (MQ=255)
                    aggCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACatgatg                                                                   <  1:859025/90‑1 (MQ=255)
                    aggCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACatgatg                                                                   >  2:1744614/1‑90 (MQ=255)
                        ggAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAG‑‑CCCCGCCTGCCAACATGATGATaa                                                               >  1:2868765/1‑90 (MQ=255)
                          agaTCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATg                                                             <  1:1136537/90‑1 (MQ=255)
                             tCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCg                                                          <  1:2963507/90‑1 (MQ=255)
                             tCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCg                                                          <  1:2821775/90‑1 (MQ=255)
                               gTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGAt                                                        >  1:419556/1‑90 (MQ=255)
                               gTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGAt                                                        <  1:2142822/90‑1 (MQ=255)
                               gTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGAt                                                        <  1:2296737/90‑1 (MQ=255)
                                 tACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCa                                                      >  1:494965/1‑90 (MQ=255)
                                 tACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCa                                                      >  1:2659076/1‑90 (MQ=255)
                                  aCTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAATATGATGATAATGGCGATCAc                                                     <  2:2276761/90‑1 (MQ=255)
                                  aCTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCAc                                                     >  2:2004004/1‑90 (MQ=255)
                                  aCTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCAc                                                     <  2:1207646/90‑1 (MQ=255)
                                  aCTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCAc                                                     <  2:3268260/90‑1 (MQ=255)
                                    tGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCAccc                                                   >  2:2851793/1‑90 (MQ=255)
                                    tGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGAGATCAccc                                                   >  2:1942413/1‑90 (MQ=255)
                                       cggcggTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACg                                                <  1:1335326/90‑1 (MQ=255)
                                        ggcggTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGc                                               >  2:1860740/1‑90 (MQ=255)
                                         gcggTGCAACCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCa                                              <  2:1818528/90‑1 (MQ=255)
                                          cggTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCATGCaa                                             <  1:2649234/90‑1 (MQ=255)
                                            gTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCaaaa                                           <  1:794446/90‑1 (MQ=255)
                                            gTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCaaaa                                           >  1:979702/1‑90 (MQ=255)
                                            gTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCaaaa                                           <  1:2583023/90‑1 (MQ=255)
                                             tGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCaaaa                                           <  2:64905/90‑1 (MQ=255)
                                              gCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAAt                                         <  2:3278228/90‑1 (MQ=255)
                                              gCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAAt                                         <  1:480042/90‑1 (MQ=255)
                                              gCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAAt                                         >  1:2855272/1‑90 (MQ=255)
                                               cAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATc                                        <  2:303479/90‑1 (MQ=255)
                                               cAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATc                                        >  1:3487216/1‑90 (MQ=255)
                                                  tCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCggg                                     <  1:2549459/90‑1 (MQ=255)
                                                  tCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCggg                                     <  2:2644286/90‑1 (MQ=255)
                                                   cGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGc                                    <  2:2338789/90‑1 (MQ=255)
                                                   cGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGc                                    >  1:733552/1‑90 (MQ=255)
                                                    gTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGcg                                   >  2:2394302/1‑90 (MQ=255)
                                                        gATATTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGcgatcg                               <  1:102350/90‑1 (MQ=255)
                                                          tatTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGcgatcgat                             <  1:2264958/90‑1 (MQ=255)
                                                           atTGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGGTCGATa                            >  1:1565563/1‑90 (MQ=255)
                                                            ttGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATaa                           <  2:2896591/90‑1 (MQ=255)
                                                            ttGCGCCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATaa                           <  1:2397588/90‑1 (MQ=255)
                                                             tGCGCCACGGGCAGTTTGAGGATCACCAG‑‑CCTCCCCCGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATaaa                          <  2:965055/90‑1 (MQ=255)
                                                              gcgcCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAg                         <  1:1431689/90‑1 (MQ=255)
                                                              gcgcCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAg                         <  1:358898/90‑1 (MQ=255)
                                                               cgcCACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGa                        >  1:735831/1‑90 (MQ=255)
                                                                 ccACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATGATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGaaa                      <  2:2417021/90‑1 (MQ=255)
                                                                 ccACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGaaa                      <  1:2130276/90‑1 (MQ=255)
                                                                 ccACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGaaa                      <  1:1452354/90‑1 (MQ=255)
                                                                  cACCGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAAt                     <  1:80975/90‑1 (MQ=255)
                                                                    ccGGCAGTTTGAGGATCGCCAGCCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAAtt                    <  2:2544556/90‑1 (MQ=255)
                                                                    ccGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTa                   >  2:247468/1‑90 (MQ=255)
                                                                     cGGCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAg                  >  2:2867676/1‑90 (MQ=255)
                                                                       gCAGTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGc                <  2:2793377/90‑1 (MQ=255)
                                                                         tgTTTGAGGGTCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCCCCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCAt              <  2:2656387/89‑1 (MQ=255)
                                                                          gTTTGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATg             <  2:2768126/90‑1 (MQ=255)
                                                                           tttGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGt            <  1:2039248/90‑1 (MQ=255)
                                                                            ttGAGGATCGCCATC‑‑CCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTc           <  2:376538/90‑1 (MQ=255)
                                                                            ttGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTc           <  1:1647464/90‑1 (MQ=255)
                                                                             tGAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCt          <  2:361337/90‑1 (MQ=255)
                                                                              gAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCtt         >  2:1645252/1‑90 (MQ=255)
                                                                              gAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCtt         >  1:1604669/1‑90 (MQ=255)
                                                                              gAGGATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCtt         <  1:1804503/90‑1 (MQ=255)
                                                                               aGGATCGCCAGCCCCCCCTGCCAACATGCTGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGTCATTTCTTa        <  2:69522/90‑1 (MQ=255)
                                                                                ggATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTaa       <  2:2087678/90‑1 (MQ=255)
                                                                                ggATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTaa       <  1:2585000/90‑1 (MQ=255)
                                                                                 gATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAAc      <  1:1835328/90‑1 (MQ=255)
                                                                                 gATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAAc      >  1:770751/1‑90 (MQ=255)
                                                                                 gATCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAAc      >  2:1461039/1‑90 (MQ=255)
                                                                                  aTCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACg     <  2:3377965/90‑1 (MQ=255)
                                                                                   tCGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACgg    <  1:2418034/90‑1 (MQ=255)
                                                                                    cGCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGc   >  2:923493/1‑90 (MQ=255)
                                                                                     gCCAGCCCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCt  <  1:1459989/90‑1 (MQ=255)
                                                                                     gCCAGCCCCCCCTGCCAACATGATGATAATGGCGACCACCCACGCAAAAATCGAGCGATCGATACAGAAATAAAACATGTCTTACCGGCt  >  2:2884629/1‑90 (MQ=255)
                                                                                           |                                                                                       
CGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCCACCGGCAGTTTGAGGATCGCCAG‑‑CCCCCCTGCCAACATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCT  >  AC_000091/483466‑483639

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: