Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,382,091 C→T 22.8% W76* (TGG→TGA frdA ← anaerobic fumarate reductase catalytic and NAD/flavoprotein subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,382,0910CT22.8% 313.1 / 88.0 166W76* (TGG→TGAfrdAanaerobic fumarate reductase catalytic and NAD/flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base C (58/70);  new base T (25/13);  total (83/83)
Fisher's exact test for biased strand distribution p-value = 4.13e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.49e-01

Only 100 of 171 total aligned reads displayed.

ATCCGGGCGACGGCTCCATGGGCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGCCCCCTTCTGCAGCAACGGTATGGCTA  >  NC_000913/4381992‑4382190
                                                                                                   |                                                                                                   
aTCCGGGCGACGGCTCCATGGGCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAAc                                                                                                     <  2:2620269/100‑1 (MQ=255)
aTCCGGGCGACGGCTCCATGGGCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAAc                                                                                                     >  1:339047/1‑100 (MQ=255)
   cGGGCGACGGCTCCATGGGCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCTACGACATCCTGCTCACACAACCAg                                                                                                  >  1:3462018/1‑100 (MQ=255)
   cGGGCGACGGCTCCATGGGCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAg                                                                                                  >  2:2036947/1‑100 (MQ=255)
   cGGGCGACGGCTCCATGGGCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAg                                                                                                  >  2:1721564/1‑100 (MQ=255)
   cGGGCGACGGCTCCATGGGCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAg                                                                                                  <  1:3297209/100‑1 (MQ=255)
   cGGGCGACGGCTCCATGGGCATCCCAACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAAACCACGACATCCTGATCACACAACCAg                                                                                                  >  1:2005780/1‑100 (MQ=255)
    gggCGACGGCTCCATGGGCATCCCCACAGTTCCAGGGGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACTACAGCCTGCTCACACAACCAGt                                                                                                 >  1:3079423/1‑100 (MQ=255)
     ggCGACGGCTCCATGGGCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTc                                                                                                <  1:3847221/100‑1 (MQ=255)
      gCGACGGCTCCATGGGCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCg                                                                                               <  1:4341848/100‑1 (MQ=255)
        gACGGCTCCATGGGCATCCCCACAGTTCCAGTGGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGcc                                                                                             <  1:1308678/100‑1 (MQ=255)
            gCTCCATGGGCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCAccc                                                                                         >  2:3233216/1‑100 (MQ=255)
             cTCCATGGGCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCg                                                                                        >  2:2162866/1‑100 (MQ=255)
               ccATGGGCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCt                                                                                      <  2:3344769/100‑1 (MQ=255)
                cATGGGCATCCCCACAGTTCCAGTTGGGTCATTTGGGTTGGACCGTCGTGGACAAAGAACTGCCCGACATCCTGCTCGCACAATCAGACGcagctcac                                                                                       >  1:3542573/1‑91 (MQ=255)
                 aTGGGCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTAc                                                                                    <  1:218974/100‑1 (MQ=255)
                  tGGGCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACt                                                                                   <  2:213223/100‑1 (MQ=255)
                  tGGGCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACt                                                                                   >  1:2680374/1‑100 (MQ=255)
                   gggCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTg                                                                                  >  2:2903194/1‑100 (MQ=255)
                   gggCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTg                                                                                  <  1:3225472/100‑1 (MQ=255)
                    ggCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGt                                                                                 >  2:2250813/1‑100 (MQ=255)
                    ggCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGt                                                                                 >  2:3481481/1‑100 (MQ=255)
                    ggCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGt                                                                                 >  2:3695045/1‑100 (MQ=255)
                     gCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTa                                                                                <  1:2442275/100‑1 (MQ=255)
                     gCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTa                                                                                <  2:1295680/100‑1 (MQ=255)
                      cATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTAt                                                                               <  2:456284/100‑1 (MQ=255)
                       aTCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATc                                                                              <  2:3369208/100‑1 (MQ=255)
                        tCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCg                                                                             >  1:1708186/1‑100 (MQ=255)
                        tCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCg                                                                             <  1:438101/100‑1 (MQ=255)
                        tCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCg                                                                             <  2:1381280/100‑1 (MQ=255)
                         ccccACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGt                                                                            <  2:2641455/100‑1 (MQ=255)
                         ccccACAGTTCCAATTGGGTCATTTCGGTTGGGCAGTGCTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGt                                                                            >  1:2723010/1‑100 (MQ=255)
                            cacaGTTGCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGGAATAATCCGCGACATCCTGCTACCACAACCAGTCGCCACGCGCTACTGTATCGTGa                                                                          >  1:2763419/1‑99 (MQ=255)
                            cacaGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGATGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGCATCGTGaa                                                                         <  1:2115213/100‑1 (MQ=255)
                              caGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGaaa                                                                        >  1:2678260/1‑99 (MQ=255)
                              caGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAg                                                                       <  1:110778/100‑1 (MQ=255)
                                 tgccACTTGGGTCATTGCGGTTGGGCATTGCTGGACGAAATAATCCACCACATCCTGCTCACACAACCAGTCACCACCCGCTACTGTATCGTTAAAGTGa                                                                    <  1:4015723/98‑1 (MQ=255)
                                  tCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTTCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGa                                                                    <  1:3747200/99‑1 (MQ=255)
                                  tCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGat                                                                   >  2:3849086/1‑100 (MQ=255)
                                   ccAGTTGGGTGATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCTGTTGCCACCCGCTACTGTATCGTGAAAGTGata                                                                  <  1:2903177/100‑1 (MQ=255)
                                    cAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGatat                                                                 >  1:2963237/1‑100 (MQ=255)
                                      gTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTc                                                               <  2:2005787/100‑1 (MQ=255)
                                        tGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGa                                                             >  2:2790875/1‑100 (MQ=255)
                                        tGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGa                                                             <  2:830351/100‑1 (MQ=255)
                                         gggTCATTTCGGTTGGGCAGTGGTGGACGAACTAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGatatn                                                                <  1:1972024/96‑2 (MQ=255)
                                          ggTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAg                                                           >  1:2485496/1‑100 (MQ=255)
                                           gTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGaa                                                            <  1:2793474/98‑1 (MQ=255)
                                            tCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGc                                                          <  1:1240398/99‑1 (MQ=255)
                                             cATTTCGGTTGGGAAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGc                                                          <  1:1503990/98‑1 (MQ=255)
                                              aTTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGt                                                       <  2:507344/100‑1 (MQ=255)
                                                 tCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCAt                                                    >  1:1874369/1‑100 (MQ=255)
                                                 tCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCAt                                                    <  2:1758599/100‑1 (MQ=255)
                                                  cGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATg                                                   >  2:110213/1‑100 (MQ=255)
                                                   ggTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGa                                                  >  2:2427528/1‑100 (MQ=255)
                                                   ggTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGa                                                  >  2:1120134/1‑100 (MQ=255)
                                                    gTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGAt                                                 >  1:2526015/1‑100 (MQ=255)
                                                       gggCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATcc                                               >  2:1803588/1‑99 (MQ=255)
                                                       gggCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCt                                              >  2:378969/1‑100 (MQ=255)
                                                            gtggtgGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTgcgag                                         >  2:2377039/1‑98 (MQ=255)
                                                            gtggtgGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGcgcg                                         <  2:2364058/100‑1 (MQ=255)
                                                                tggACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACa                                      >  1:2772595/1‑99 (MQ=255)
                                                                  gACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAgcg                                   <  1:3650483/100‑1 (MQ=255)
                                                                    cGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAgcggc                                 >  2:2206350/1‑100 (MQ=255)
                                                                     gAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAgcggcg                                <  2:390565/100‑1 (MQ=255)
                                                                      aaaTAATCCACGACATCCTGCTCACACAATCAGTCGCTACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGcggcgg                               >  1:2357709/1‑100 (MQ=255)
                                                                      aaaTAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGcggcgg                               >  2:1377981/1‑100 (MQ=255)
                                                                      aaaTAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGcggcgg                               <  2:2755931/100‑1 (MQ=255)
                                                                       aagaatCCACGTCATCCTGCGCACACAACCAGTCGCCACCCGCGACTGTATCGTGAAAGTGATATTCGAAGCTGGCATGAGCCTGCGCGACAGcggcgg                               <  1:394094/96‑1 (MQ=255)
                                                                        ataatCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAg                             <  2:3289103/100‑1 (MQ=255)
                                                                         taatCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGc                            <  1:1005364/100‑1 (MQ=255)
                                                                          aatCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGc                            <  1:1046654/99‑1 (MQ=255)
                                                                           atCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGccc                          <  1:2194030/100‑1 (MQ=255)
                                                                           atCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGccc                          <  2:4047252/100‑1 (MQ=255)
                                                                            tCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGcccc                         >  1:2874527/1‑100 (MQ=255)
                                                                            tCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTACCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGcccc                         >  1:1800380/1‑100 (MQ=255)
                                                                             ccACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGccccc                        <  1:4317138/100‑1 (MQ=255)
                                                                              cACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGCCCCCt                       >  1:2470218/1‑100 (MQ=255)
                                                                                cGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCTTGCGCGACAGCGGCGGAGCCCCCTTc                     >  2:2576398/1‑100 (MQ=255)
                                                                                cGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGCCCCCTTc                     <  2:1954425/100‑1 (MQ=255)
                                                                                 gACATCCTGCTCACACAATCAGTCTCCAGCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGCCCCCTTCt                    <  1:4171824/100‑1 (MQ=255)
                                                                                  aCATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGCCCCCTTCTg                   >  1:3623127/1‑100 (MQ=255)
                                                                                   cATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTGTCGTGAAAGTGTTAGTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGcccc                         <  1:1986348/93‑1 (MQ=255)
                                                                                     tCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGCCCCCTTCTgcag                <  2:1603902/100‑1 (MQ=255)
                                                                                       cTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGCCCCCTTCTgcagca              >  1:1772596/1‑100 (MQ=255)
                                                                                       cTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGCCCCCTTCTgcagca              >  1:138839/1‑100 (MQ=255)
                                                                                        tGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGCCCCCTTCTGCAGCaa             >  1:3791583/1‑100 (MQ=255)
                                                                                         gCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGCCCCCTTCTGCAGCAAc            >  1:799604/1‑100 (MQ=255)
                                                                                         gCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGCCCCCTTCTGCAGCAAc            <  2:3517114/100‑1 (MQ=255)
                                                                                         gCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGCCCCCTTCTGCAGCAAc            <  2:1263320/100‑1 (MQ=255)
                                                                                         gCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGCCCCCTTCTGCAGCAAc            <  1:3494240/100‑1 (MQ=255)
                                                                                          cTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGCCCCCTTCTGCAGCAACg           <  2:3556912/100‑1 (MQ=255)
                                                                                           tCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGCCCCCTTCTGCAGCAACgg          <  1:577307/100‑1 (MQ=255)
                                                                                           tCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCCGCGCGACAGCGGCGGAGCCCCCTTCTGCAGCAACgn          <  2:3101853/100‑2 (MQ=255)
                                                                                            cacacaACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGCCCCCTTCTGCAGCAACgg          >  1:1422283/1‑99 (MQ=255)
                                                                                            cacacaACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGCCCCCTTCTGCAGCAACGGt         >  2:4045040/1‑100 (MQ=255)
                                                                                             acacaATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGCCCCCTTCTGCAGCAACGGTa        >  2:4324439/1‑100 (MQ=255)
                                                                                              cacaACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGCCCCCTTCTGCAGCAACGGTAt       >  1:1154336/1‑100 (MQ=255)
                                                                                              cacaACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGCCCCCTTCTGCAGCAACGGTAt       >  2:594645/1‑100 (MQ=255)
                                                                                              cacaACCAGTCGCCACCCGCTACTGGATCGTGAAAGTGATATTCGAAACTGTCATGATCCTGCGCGACAGCGGCGGAGCCCCCTTCTGCAGCAACGGTAt       >  1:320493/1‑100 (MQ=255)
                                                                                                caACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGCCCCCTTCTGCAGCAACGGTATgg     >  2:1987081/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
ATCCGGGCGACGGCTCCATGGGCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGCCCCCTTCTGCAGCAACGGTATGGCTA  >  NC_000913/4381992‑4382190

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: