Predicted mutation
evidence seq id position mutation annotation gene description
RA AC_000091 3,266,304 T→C Y206C (TAC→TGC)  tdcA ← DNA‑binding transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AC_0000913,266,3040TC99.2% 495.0 / NA 132Y206C (TAC→TGC) tdcADNA‑binding transcriptional activator
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (68/63);  total (68/64)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 138 total aligned reads displayed.

TTATAAATTGTCACGACTGAGTCGGTTTTAACGATGTTTTCAATACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGTACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAGCGTGGTGGTGCCGGTGCATGTTCGGGACTTACTG  >  AC_000091/3266217‑3266390
                                                                                       |                                                                                      
ttATAAATTGTCACGACTGAGTCGGTTTTAACGATGTTTTCAATACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCAc                                                                                      <  2:2865921/90‑1 (MQ=255)
 tataAATTGTCACGACTGAGTCGGTTCTAACGATGTTTTCAATACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCAcc                                                                                     <  1:1777650/90‑1 (MQ=255)
  ataAATTGTCACGACTGAGTCGGTTTTAACGATGTTTTCAATACTGATTCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCAccc                                                                                    >  1:1423861/1‑90 (MQ=255)
  ataAATTGTCACGACTGAGTCGGTTTTAACGATGTTTTCAATACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCAccc                                                                                    >  2:1696318/1‑90 (MQ=255)
    aaaTTGTCACGACTGAGTCGGTTTTAACGATGTTTTCAATACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCa                                                                                  >  1:3483479/1‑90 (MQ=255)
    aaaTTGTCACGACTGAGTCGGTTTTAACGATGTTTTCAATACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCa                                                                                  >  2:438965/1‑90 (MQ=255)
    aaaTTGTCACGACTGAGTCGGTTTTAACGATGTTTTCAATACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCa                                                                                  <  1:1834756/90‑1 (MQ=255)
     aaTTGTCACGACTGAGTCGGTTTTAACGATGTTTTCAATACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCat                                                                                 <  2:1419789/90‑1 (MQ=255)
      aTTGTCACGACTGAGTCGGTTTTAACGATGTTTTCAATACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCata                                                                                <  1:2670309/90‑1 (MQ=255)
        tGTCACGACTGAGTCGGTTTTAACGATGTTTTCAATACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATAtt                                                                              <  1:554287/90‑1 (MQ=255)
        tGTCACGACTGAGTCGGTTTTAACGATGTTTTCAATACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATAtt                                                                              <  1:3127286/90‑1 (MQ=255)
          tCACGACTGAGTCGGTTTTAACGATGTTTTCAATACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAg                                                                            >  2:1284938/1‑90 (MQ=255)
           cACGACTGAGTCGGTTTTAACGATGTTTTCAATACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGt                                                                           <  2:1054028/90‑1 (MQ=255)
           cACGACTGAGTCGGTTTTAACGATGTTTTCAATACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGt                                                                           <  1:3056701/90‑1 (MQ=255)
           cACGACTGAGTCGGTTTTAACGATGTTTTCAATACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGt                                                                           <  1:281203/90‑1 (MQ=255)
           cACGACTGAGTCGGTTTTAACGATGTTTTCAATACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGt                                                                           >  2:2688569/1‑90 (MQ=255)
               aCTGAGTCGGTTTTAACGATGTTTTCAATACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTtgt                                                                       <  2:818393/90‑1 (MQ=255)
                cTGAGTCGGTTTTAACGATGTTTTCAATACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTtgtg                                                                      <  2:3456346/90‑1 (MQ=255)
                 tGAGTCGGTTTTAACGATGTTTTCAATACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTgg                                                                     >  2:1355156/1‑90 (MQ=255)
                  gAGTCGGTTTTAACGATGTTTTCAATACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGc                                                                    >  1:1585778/1‑90 (MQ=255)
                   aGTCGGTTTTAACGATGTTTTCAATACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCa                                                                   >  1:3180922/1‑90 (MQ=255)
                   aGTCGGTTTTAACGATGTTTTCAATACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCa                                                                   <  1:3188898/90‑1 (MQ=255)
                    gTCGGTTTTAACGATGTTTTCAATACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCaa                                                                  <  2:3433346/90‑1 (MQ=255)
                     tCGGTTTTAACGATGTTTTCAATACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAac                                                                 <  1:1738901/90‑1 (MQ=255)
                     tCGGTTTTAACGATGTTTTCAATACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAac                                                                 <  2:219315/90‑1 (MQ=255)
                       ggTTTTAACGATGTTTTCAATACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAacac                                                               <  2:2892934/90‑1 (MQ=255)
                       ggTTTTAACGATGTTTTCAATACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAacac                                                               <  2:1013468/90‑1 (MQ=255)
                        gTTTTAACGATGTTTTCAATACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACAcc                                                              >  1:3461288/1‑90 (MQ=255)
                        gTTTTAACGATGTTTTCAATACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGCGGCAAAAcc                                                              >  1:3125208/1‑90 (MQ=255)
                         ttttAACGATGTTTTCAATACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACAccc                                                             >  1:3224666/1‑90 (MQ=255)
                              aCGATGTTTTCAATACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTgt                                                        >  1:1983309/1‑90 (MQ=255)
                                gATGTTTTCAATACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTgttc                                                      <  2:2022322/90‑1 (MQ=255)
                                 aTGTTTTCAATACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTgttcg                                                     >  1:1087093/1‑90 (MQ=255)
                                 aTGTTTTCAATACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTgttcg                                                     <  1:2905622/90‑1 (MQ=255)
                                  tGTTTTCAATACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTgttcgt                                                    >  2:791795/1‑90 (MQ=255)
                                   gTTTTCAATACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTgttcgtt                                                   <  2:50566/90‑1 (MQ=255)
                                    ttttCAATACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACCCCCACTgttcgttc                                                  <  2:1157669/90‑1 (MQ=255)
                                    ttttCAATACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTgttcgttc                                                  >  1:1156689/1‑90 (MQ=255)
                                    ttttCAATACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTgttcgttc                                                  >  1:364306/1‑90 (MQ=255)
                                    ttttCAATACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTgttcgttc                                                  >  1:817185/1‑90 (MQ=255)
                                    ttttCAATACTGATGCCATTTCTTTGTAACGAAGTAAGCAGTTCGCTGTAGCAACCCATATTAGTTTGTGGCAACACCCACTgttcgttc                                                  >  2:1818301/1‑90 (MQ=255)
                                       tCAATACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGttcttc                                               >  2:2129903/1‑90 (MQ=255)
                                       tCAATACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGttcttc                                               >  2:2495440/1‑90 (MQ=255)
                                        cAATACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCa                                              >  2:770951/1‑90 (MQ=255)
                                         aaTACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCaa                                             >  2:522873/1‑90 (MQ=255)
                                         aaTACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCaa                                             >  1:985798/1‑90 (MQ=255)
                                         aaTACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCaa                                             >  2:973511/1‑90 (MQ=255)
                                         aaTACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCaa                                             >  1:1191903/1‑90 (MQ=255)
                                          aTACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAAc                                            <  1:632254/90‑1 (MQ=255)
                                           tACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACg                                           <  1:3207680/90‑1 (MQ=255)
                                           tACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACg                                           <  1:1473345/90‑1 (MQ=255)
                                           tACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACg                                           >  1:185880/1‑90 (MQ=255)
                                            aCTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGa                                          >  1:2694155/1‑90 (MQ=255)
                                              tGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACt                                        >  1:1178054/1‑90 (MQ=255)
                                              tGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACt                                        >  2:2642420/1‑90 (MQ=255)
                                              tGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATAGTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACt                                        >  1:826454/1‑90 (MQ=255)
                                                 tGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCa                                     <  1:2107624/90‑1 (MQ=255)
                                                  gCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAg                                    >  2:1293876/1‑90 (MQ=255)
                                                   ccATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAGc                                   <  2:2397247/90‑1 (MQ=255)
                                                      tttCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAGCgtg                                <  1:2318193/90‑1 (MQ=255)
                                                       ttCTTTGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAGCgtgg                               <  2:2310280/90‑1 (MQ=255)
                                                          tttGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAGCgtggtgg                            <  1:1364226/90‑1 (MQ=255)
                                                           ttGTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAGCgtggtggt                           <  2:1264309/90‑1 (MQ=255)
                                                             gTAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAGCGTGGTGGTGc                         >  2:1634814/1‑90 (MQ=255)
                                                              tAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAGCGTGGTGGTGcc                        <  1:102392/90‑1 (MQ=255)
                                                              tAACGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAGCGTGGTGGTGcc                        >  1:2174402/1‑90 (MQ=255)
                                                               aaCGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAGCGTGGTGGTGCCg                       >  1:2190541/1‑90 (MQ=255)
                                                               aaCGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAGCGTGGTGGTGCCg                       <  2:572610/90‑1 (MQ=255)
                                                                aCGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAGCGTGGTGGTGCCgg                      >  2:967650/1‑90 (MQ=255)
                                                                aCGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAGCGTGGTGGTGCCgg                      <  2:1496363/90‑1 (MQ=255)
                                                                 cGTAGTAAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAGCGTGGTGGTGCCGGt                     >  2:1021407/1‑90 (MQ=255)
                                                                  gtagtaAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAGCGTGGTGGTGCCGGTg                    >  2:2997866/1‑90 (MQ=255)
                                                                   tagtaAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAGCGTGGTGGTGCCGGTGc                   <  1:2636025/90‑1 (MQ=255)
                                                                   tagtaAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAGCGTGGTGGTGCCGGTGc                   >  2:2663781/1‑90 (MQ=255)
                                                                    agtaAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAGCGTGGTGGTGCCGGTGCa                  <  2:3270965/90‑1 (MQ=255)
                                                                    agtaAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAGCGTGGTGGTGCCGGTGCa                  >  1:2957818/1‑90 (MQ=255)
                                                                    agtaAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAGCGTGGTGGTGCCGGTGCa                  <  1:1928111/90‑1 (MQ=255)
                                                                     gtaAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAGCGTGGTGGTGCCGGTGCAt                 <  2:3109314/90‑1 (MQ=255)
                                                                     gtaAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAGCGTGGTGGTGCCGGTGCAt                 <  2:3068959/90‑1 (MQ=255)
                                                                      taAGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAGCGTGGTGGTGCCGGTGCATg                <  1:3280262/90‑1 (MQ=255)
                                                                       aaGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAGCGTGGTGGTGCCGGTGCATGt               <  2:2375322/90‑1 (MQ=255)
                                                                       aaGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAGCGTGGTGGTGCCGGTGCATGt               <  2:636014/90‑1 (MQ=255)
                                                                       aaGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAGCGTGGTGGTGCCGGTGCATGt               >  1:1100963/1‑90 (MQ=255)
                                                                        aGCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAGCGTGGTGGTGCCGGTGCATGtt              >  2:2584699/1‑90 (MQ=255)
                                                                         gCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAGCGTGGTGGTGCCGGTGCATGTTc             <  1:2979811/90‑1 (MQ=255)
                                                                         gCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAGCGTGGTGGTGCCGGTGCATGTTc             <  1:1306886/90‑1 (MQ=255)
                                                                         gCAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAGCGTGGTGGTGCCGGTGCATGTTc             <  2:1518945/90‑1 (MQ=255)
                                                                          cAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAGCGTGGTGGTGCCGGTGCATGTTCg            <  1:1725109/90‑1 (MQ=255)
                                                                          cAGTTCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAGCGTGGTGGTGCCGGTGCATGTTCg            >  2:2590374/1‑90 (MQ=255)
                                                                             ttCGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAGCGTGGTGGTGCCGGTGCATGTTCGGGa         <  1:649612/90‑1 (MQ=255)
                                                                               cGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAGCGTGGTGGTGCCGGTGCATGTTCGGGACt       >  1:1891304/1‑90 (MQ=255)
                                                                               cGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAGCGTGGTGGTGCCGGTGCATGTTCGGGACt       >  2:3188339/1‑90 (MQ=255)
                                                                               cGCTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAGCGTGGTGGTGCCGGTGCATGTTCGGGACt       >  2:406422/1‑90 (MQ=255)
                                                                                 cTGTAGCACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAGCGTGGTGGTGCCGGTGCATGTTCGCGACTTa     <  2:2144817/90‑1 (MQ=255)
                                                                                   gtagcaCCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAGCGTGGTGGTGCCGGTGCATGTTCGGGACTTACt   <  1:2090066/90‑1 (MQ=255)
                                                                                   gtagcaCCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAGCGTGGTGGTGCCGGTGCATGTTCGGGACTTACt   >  1:955764/1‑90 (MQ=255)
                                                                                   gtagcaCCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAGCGTGGTGGTGCCGGTGCATGTTCGGGACTTACt   >  1:2410378/1‑90 (MQ=255)
                                                                                    tagcaCCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAGCGTGGTGGTGCCGGTGCATGTTCGGGACTTACTg  >  2:1349071/1‑90 (MQ=255)
                                                                                    tagcaCCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAGCGTGGTGGTGCCGGTGCATGTTCGGGACTTACTg  <  1:2142414/90‑1 (MQ=255)
                                                                                    tagcaCCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAGCGTGGTGGTGCCGGTGCATGTTCGGGACTTACTg  >  1:2763515/1‑90 (MQ=255)
                                                                                       |                                                                                      
TTATAAATTGTCACGACTGAGTCGGTTTTAACGATGTTTTCAATACTGATGCCATTTCTTTGTAACGTAGTAAGCAGTTCGCTGTAGTACCCCATATTAGTTTGTGGCAACACCCACTGTTCGTTCTTCAACGACTCCAGCGTGGTGGTGCCGGTGCATGTTCGGGACTTACTG  >  AC_000091/3266217‑3266390

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