Predicted mutation
evidence seq id position mutation annotation gene description
RA AC_000091 51,338 +G coding (269/378 nt) apaG ← protein associated with Co2+ and Mg2+ efflux

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AC_00009151,3341.G96.0% 379.7 / 10.0 101coding (273/378 nt)apaGprotein associated with Co2+ and Mg2+ efflux
Reads supporting (aligned to +/- strand):  ref base . (3/1);  new base G (50/47);  total (53/48)
Fisher's exact test for biased strand distribution p-value = 6.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.74e-01

Only 100 of 111 total aligned reads displayed.

GAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGC  >  AC_000091/51246‑51416
                                                                                         |                                                                                  
gAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCg                                                                                   <  2:361190/90‑1 (MQ=255)
 aaCGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCgg                                                                                  >  1:1787705/1‑90 (MQ=255)
  aCGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCggg                                                                                 <  1:1491331/90‑1 (MQ=255)
   cGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCgggg                                                                                >  1:2166327/1‑90 (MQ=255)
    ggCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCggggg                                                                               <  1:3249569/90‑2 (MQ=255)
     gCGAGTCGGAATACGGGAATGTCGATGCTGAAATGTACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCgggggt                                                                              <  2:1713613/90‑3 (MQ=255)
      cGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGgtt                                                                             >  2:220302/1‑90 (MQ=255)
       gAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGgttt                                                                            >  1:2405343/1‑89 (MQ=255)
        aGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGgtttc                                                                           >  2:930137/1‑88 (MQ=255)
        aGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGgtttc                                                                           <  1:1330641/90‑3 (MQ=255)
           cGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAAt                                                                         >  2:1773805/1‑90 (MQ=255)
             gAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGa                                                                       <  1:1191533/90‑1 (MQ=255)
              aaTACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGAt                                                                      >  1:443867/1‑90 (MQ=255)
               aTACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGAtt                                                                     <  2:2396380/90‑1 (MQ=255)
                tACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTg                                                                    <  2:1981802/90‑1 (MQ=255)
                 aCGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGc                                                                   <  1:1957836/90‑1 (MQ=255)
                  cGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCa                                                                  >  2:1110858/1‑90 (MQ=255)
                  cGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCa                                                                  <  2:2988928/90‑1 (MQ=255)
                   gggAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCAc                                                                 <  2:266233/90‑1 (MQ=255)
                    ggAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCAcc                                                                >  2:3442978/1‑90 (MQ=255)
                    ggAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCAcc                                                                >  1:1414348/1‑90 (MQ=255)
                     gAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCg                                                               <  1:1177955/90‑1 (MQ=255)
                       aTGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCt                                                             <  1:1977619/90‑1 (MQ=255)
                       aTGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCt                                                             <  2:3095848/90‑1 (MQ=255)
                       aTGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCt                                                             <  1:2567720/90‑1 (MQ=255)
                         gTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTgg                                                           <  1:1202239/90‑1 (MQ=255)
                          tCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGgt                                                          <  1:1937431/90‑1 (MQ=255)
                          tCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGgt                                                          >  1:2661882/1‑90 (MQ=255)
                          tCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGgt                                                          >  1:26146/1‑90 (MQ=255)
                           cGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGgtg                                                         <  2:2722130/90‑1 (MQ=255)
                           cGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGgtg                                                         >  1:3016191/1‑90 (MQ=255)
                           cGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGgtg                                                         >  1:1657664/1‑90 (MQ=255)
                              tGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTAc                                                      >  1:3456155/1‑90 (MQ=255)
                               gCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACt                                                     <  2:1994574/90‑1 (MQ=255)
                               gCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACt                                                     <  1:1566203/90‑1 (MQ=255)
                               gCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACt                                                     <  1:2388482/90‑1 (MQ=255)
                               gCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACt                                                     <  1:223469/90‑1 (MQ=255)
                                cTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTg                                                    <  2:668616/90‑1 (MQ=255)
                                   aaaGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCGCAGCGGGGGATACAATGATTGCACCGCTGGGGTACTGGTa                                                 >  2:1392832/1‑90 (MQ=255)
                                    aaGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTAc                                                >  2:1658526/1‑90 (MQ=255)
                                    aaGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTAc                                                >  2:1189507/1‑90 (MQ=255)
                                     aGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTAct                                               <  2:2551809/90‑1 (MQ=255)
                                        gACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCtt                                            >  2:2316473/1‑90 (MQ=255)
                                        gACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCtt                                            >  1:1981053/1‑90 (MQ=255)
                                        gACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCtt                                            <  1:1935527/90‑1 (MQ=255)
                                         aCGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTc                                           >  2:887136/1‑90 (MQ=255)
                                         aCGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTc                                           <  2:114780/90‑1 (MQ=255)
                                           gCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGc                                         <  1:2894371/90‑1 (MQ=255)
                                           gCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGc                                         <  2:3488266/90‑1 (MQ=255)
                                            ccATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGcc                                        <  1:3480041/90‑1 (MQ=255)
                                             cATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCa                                       <  2:3361194/90‑1 (MQ=255)
                                              aTTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAg                                      >  1:2707589/1‑90 (MQ=255)
                                               ttttCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAgg                                     <  2:2193495/90‑1 (MQ=255)
                                               ttttCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGAACTCTTCGCCAgg                                     >  1:3136093/1‑90 (MQ=255)
                                                tttCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGgc                                    >  2:470731/1‑90 (MQ=255)
                                                 ttCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGgcg                                   >  1:1842182/1‑90 (MQ=255)
                                                 ttCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGgcg                                   <  2:485762/90‑1 (MQ=255)
                                                 ttCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGgcg                                   <  2:519858/90‑1 (MQ=255)
                                                  tCATCTATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCAGGTGTACTGGTACTCTTCGCCAGgcgc                                  >  2:2840743/1‑90 (MQ=255)
                                                  tCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGgcgc                                  >  2:866754/1‑90 (MQ=255)
                                                   cATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGcgcg                                 >  2:1273661/1‑90 (MQ=255)
                                                   cATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGcgcg                                 >  2:2359926/1‑90 (MQ=255)
                                                    aTCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGa                                >  2:2186845/1‑90 (MQ=255)
                                                     tCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGAt                               >  1:2867909/1‑90 (MQ=255)
                                                      cGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTGAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATa                              >  2:1685929/1‑90 (MQ=255)
                                                      cGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATa                              >  1:3195922/1‑90 (MQ=255)
                                                      cGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATa                              <  2:2539537/90‑1 (MQ=255)
                                                       gATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATaa                             >  1:548466/1‑90 (MQ=255)
                                                       gATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATaa                             >  1:2415386/1‑90 (MQ=255)
                                                        aTCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAg                            >  2:3299080/1‑90 (MQ=255)
                                                        aTCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAg                            >  2:2267369/1‑90 (MQ=255)
                                                        aTCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAg                            >  2:1514301/1‑90 (MQ=255)
                                                        aTCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGGTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAg                            >  1:2941311/1‑90 (MQ=255)
                                                         tCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGt                           >  1:2894726/1‑90 (MQ=255)
                                                          cATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTg                          >  1:1090859/1‑90 (MQ=255)
                                                          cATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTg                          >  2:2048109/1‑90 (MQ=255)
                                                            tttCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGc                        >  2:139660/1‑90 (MQ=255)
                                                            tttCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGc                        >  2:1636273/1‑90 (MQ=255)
                                                             ttCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCt                       <  2:368168/90‑1 (MQ=255)
                                                              tCGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTg                      <  2:560459/90‑1 (MQ=255)
                                                               cGTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTgg                     <  1:462930/90‑1 (MQ=255)
                                                                gTAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGa                    >  1:1505221/1‑90 (MQ=255)
                                                                 tAGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGAc                   <  1:2527506/90‑1 (MQ=255)
                                                                  aGTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACg                  >  1:3473422/1‑90 (MQ=255)
                                                                   gTGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGc                 >  2:3480368/1‑90 (MQ=255)
                                                                    tGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGcc                <  2:729428/90‑1 (MQ=255)
                                                                    tGACCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGcc                <  2:3373872/90‑1 (MQ=255)
                                                                      aCCCTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCaa              >  1:399206/1‑90 (MQ=255)
                                                                       cccTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAAc             >  2:2137607/1‑90 (MQ=255)
                                                                        ccTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAAcc            >  1:738072/1‑90 (MQ=255)
                                                                         cTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCa           >  2:1873127/1‑90 (MQ=255)
                                                                         cTGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCAAGGCGCGATAAGTGGCTGGACGCCAACCa           <  1:2632201/90‑1 (MQ=255)
                                                                          tGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCAc          <  2:3271787/90‑1 (MQ=255)
                                                                          tGCATGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCAc          >  2:2966411/1‑90 (MQ=255)
                                                                              tGGTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCt      <  2:2413528/90‑1 (MQ=255)
                                                                               ggTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCtt     <  1:3128023/90‑1 (MQ=255)
                                                                                gTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTc    <  1:8030/90‑1 (MQ=255)
                                                                                gTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTc    >  2:3086674/1‑90 (MQ=255)
                                                                                gTGCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTc    <  1:2349286/90‑1 (MQ=255)
                                                                                  gCCCAGCGGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGc  <  1:1004835/90‑1 (MQ=255)
                                                                                         |                                                                                  
GAACGGCGAGTCGGAATACGGGAATGTCGATGCTGAAAGGGACGCCATTTTCATCGATCATTTCGTAGTGACCCTGCATGGTGCCCAGCGGGGTTTCAATGATTGCACCGCTGGTGTACTGGTACTCTTCGCCAGGCGCGATAAGTGGCTGGACGCCAACCACTCCTTCGC  >  AC_000091/51246‑51416

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