Predicted mutation
evidence seq id position mutation annotation gene description
RA AC_000091 2,919,665 G→A G364G (GGC→GGT)  gudP ← predicted D‑glucarate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AC_0000912,919,6650GA99.2% 451.0 / NA 123G364G (GGC→GGT) gudPpredicted D‑glucarate transporter
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (63/59);  total (64/59)
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.95e-01

Only 100 of 124 total aligned reads displayed.

CGAACATGTTGAACAGGCCACCGGAAAGACCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACA  >  AC_000091/2919580‑2919750
                                                                                     |                                                                                     
cGAACATGTTGAACAGGCCACCGGAAAGACCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCt                                                                                   >  1:1059697/1‑90 (MQ=255)
cGAACATGTTGAACAGGCCACCGGAAAGACCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCt                                                                                   <  2:3271753/90‑1 (MQ=255)
  aaCATGTTGAACAGGCCACCGGAAAGACCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTa                                                                                 <  1:2862700/90‑1 (MQ=255)
     aTGTTGAACAGGCCACCGGAAAGACCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCg                                                                              <  2:916583/90‑1 (MQ=255)
     aTGTTGAACAGGCCACCGGAAAGACCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCg                                                                              >  1:3073747/1‑90 (MQ=255)
     aTGTTGAACAGGCCACCGGAAAGACCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCg                                                                              <  2:2298360/90‑1 (MQ=255)
       gTTGAACAGGCCACCGGAAAGACCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCgaa                                                                            >  2:168122/1‑90 (MQ=255)
         tGAACAGGCCACCGGAAAGACCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCgaaga                                                                          >  1:2819825/1‑90 (MQ=255)
           aaCAGGCCACCGGAAAGACCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGaagaag                                                                        >  1:1028743/1‑90 (MQ=255)
           aaCAGGCCACCGGAAAGACCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGaagaag                                                                        >  1:2342793/1‑90 (MQ=255)
           aaCAGGCCACCGGAAAGACCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGaagaag                                                                        >  1:1267696/1‑90 (MQ=255)
           aaCAGGCCACCGGAAAGACCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGaagaag                                                                        <  2:916735/90‑1 (MQ=255)
             cAGGCCACCGGAAAGACCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGc                                                                      >  2:221741/1‑90 (MQ=255)
             cAGGCCACCGGAAAGACCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGc                                                                      <  1:2216103/90‑1 (MQ=255)
              aGGCCACCGGAAAGACCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGcc                                                                     <  2:3117787/90‑1 (MQ=255)
              aGGCCACCGGAAAGACCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGcc                                                                     >  2:713421/1‑90 (MQ=255)
               ggCCACCGGAAAGACCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCa                                                                    <  1:2482881/90‑1 (MQ=255)
                 ccACCGGAAAGACCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAgc                                                                  >  1:543045/1‑90 (MQ=255)
                 ccACCGGAAAGACCGCTGAGCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAgc                                                                  <  1:2848794/90‑1 (MQ=255)
                  cACCGGAAAGACCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAgcg                                                                 >  2:1432763/1‑90 (MQ=255)
                   aCCGGAAAGACCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAgcgc                                                                >  1:913079/1‑90 (MQ=255)
                     cGGAAAGACCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGCAGAAGGCCAGCGCCa                                                              >  1:2092584/1‑90 (MQ=255)
                     cGGAAAGACCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCa                                                              >  1:2845187/1‑90 (MQ=255)
                      ggAAAGACCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCAt                                                             >  2:2517277/1‑90 (MQ=255)
                      ggAAAGACCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCAt                                                             >  1:260172/1‑90 (MQ=255)
                      ggAAAGACCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCAt                                                             >  2:605074/1‑90 (MQ=255)
                       gAAAGACCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATa                                                            <  2:399776/90‑1 (MQ=255)
                       gAAAGACCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATa                                                            <  2:1089213/90‑1 (MQ=255)
                       gAAAGACCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATa                                                            >  1:2121491/1‑90 (MQ=255)
                        aaaGACCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATaa                                                           <  1:888620/90‑1 (MQ=255)
                        aaaGACCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATaa                                                           >  2:3480561/1‑90 (MQ=255)
                         aaGACCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAGCCCAGCGCCCCGATACCCTTACCAAAGAAGGCCAGCGCCATaaa                                                          >  2:1159082/1‑90 (MQ=255)
                         aaGACCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATaaa                                                          >  1:3096342/1‑90 (MQ=255)
                         aaGACCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATaaa                                                          <  1:798900/90‑1 (MQ=255)
                         aaGACCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATaaa                                                          <  2:13224/90‑1 (MQ=255)
                          aGACCGCTGTTCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAg                                                         <  2:2576886/90‑1 (MQ=255)
                           gACCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGc                                                        <  1:2890602/90‑1 (MQ=255)
                           gACCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGc                                                        <  1:1619493/90‑1 (MQ=255)
                            aCCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGcc                                                       <  2:987216/90‑1 (MQ=255)
                             ccGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACGCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCg                                                      >  2:2044255/1‑90 (MQ=255)
                             ccGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCg                                                      >  1:1807167/1‑90 (MQ=255)
                               gCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCgat                                                    >  2:2860489/1‑90 (MQ=255)
                               gCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCgat                                                    <  2:1394223/90‑1 (MQ=255)
                                 tGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCgatga                                                  <  2:2208542/90‑1 (MQ=255)
                                 tGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCgatga                                                  <  2:160619/90‑1 (MQ=255)
                                 tGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCgatga                                                  <  1:2996147/90‑1 (MQ=255)
                                  taTCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCgatgat                                                 >  2:1989341/2‑90 (MQ=255)
                                  gATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCgatgat                                                 <  1:200822/90‑1 (MQ=255)
                                  gATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCgatgat                                                 >  1:801387/1‑90 (MQ=255)
                                   aTCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGatc                                                >  2:1146379/1‑90 (MQ=255)
                                    tctcTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGatca                                               <  2:1045966/90‑1 (MQ=255)
                                    tctcTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGatca                                               <  2:3173159/90‑1 (MQ=255)
                                    tctcTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGatca                                               <  2:2889245/90‑1 (MQ=255)
                                     ctctTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGatcat                                              <  1:3408493/90‑1 (MQ=255)
                                     ctctTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGatcat                                              <  1:2471676/90‑1 (MQ=255)
                                      tctTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGatcatc                                             >  1:3386844/1‑90 (MQ=255)
                                       ctTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATcc                                            >  2:260217/1‑90 (MQ=255)
                                              cgcgGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAAc                                     <  1:975087/90‑1 (MQ=255)
                                                cgGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGt                                   >  1:233836/1‑90 (MQ=255)
                                                cgGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGt                                   <  1:101047/90‑1 (MQ=255)
                                                cgGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGt                                   >  2:2015788/1‑90 (MQ=255)
                                                 ggTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGtt                                  <  1:2528259/90‑1 (MQ=255)
                                                  gTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTg                                 >  1:1109245/1‑90 (MQ=255)
                                                  gTATCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTg                                 >  2:2338022/1‑90 (MQ=255)
                                                     tCTGCCATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACg                              <  2:1796728/90‑1 (MQ=255)
                                                         ccATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGt                          >  1:2558408/1‑90 (MQ=255)
                                                         ccATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGt                          >  2:3255961/1‑90 (MQ=255)
                                                          cATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGtt                         >  1:1596562/1‑90 (MQ=255)
                                                          cATTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGtt                         <  1:555118/90‑1 (MQ=255)
                                                           aTTACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTg                        <  1:2255238/90‑1 (MQ=255)
                                                            ttACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGc                       <  1:2019674/90‑1 (MQ=255)
                                                             tACTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCa                      >  2:2710628/1‑90 (MQ=255)
                                                              aCTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAg                     <  2:3277966/90‑1 (MQ=255)
                                                               cTGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGa                    >  2:1609600/1‑90 (MQ=255)
                                                                tGCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGaa                   >  1:1166945/1‑90 (MQ=255)
                                                                 gCCCAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAac                  >  2:3256168/1‑90 (MQ=255)
                                                                  cccAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAaca                 <  2:109632/90‑1 (MQ=255)
                                                                   ccAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAacac                <  1:308550/90‑1 (MQ=255)
                                                                   ccAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAacac                >  2:1642378/1‑90 (MQ=255)
                                                                   ccAACCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAacac                <  1:2106452/90‑1 (MQ=255)
                                                                     aaCCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAAcacca              <  1:1831636/90‑1 (MQ=255)
                                                                      aCCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAAcaccat             >  1:2094663/1‑90 (MQ=255)
                                                                      aCCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAAcaccat             <  1:2640011/90‑1 (MQ=255)
                                                                      aCCCAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAAcaccat             >  1:765004/1‑90 (MQ=255)
                                                                        ccAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAAcaccatca           <  2:764328/90‑1 (MQ=255)
                                                                        ccAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAAcaccatca           >  2:653621/1‑90 (MQ=255)
                                                                         cAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAAcaccatcac          >  1:2858071/1‑90 (MQ=255)
                                                                         cAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAAcaccatcac          >  2:2188878/1‑90 (MQ=255)
                                                                         cAGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAAcaccatcac          <  2:2269905/90‑1 (MQ=255)
                                                                          aGCGCCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGGAGTTGCCGAAcaccatcacc         <  2:1446170/90‑1 (MQ=255)
                                                                           gcgcCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAAcaccatcacca        <  2:2048263/90‑1 (MQ=255)
                                                                            cgcCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAAcaccatcaccat       <  1:2862775/90‑1 (MQ=255)
                                                                             gcCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACaccatcaccatc      <  2:3411879/90‑1 (MQ=255)
                                                                             gcCCCGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACaccatcaccatc      <  1:735599/90‑1 (MQ=255)
                                                                              ccccGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCg     <  2:3399879/90‑1 (MQ=255)
                                                                              ccccGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCg     <  2:95683/90‑1 (MQ=255)
                                                                               cccGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGa    >  2:109819/1‑90 (MQ=255)
                                                                               cccGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGa    >  1:1735720/1‑90 (MQ=255)
                                                                               cccGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGa    >  2:3197793/1‑90 (MQ=255)
                                                                               cccGATACCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGa    >  2:3158413/1‑90 (MQ=255)
                                                                                     |                                                                                     
CGAACATGTTGAACAGGCCACCGGAAAGACCGCTGATCTCTTTTGGCGCGGTATCTGCCATTACTGCCCAACCCAGCGCCCCGATGCCCTTACCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACA  >  AC_000091/2919580‑2919750

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: ‑