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,9780GA96.7% 660.6 / 3.9 213D42N (GAC→AAC) cysBN‑acetylserine‑responsive cysteine regulon transcriptional activator; autorepressor
Reads supporting (aligned to +/- strand):  ref base G (3/4);  new base A (98/108);  total (101/112)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 213 total aligned reads displayed.

GTTGAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGAAGTC  >  NC_000913/1333882‑1334073
                                                                                                |                                                                                               
gttgAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAAc                                                                                               <  1:4869666/99‑1 (MQ=255)
gttgAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAAc                                                                                               <  2:1452190/99‑1 (MQ=255)
gttgAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAAc                                                                                               <  2:1617525/99‑1 (MQ=255)
  tgAGGTGTTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGa                                                                                             <  2:818047/99‑1 (MQ=255)
  tgAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGa                                                                                             <  2:91323/99‑1 (MQ=255)
  tgAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGa                                                                                             <  2:622481/99‑1 (MQ=255)
    aGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACAGCACAACCCAGGAACAGAAAACAAGTCAGAATGCTGGAAAACGAGc                                                                                           >  1:4467318/1‑99 (MQ=255)
      gtggtCAGTCATAACCTGAATGTCTCATCAACAACGGAAGAGCTATACACATCACAACGCGGTACGAGTAAACAAGTCATAATGATGGAAAACGAGCTa                                                                                         >  2:1825411/1‑99 (MQ=255)
      gtggtCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTa                                                                                         <  2:4501734/99‑1 (MQ=255)
      gtggtCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTa                                                                                         >  1:3232056/1‑99 (MQ=255)
      gtggtCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTa                                                                                         >  1:977618/1‑99 (MQ=255)
         gtCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGc                                                                                      >  2:4620536/1‑99 (MQ=255)
            aataatACCCTGAATATATCATCAAAAGAGGAAGGACCTTACACACCACCACACGGGATCAGTAAACAATGCAGAATGCTGGAAGCCGAGCTAGGAAtt                                                                                   >  2:4398738/1‑99 (MQ=255)
            aaTCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATAAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCAtt                                                                                   >  1:4291538/1‑99 (MQ=255)
             aTCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTc                                                                                  >  1:542122/1‑99 (MQ=255)
              tCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCa                                                                                 >  1:3984629/1‑99 (MQ=255)
                aTAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAATCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCaaa                                                                               <  1:691085/99‑1 (MQ=255)
                aTAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCaaa                                                                               <  1:903584/99‑1 (MQ=255)
                 tAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAAt                                                                              <  2:530760/99‑1 (MQ=255)
                 tAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAAt                                                                              <  1:4084980/99‑1 (MQ=255)
                 tAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAAt                                                                              <  1:3077714/99‑1 (MQ=255)
                    ccTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAAtttt                                                                           <  2:3271291/99‑1 (MQ=255)
                    ccTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAAtttt                                                                           <  1:3355618/99‑1 (MQ=255)
                    ccTGAATGTCTCATCAACACCGGAAGGACTTTAAACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAAtttt                                                                           <  2:2434060/99‑1 (MQ=255)
                     cTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAAttttt                                                                          <  2:3925795/99‑1 (MQ=255)
                       gAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTc                                                                        <  1:186119/99‑1 (MQ=255)
                       gAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTc                                                                        <  1:1845591/99‑1 (MQ=255)
                         aTGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTccc                                                                      <  1:5058826/99‑1 (MQ=255)
                          tGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCg                                                                     >  1:3592650/1‑99 (MQ=255)
                              tcatcaACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGc                                                                 <  2:138722/99‑1 (MQ=255)
                              tcatcaACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGc                                                                 <  1:1720557/99‑1 (MQ=255)
                               catcaACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCg                                                                <  1:1083776/99‑1 (MQ=255)
                               catcaACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCg                                                                >  1:5218676/1‑99 (MQ=255)
                               catcaACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCg                                                                >  2:3698324/1‑99 (MQ=255)
                                atcaACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCgg                                                               >  2:83507/1‑99 (MQ=255)
                                atcaACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCgg                                                               >  2:2243853/1‑99 (MQ=255)
                                atcaACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCgg                                                               >  1:255740/1‑99 (MQ=255)
                                  caACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCa                                                             >  1:1242440/1‑99 (MQ=255)
                                   aaCAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCaa                                                            <  1:2416835/99‑1 (MQ=255)
                                    aCAGCGGAAGGACTTTACACATCACAACCCGGGATCATTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAg                                                           <  1:104734/99‑1 (MQ=255)
                                    aCAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAg                                                           <  1:4298192/99‑1 (MQ=255)
                                    aCAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAg                                                           <  1:2892569/99‑1 (MQ=255)
                                     cAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGTAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGc                                                          <  1:894190/99‑1 (MQ=255)
                                     cAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGc                                                          <  2:3021069/99‑1 (MQ=255)
                                      aGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCa                                                         >  2:4233884/1‑99 (MQ=255)
                                       gCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCAc                                                        <  1:3483569/99‑1 (MQ=255)
                                       gCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCAc                                                        <  1:154160/99‑1 (MQ=255)
                                        cGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCAcc                                                       <  1:3467444/99‑1 (MQ=255)
                                          gAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTg                                                     >  2:2749274/1‑99 (MQ=255)
                                          gAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTCCCCGAAGCGGCAAGCACCta                                                     >  1:5019417/1‑98 (MQ=255)
                                            aGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGAc                                                   >  1:623048/1‑99 (MQ=255)
                                            aGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGAc                                                   >  1:4684779/1‑99 (MQ=255)
                                            aGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACAAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGAc                                                   <  1:3022174/99‑1 (MQ=255)
                                             ggACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACg                                                  <  2:1776733/99‑1 (MQ=255)
                                              gACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGc                                                 <  1:1280359/99‑1 (MQ=255)
                                              gACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGc                                                 <  2:65577/99‑1 (MQ=255)
                                                cTTTACACATCACAACCCGGGAACAGTAATCAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAg                                               <  2:4559447/99‑1 (MQ=255)
                                                 tttACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAgg                                              >  2:4656522/1‑99 (MQ=255)
                                                  ttACACATCACAACCCGGGAGCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGGCGCGGCAAGCACCTGACGCAGGt                                             <  2:3122826/99‑1 (MQ=255)
                                                    acacATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTaa                                           <  1:1476856/99‑1 (MQ=255)
                                                    acacATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTaa                                           <  1:4997232/99‑1 (MQ=255)
                                                     cacaTCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAAc                                          <  2:34356/99‑1 (MQ=255)
                                                     cacaTCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAAc                                          <  2:1160447/99‑1 (MQ=255)
                                                      acaTCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACg                                         >  1:111091/1‑99 (MQ=255)
                                                      acaTCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACg                                         >  1:3250825/1‑99 (MQ=255)
                                                      aaaacacaACCCGGGATCAGTAACAAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGCTAACg                                         <  2:1033247/95‑1 (MQ=255)
                                                        aTCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGcc                                       >  2:129178/1‑99 (MQ=255)
                                                        aTCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGcc                                       >  1:1639896/1‑99 (MQ=255)
                                                        aTCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGcc                                       <  2:2957572/99‑1 (MQ=255)
                                                        aTCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGcc                                       >  2:2675060/1‑99 (MQ=255)
                                                        aTCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAAGTTTTTCCCGAAGCGGCAAGCTCCTGACGCAGGTAACGcc                                       >  1:4939349/1‑99 (MQ=255)
                                                         tCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCAAGAAGCGGCAAGCACCTGACGCAGGTAACGCca                                      >  1:3837741/1‑99 (MQ=255)
                                                          cacaACCCGGGATCAGTAAACAAGTCATAATCCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCcag                                     >  2:4004029/1‑99 (MQ=255)
                                                          cacaACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCcag                                     <  2:494239/99‑1 (MQ=255)
                                                               cccGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAggg                                <  2:3180789/99‑1 (MQ=255)
                                                               cccGGGATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAACGGCAAGCACCTGACGCAGGTAACGCCAGCAggg                                >  2:1829117/1‑99 (MQ=255)
                                                                   ggATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAg                            <  1:4263026/99‑1 (MQ=255)
                                                                   ggATCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAg                            <  1:2072554/99‑1 (MQ=255)
                                                                      tCAGTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGaaa                         <  1:2510827/99‑1 (MQ=255)
                                                                         gTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAaataa                      <  1:3928673/99‑1 (MQ=255)
                                                                         gTAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAaataa                      <  1:2962707/99‑1 (MQ=255)
                                                                          tAAACAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAaataat                     >  2:3500096/1‑99 (MQ=255)
                                                                           aaaCAAGTCAGAATGCTGGAAAACGAGCCAGGCACTCAAATTTTTTCCCGCCGCGGCAAGCCCCTGACGCAGGTCACGCCCGCCGGGCCAGAcataaat                    >  2:1496664/1‑97 (MQ=255)
                                                                            aaCAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACACCAGCAGGGCAAGAAATAATTc                   <  1:2795028/99‑1 (MQ=255)
                                                                            aaCAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTc                   >  1:187509/1‑99 (MQ=255)
                                                                            aaCAAGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTc                   <  2:1203310/99‑1 (MQ=255)
                                                                            aaCAAGTAAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTc                   <  2:1643432/99‑1 (MQ=255)
                                                                               aaGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTa                <  1:1954177/99‑1 (MQ=255)
                                                                                tgTCACAATGCTGTAAAACGAGCTAGGCATTCAAATTTTTTCCCGCAGCGGCCAGCACCTGACGCAGGTAACGTCAGCAGGGCAAGAAATAATTCGTAt               <  2:4817728/98‑1 (MQ=255)
                                                                                aGTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTAt               <  2:2159879/99‑1 (MQ=255)
                                                                                 gTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTAtc              >  2:4117673/1‑99 (MQ=255)
                                                                                 gTCAGAATGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTAtc              >  1:2544819/1‑99 (MQ=255)
                                                                                     gAATGCTCTAAAACTATCTATTCATTCAAATTTTTTCCCCAATCTTCAAACACCTAACCCAAATAACACCAGCAGGGCAAGAAATAATTCGTAtcgctc          <  1:224949/99‑1 (MQ=255)
                                                                                       tctcTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTAtcgctcgc        <  2:769551/96‑1 (MQ=255)
                                                                                       aTGCTGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTAtcgctcgc        <  2:2124142/99‑1 (MQ=255)
                                                                                           tGGAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGaga    >  1:788026/1‑97 (MQ=255)
                                                                                            ggAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGAAGt   >  1:2463541/1‑99 (MQ=255)
                                                                                             gAAAACGAGCTAGGCATTCAAATTTTTTTTTTTTTTTTTAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCCAAGTc  >  2:3737997/1‑99 (MQ=255)
                                                                                             gAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGAAGTc  <  1:1502348/99‑1 (MQ=255)
                                                                                             gAAAACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGAAGTc  >  2:4263513/1‑99 (MQ=255)
                                                                                                |                                                                                               
GTTGAGGTGGTCAATCATAACCTGAATGTCTCATCAACAGCGGAAGGACTTTACACATCACAACCCGGGATCAGTAAACAAGTCAGAATGCTGGAAGACGAGCTAGGCATTCAAATTTTTTCCCGAAGCGGCAAGCACCTGACGCAGGTAACGCCAGCAGGGCAAGAAATAATTCGTATCGCTCGCGAAGTC  >  NC_000913/1333882‑1334073

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