Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,382,091 C→T 81.3% 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,0910CT81.3% 287.9 / 63.1 144W76* (TGG→TGAfrdAanaerobic fumarate reductase catalytic and NAD/flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base C (13/14);  new base T (60/57);  total (73/71)
Fisher's exact test for biased strand distribution p-value = 8.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.86e-01

Only 100 of 146 total aligned reads displayed.

TCCGGGCGACGGCTCCATGGGCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGCCCCCTTCTGCAGCAACGGTATGGCTA  >  NC_000913/4381993‑4382190
                                                                                                  |                                                                                                   
  cGGGGGACGGCTCCAGGGGCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAg                                                                                                  <  1:2205633/100‑1 (MQ=255)
  cGGGCGACGGCTCCATGGGCATCCTCACAGTTCCAGTTGTGTCATTTCGGGTGGGCAGGTGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAg                                                                                                  >  1:2328937/1‑100 (MQ=255)
  cGGGCGACGGCTCCATGGGCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAg                                                                                                  >  2:1227530/1‑100 (MQ=255)
  cGGGCGACGGCTCCATGGGCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAg                                                                                                  >  1:2199330/1‑100 (MQ=255)
   gggCGACGGCTCCAGGGGCATCCCCACCGTTCCAGTTGGGTCATTTCGGTTTGGCAGTGGTGGACGAAATAATCCGCGACATCCTGCTCACACAATCa                                                                                                   <  1:458347/98‑1 (MQ=255)
     gCGACGGCTCCATGGGCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCg                                                                                               >  2:840422/1‑100 (MQ=255)
       gACGGCTCCATGGGCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGc                                                                                              >  1:1784886/1‑99 (MQ=255)
         cGGCTCCATGGGCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCAc                                                                                           >  1:1942580/1‑100 (MQ=255)
         cGGCTCCATGGGCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCAc                                                                                           >  2:961726/1‑100 (MQ=255)
          ggCTCCATGGGCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCACACACAATCAGTCGCCAcc                                                                                          >  1:1659635/1‑100 (MQ=255)
           gCTCCATGGGCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCAccc                                                                                         >  1:2265772/1‑100 (MQ=255)
           gCTCCATGGGCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCAccc                                                                                         <  2:1049314/100‑1 (MQ=255)
           gCTCCATGGGCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCAccc                                                                                         <  2:1048890/100‑1 (MQ=255)
             tCCATGGGCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGc                                                                                       <  1:1681009/100‑1 (MQ=255)
              ccATGGGCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCt                                                                                      >  1:1408127/1‑100 (MQ=255)
                aTGGGCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTAc                                                                                    >  1:1532756/1‑100 (MQ=255)
                 tGGGCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACt                                                                                   >  1:1451670/1‑100 (MQ=255)
                 tGGGCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACt                                                                                   >  2:2314918/1‑100 (MQ=255)
                 tGGGCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGAGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACt                                                                                   >  1:1592970/1‑100 (MQ=255)
                  gggCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTg                                                                                  <  1:1625780/100‑1 (MQ=255)
                   ggCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGt                                                                                 >  2:488259/1‑100 (MQ=255)
                   ggCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGt                                                                                 >  2:1302709/1‑100 (MQ=255)
                   ggCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGt                                                                                 <  1:852733/100‑1 (MQ=255)
                   ggCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGt                                                                                 >  2:1959505/1‑100 (MQ=255)
                      aTCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATc                                                                              >  1:925356/1‑100 (MQ=255)
                       tCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCCGTGGTCGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCg                                                                             >  1:1792408/1‑100 (MQ=255)
                       tCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCg                                                                             >  1:1875113/1‑100 (MQ=255)
                       tCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCg                                                                             <  2:1646015/100‑1 (MQ=255)
                        ccccccAGTGCCAGTTGGGTCATGTCGTTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCg                                                                             <  1:63337/99‑1 (MQ=255)
                        ccccACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCg                                                                             <  1:1751963/99‑1 (MQ=255)
                        ccccACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGt                                                                            >  2:75694/1‑100 (MQ=255)
                        ccccACAGTTCCAGTTGGGTCATTTCGGTTCGGCAGTGGCGGACGAAATAATCCACGACATGCTGCTCACACAATCAGTAGCCAACCGCTACTGTATCGt                                                                            >  1:305199/1‑100 (MQ=255)
                         cccACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTg                                                                           >  1:1404254/1‑100 (MQ=255)
                          ccACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGa                                                                          >  2:1089609/1‑100 (MQ=255)
                           cacaGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGaa                                                                         <  2:2530880/100‑1 (MQ=255)
                             caGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAg                                                                       <  2:1188884/100‑1 (MQ=255)
                              aGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGt                                                                      <  2:2507097/100‑1 (MQ=255)
                                ttCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGa                                                                    <  1:441701/100‑1 (MQ=255)
                                ttCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGa                                                                    <  2:1172502/100‑1 (MQ=255)
                                ttCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGa                                                                    <  2:2376299/100‑1 (MQ=255)
                                 tCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGat                                                                   >  2:2127580/1‑100 (MQ=255)
                                 tCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGat                                                                   <  1:2525221/100‑1 (MQ=255)
                                  ccAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGata                                                                  >  2:2594237/1‑100 (MQ=255)
                                     gTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACCGTATCGTGAAAGTGATATTc                                                               <  1:319509/100‑1 (MQ=255)
                                       tGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGa                                                             >  2:656498/1‑100 (MQ=255)
                                       tGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGa                                                             >  2:1110813/1‑100 (MQ=255)
                                         ggTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAg                                                           >  2:1646935/1‑100 (MQ=255)
                                          gTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGc                                                          >  2:752697/1‑100 (MQ=255)
                                          gTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAg                                                           >  1:1007724/1‑99 (MQ=255)
                                             aTTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGt                                                       >  2:1600307/1‑100 (MQ=255)
                                              tttCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTc                                                      <  2:1464484/100‑1 (MQ=255)
                                              tttCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTc                                                      >  2:635922/1‑100 (MQ=255)
                                                  ggTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGa                                                  <  2:2034848/100‑1 (MQ=255)
                                                  ggTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGa                                                  <  2:2362278/100‑1 (MQ=255)
                                                  ggTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGa                                                  >  2:444505/1‑100 (MQ=255)
                                                     tcggCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGCCGCCACCCGCTACTGTCTCGGGAAAGTGCTATTCGAAGCTGTCATGCTcc                                               <  1:68733/98‑1 (MQ=255)
                                                     tGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATcc                                               <  2:1173789/100‑1 (MQ=255)
                                                      gggCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCt                                              >  2:1014743/1‑100 (MQ=255)
                                                       ggCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTg                                             >  1:714557/1‑100 (MQ=255)
                                                        gCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTgc                                            >  2:1803889/1‑100 (MQ=255)
                                                          aGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTgcgc                                          >  2:1103400/1‑100 (MQ=255)
                                                           gtggtgGACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGcgcg                                         <  2:468917/100‑1 (MQ=255)
                                                           gtggtgCACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGcgcg                                         <  1:1886160/100‑1 (MQ=255)
                                                             ggtggACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGAc                                       >  2:970756/1‑100 (MQ=255)
                                                             ggtggACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGAc                                       >  2:970183/1‑100 (MQ=255)
                                                               tggACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAg                                     <  2:1348991/100‑1 (MQ=255)
                                                                ggACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAgc                                    >  2:2000545/1‑100 (MQ=255)
                                                                ggACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAgc                                    >  2:2105052/1‑100 (MQ=255)
                                                                ggACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAgc                                    <  1:1086221/100‑1 (MQ=255)
                                                                 gACGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAgc                                    <  1:1768853/99‑1 (MQ=255)
                                                                  aCGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAgcgg                                  <  1:409154/100‑1 (MQ=255)
                                                                   cGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAgcggt                                 <  1:1103396/100‑2 (MQ=255)
                                                                   cGAAATAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAgcggc                                 <  1:2033325/100‑1 (MQ=255)
                                                                     aaaTAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGcggcgg                               <  2:1408135/100‑1 (MQ=255)
                                                                     aaaTAATCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGcggcgg                               >  2:160269/1‑100 (MQ=255)
                                                                      aataatCCACTACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGa                              <  1:652518/100‑1 (MQ=255)
                                                                      aataatCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGcggcgg                               <  1:1425097/99‑1 (MQ=255)
                                                                        taatCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGc                            <  1:1014740/100‑1 (MQ=255)
                                                                          atCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGccc                          <  2:2442514/100‑1 (MQ=255)
                                                                          atCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGccc                          <  2:915823/100‑1 (MQ=255)
                                                                          atCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGccc                          <  2:3103/100‑1 (MQ=255)
                                                                           tCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGcccc                         >  2:373185/1‑100 (MQ=255)
                                                                           tCCACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGGGAAAGTGATATGCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGcc                           <  1:204215/98‑1 (MQ=255)
                                                                            ccACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGccccc                        <  1:1330357/100‑1 (MQ=255)
                                                                            ccACGACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGccc                          <  1:2521490/98‑1 (MQ=255)
                                                                                gACATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGCCCCCTTCt                    >  2:270216/1‑100 (MQ=255)
                                                                                  cATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGCCCCCTTCTgc                  <  1:2482735/100‑1 (MQ=255)
                                                                                  cATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGCCCCCTTCTgc                  <  2:226914/100‑1 (MQ=255)
                                                                                  cATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGCCCCCTTCTgc                  <  1:1271193/100‑1 (MQ=255)
                                                                                  cATCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGCCCCCTTCTgc                  >  1:1534865/1‑100 (MQ=255)
                                                                                    tCCTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGCCCCCTTCTgca                 <  1:21244/99‑1 (MQ=255)
                                                                                     ccTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGCCCCCTTCTgtagc               <  2:1263486/100‑1 (MQ=255)
                                                                                     ccTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGCCCCCTTCTgcagc               <  1:1009772/100‑1 (MQ=255)
                                                                                      cTGCTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGCCCCCTTCTgcagca              <  2:1875123/100‑1 (MQ=255)
                                                                                        tcTCACACAATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGCCCCCTTCTGCAGCAAc            <  2:1815432/99‑1 (MQ=255)
                                                                                           cacacaATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGCCCCCTTCTGCAGCAACGGt         >  2:2244785/1‑100 (MQ=255)
                                                                                           cacacaATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATCATACTGCGCCACAGCGGCGGAGCCCCCTTCTGCAGCAAcaga         >  1:1815422/1‑97 (MQ=255)
                                                                                            acacaATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGCCCCCTTCTGCAGCAACGGTa        >  1:799391/1‑100 (MQ=255)
                                                                                             cacaATCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGCCCCCTTCTGCAGCAACGGTAt       >  2:794109/1‑100 (MQ=255)
                                                                                                  ccAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGTTGTCATGATCCTGCGCGACAGCGGCGGAGCCCCCTTCTGCAGCAACGGTATGGCTa  >  2:2539372/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
TCCGGGCGACGGCTCCATGGGCATCCCCACAGTTCCAGTTGGGTCATTTCGGTTGGGCAGTGGTGGACGAAATAATCCACGACATCCTGCTCACACAACCAGTCGCCACCCGCTACTGTATCGTGAAAGTGATATTCGAAGCTGTCATGATCCTGCGCGACAGCGGCGGAGCCCCCTTCTGCAGCAACGGTATGGCTA  >  NC_000913/4381993‑4382190

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