Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,333,978 G→A D42N (GAC→AAC)  cysB → N‑acetylserine‑responsive cysteine regulon transcriptional activator; autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,333,9780GA97.8% 438.6 / ‑2.3 136D42N (GAC→AAC) cysBN‑acetylserine‑responsive cysteine regulon transcriptional activator; autorepressor
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (80/53);  total (82/54)
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.84e-01

Only 100 of 136 total aligned reads displayed.

GTTGAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGAAGTCCT  >  NC_000913/1333882‑1334075
                                                                                                |                                                                                                 
gttgAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAAc                                                                                                 >  2:1775613/1‑99 (MQ=255)
    aGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGc                                                                                             >  2:370882/1‑99 (MQ=255)
    aGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGc                                                                                             <  2:1679430/99‑1 (MQ=255)
      gtggtCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTa                                                                                           >  2:3144299/1‑99 (MQ=255)
       tggtCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAg                                                                                          <  1:2459411/99‑1 (MQ=255)
         gtCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGc                                                                                        >  1:2765199/1‑99 (MQ=255)
          tCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCa                                                                                       >  2:1443204/1‑99 (MQ=255)
          tCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCa                                                                                       >  2:270486/1‑99 (MQ=255)
           cAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCAt                                                                                      >  1:704113/1‑99 (MQ=255)
            aaTCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCAtt                                                                                     >  2:220139/1‑99 (MQ=255)
            aaTCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCAtt                                                                                     <  1:16178/99‑1 (MQ=255)
            aaTCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCAtt                                                                                     >  1:1669499/1‑99 (MQ=255)
            aaTCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCAtt                                                                                     >  1:789765/1‑99 (MQ=255)
              tCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGCAAACGAGCTAGGCATTCa                                                                                   >  2:935430/1‑99 (MQ=255)
              tCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCa                                                                                   >  2:1097560/1‑99 (MQ=255)
              tCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCa                                                                                   >  1:3402424/1‑99 (MQ=255)
              tCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCa                                                                                   >  1:3359573/1‑99 (MQ=255)
              tCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCa                                                                                   >  1:3278769/1‑99 (MQ=255)
              tCATAACCTGAATGTCTCATCAACAGCGGAAGCACTTTACACATCACTACCCTCCATCAGTAACCAAATAAGAATACTCTAAAACGAGCTATGCACTCa                                                                                   >  2:3455999/1‑99 (MQ=255)
                aTAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAAAAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCaaa                                                                                 >  1:2911152/1‑99 (MQ=255)
                  aaCCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAAtt                                                                               >  2:257683/1‑99 (MQ=255)
                   aCCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAAttt                                                                              >  2:3882101/1‑99 (MQ=255)
                    ccTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAAtttt                                                                             <  1:1877824/99‑1 (MQ=255)
                    ccTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGAACAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAAtttt                                                                             >  2:1099422/1‑99 (MQ=255)
                       gAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAACCGAGCTAGGCATTCAAATTTTTTc                                                                          <  2:415168/99‑1 (MQ=255)
                       gAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTc                                                                          <  1:93628/99‑1 (MQ=255)
                        aaTGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTcc                                                                         <  2:307251/99‑1 (MQ=255)
                        aaTGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTcc                                                                         <  1:1047908/99‑1 (MQ=255)
                         aTGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTccc                                                                        <  2:2417325/99‑1 (MQ=255)
                         aTGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTccc                                                                        <  2:3270369/99‑1 (MQ=255)
                         aTGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTccc                                                                        <  1:853888/99‑1 (MQ=255)
                            tctcATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGaa                                                                     >  2:1512565/1‑99 (MQ=255)
                            tctcATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGaa                                                                     >  1:2729010/1‑99 (MQ=255)
                              tcatcaACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGc                                                                   >  1:1848312/1‑99 (MQ=255)
                               catcaACAGCGGAAGGACTTTACACATCACAACCCCGGATCAGTAAACAACTCAGAATGCTGGAAAACGAGCTAGGCCTCCAAATTTTTGCCCCaaccc                                                                  >  2:56942/1‑96 (MQ=255)
                                atcaACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCgg                                                                 >  2:2769456/1‑99 (MQ=255)
                                atcaACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCgg                                                                 <  2:2116724/99‑1 (MQ=255)
                                atcaACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCgg                                                                 <  1:3292867/99‑1 (MQ=255)
                                atcaACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCgg                                                                 >  1:2100910/1‑99 (MQ=255)
                                atcaACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCgg                                                                 >  1:3149696/1‑99 (MQ=255)
                                 tcaACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGc                                                                <  1:2889512/99‑1 (MQ=255)
                                   aaCAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCaa                                                              <  2:709998/99‑1 (MQ=255)
                                     cAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGc                                                            >  1:1480464/1‑99 (MQ=255)
                                       gCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCAc                                                          <  1:1643723/99‑1 (MQ=255)
                                       gCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCAc                                                          <  1:1692998/99‑1 (MQ=255)
                                        cGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCAcc                                                         <  2:773832/99‑1 (MQ=255)
                                        cGGAAGGACTTCACACATGCCAACCCGGGATCAGTAAACAACTCAGAATCCTGCAAAACGAGCAAGACACGCAGATTTTTTCCCGAAGCCGGAAGCAcc                                                         >  2:901230/1‑99 (MQ=255)
                                          gAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTg                                                       >  2:253012/1‑99 (MQ=255)
                                           aaGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGTTGCAGCAAGCACCTGa                                                      >  2:2400467/1‑99 (MQ=255)
                                           aaGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGa                                                      >  2:1229077/1‑99 (MQ=255)
                                           aaGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGa                                                      <  1:1537031/99‑1 (MQ=255)
                                             ggACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACg                                                    >  1:4005917/1‑99 (MQ=255)
                                              gACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCACGCACCTGACGc                                                   >  2:3970332/1‑99 (MQ=255)
                                               aCTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCa                                                  <  1:800102/99‑1 (MQ=255)
                                                cTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAg                                                 >  2:3976811/1‑99 (MQ=255)
                                                 tttACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAgg                                                <  2:3661754/99‑1 (MQ=255)
                                                 tttACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAgg                                                <  1:2565671/99‑1 (MQ=255)
                                                    acacATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTaa                                             >  1:3999783/1‑99 (MQ=255)
                                                    acacATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTaa                                             <  2:1401407/99‑1 (MQ=255)
                                                    acacAACACGACCCGGGAGCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGAATACAAATTTTTTTCCGAAGCGGCAAGCACGTGACGCAAGTaa                                             >  2:472821/1‑99 (MQ=255)
                                                       caTCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGc                                          <  1:307055/99‑1 (MQ=255)
                                                        aTCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGcc                                         >  1:178735/1‑99 (MQ=255)
                                                         tCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCca                                        >  1:1041533/1‑99 (MQ=255)
                                                          cacaACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCcag                                       >  2:2449411/1‑99 (MQ=255)
                                                          cacaACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCcag                                       >  1:3458758/1‑99 (MQ=255)
                                                           acaACCCGGGATCAGTAAACAAGTCAGAATCCTGGAAAACGAGCTAGGCATTCAACTTTTTTCCCGAAGCGGCCCCCACCTGACGCAGGCACCGCcagc                                      >  1:3317435/1‑99 (MQ=255)
                                                            tctcccGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCcagca                                     <  1:2608619/96‑1 (MQ=255)
                                                              aCCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAgg                                   >  1:762008/1‑99 (MQ=255)
                                                               cccGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTTTTTTTTTTTCAAGCACCTGACGCAGGTAACGCcaccagtg                                  >  2:3965727/1‑97 (MQ=255)
                                                               cccGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAggg                                  <  2:2239574/99‑1 (MQ=255)
                                                               cccGGGATAAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAggg                                  >  1:3744619/1‑99 (MQ=255)
                                                                ccGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCcagcagaga                                 >  2:2571172/1‑96 (MQ=255)
                                                                ccGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGc                                 >  1:3462621/1‑99 (MQ=255)
                                                                ccGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGc                                 <  2:3011259/99‑1 (MQ=255)
                                                                 cGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCa                                >  2:1524188/1‑99 (MQ=255)
                                                                        aGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAaata                         >  1:2422672/1‑99 (MQ=255)
                                                                        aGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAaata                         >  2:3701288/1‑99 (MQ=255)
                                                                        aGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAaata                         >  1:2350383/1‑99 (MQ=255)
                                                                         gTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAaataa                        <  1:1467755/99‑1 (MQ=255)
                                                                         gTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAaataa                        >  1:4034855/1‑99 (MQ=255)
                                                                          tAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAaataat                       <  1:2954942/99‑1 (MQ=255)
                                                                           aaaCAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAATCTTCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAAtt                      >  2:1368732/1‑99 (MQ=255)
                                                                           aaaCAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAAtt                      >  1:2805404/1‑99 (MQ=255)
                                                                           aaaCAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAAtt                      <  2:2634899/99‑1 (MQ=255)
                                                                            aaCAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTc                     <  2:3568978/99‑1 (MQ=255)
                                                                            aaCAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTc                     <  1:3594490/99‑1 (MQ=255)
                                                                            aaCAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTc                     <  2:2328832/99‑1 (MQ=255)
                                                                            aaCAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGACAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTc                     >  2:1830379/1‑99 (MQ=255)
                                                                               aaGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTa                  >  1:212718/1‑99 (MQ=255)
                                                                                 gTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTAtc                >  1:692796/1‑99 (MQ=255)
                                                                                    aGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTAtcgct             >  2:2707019/1‑99 (MQ=255)
                                                                                       aTGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTAtcgctcgc          <  2:2664621/99‑1 (MQ=255)
                                                                                       aTGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTAtcgctcgc          <  1:1546149/99‑1 (MQ=255)
                                                                                         gCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGa        <  2:1924829/99‑1 (MQ=255)
                                                                                          cTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGaa       >  2:3927659/1‑99 (MQ=255)
                                                                                          cTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGaa       >  1:645641/1‑99 (MQ=255)
                                                                                          cTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGaa       <  1:1315113/99‑1 (MQ=255)
                                                                                          cTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGaa       <  1:1013819/99‑1 (MQ=255)
                                                                                             gAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGAAGTc    >  2:2123211/1‑99 (MQ=255)
                                                                                               aGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCATGCAACTGACGCAGGTAACGCCAGCAGGGCACGACATAATTCGTATGGCTCGCGAAGTCCt  >  2:518911/1‑99 (MQ=255)
                                                                                                |                                                                                                 
GTTGAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGAAGTCCT  >  NC_000913/1333882‑1334075

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