Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 5,010 C→G P358A (CCT→GCT)  SAUSA300_RS00025 → DNA recombination protein RecF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077935,0100CG99.7% 867.4 / ‑5.9 328P358A (CCT→GCT) SAUSA300_RS00025DNA recombination protein RecF
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (314/12);  minor base A (1/0);  total (316/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.74e-01

Only 100 of 328 total aligned reads displayed.

TCAAACGCATTTATTAAGTACGATTCAGCATAAAGTACAAACATTTGTCACTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACCCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCATTGTCAGATGTAAACAACACGGATAATTATGGTG  >  NC_007793/4916‑5106
                                                                                              |                                                                                                
 gcaaCGCATTTATTAAGTACGATTCAGCATAAAGTACAAACATTTGTCACTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAAt                                                                                            >  1:3448926/3‑100 (MQ=255)
 cAAACGCATTTATTAAGTACGATTCAGCATAAAGTACAAACATTTGTCACTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTACAt                                                                                            >  1:1590839/1‑100 (MQ=255)
 cAAACGCATTTATTAAGTACGATTCAGCATAAAGTACAAACATTTGTCACTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAAt                                                                                            >  1:6034253/1‑100 (MQ=255)
 cAAACGCATTTATTAAGTACGATTCAGCATAAAGTACAAACATTTGTCACTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAAt                                                                                            >  1:656608/1‑100 (MQ=255)
   aaCGCATTTATTAAGTACGATTCAGCATAACGTACAAACATTTGTCACTACGACATCTGTCGATGGTATTGATCATGAAATCATGGATAACGCTAAAtcg                                                                                          >  1:4933403/1‑98 (MQ=255)
   aaCGCATTTATTAAGTACGATTCAGCATAAAGTACAAACATTTGTCACTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTg                                                                                          >  1:947799/1‑100 (MQ=255)
   aaCGCATTTATTAAGTACGATTCAGCATAAAGTACAAACATTTGTCACTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTg                                                                                          >  1:2474986/1‑100 (MQ=255)
               aaGTACGATTCAGCATAAAGTACAAACATTTGTCACTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAAt                                                                              >  1:6339146/1‑100 (MQ=255)
                   aCGATTCAGCATAAAGTACAAACATTTGTCACTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAg                                                                          >  1:1453411/1‑100 (MQ=255)
                      aTTCAGCATAAAGTACAAACATTTGTCACTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTg                                                                       >  1:754429/1‑100 (MQ=255)
                      aTTCAGCATAAAGTACAAACATTTGTCACTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTg                                                                       >  1:4150924/1‑100 (MQ=255)
                      aTTCAGCATAAAGTACAAACATTTGTCACTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTg                                                                       >  1:4573925/1‑100 (MQ=255)
                       ttCAGCATAAAGTACAAACATTTGTCACTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGa                                                                      >  1:1758631/1‑100 (MQ=255)
                       ttCAGCATAAAGTACAAACATTTGTCACTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGa                                                                      >  1:685354/1‑100 (MQ=255)
                        tCAGCATAAAGTACAAACATTTGTCACTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGaa                                                                     >  1:6490858/1‑100 (MQ=255)
                        tCAGCATAAAGTACAAACATTTGTCACTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGaa                                                                     >  1:267130/1‑100 (MQ=255)
                        tCAGCATAAAGTACAAACATTTGTCACTACGACATCCGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGaa                                                                     >  1:680762/1‑100 (MQ=255)
                        tCAGCATAAAGAACAAACATTTGTCACTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGGGaa                                                                     >  1:5219862/1‑100 (MQ=255)
                        acAGCATAAAGTACAAACATTTGTCACTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGaa                                                                     >  1:2410703/2‑100 (MQ=255)
                         cAGCATAAAGTACAAACATTTGTCACTACGACATCTGTAGATGGTATTGATTATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGaaa                                                                    >  1:5382765/1‑100 (MQ=255)
                         cAGCATAAAGTACAAACATTTGTCACTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGaaa                                                                    >  1:4004879/1‑100 (MQ=255)
                         cAGCATAAAGTACAAACATTTGTCACTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGaaa                                                                    >  1:3792384/1‑100 (MQ=255)
                         cAGCATAAAGTACAAACATTTGTCACTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGaaa                                                                    >  1:373843/1‑100 (MQ=255)
                         cAGCATAAAGTACAAACATTTGTCACTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGaaa                                                                    >  1:4076012/1‑100 (MQ=255)
                         cAGCATAAAGTACAAACATTTGTCACTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGaaa                                                                    >  1:4356653/1‑100 (MQ=255)
                         cAGCATAAAGTACAAACATTTGTCACTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGaaa                                                                    >  1:4502835/1‑100 (MQ=255)
                         cAGCATAAAGTACAAACATTTGTCACTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGaaa                                                                    >  1:4642451/1‑100 (MQ=255)
                         cAGCATAAAGTACAAACATTTGTCACTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGaaa                                                                    >  1:5046137/1‑100 (MQ=255)
                         cAGCATAAAGTACAAACATTTGTCACTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGaaa                                                                    >  1:5655760/1‑100 (MQ=255)
                         cAGCATAAAGTACAAACATTTGTCACTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGaaa                                                                    >  1:6296578/1‑100 (MQ=255)
                         cAGCATAAAGTACAAACATTTGTCACTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGaaa                                                                    >  1:6375039/1‑100 (MQ=255)
                         cAGCATAAAGTACAAACATTTGTCACTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGaaa                                                                    >  1:217785/1‑100 (MQ=255)
                         cAGCATAAAGTACAAACATTTGTCACTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGaaa                                                                    >  1:126733/1‑100 (MQ=255)
                         cAGCATAAAGTACAAACATTTGTCACTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGaaa                                                                    >  1:1152847/1‑100 (MQ=255)
                          aGCATAAAGTACAAACATTTGTCACTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAAt                                                                   >  1:1633849/1‑100 (MQ=255)
                          aGCATAAAGTACAAACATTTGTCACTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAAt                                                                   >  1:4770703/1‑100 (MQ=255)
                          aGCATAAAGTACAAACATTTGTCACTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAAt                                                                   >  1:575249/1‑100 (MQ=255)
                               aaaGTACAAACATTTGTCACTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAGCGCTAAATTGTATCGTATTAATCAAGGTGAAATTATaa                                                              >  1:5471035/1‑100 (MQ=255)
                                taGTACAAACATTTGTCACTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATaaa                                                             >  1:585326/2‑100 (MQ=255)
                                  gTACAAACATTTGTCACTACGACATCTGTAGATGGGATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGt                                                           >  1:6150281/1‑100 (MQ=255)
                                    aCAAACATTTGTCACTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTaa                                                         >  1:1021457/1‑100 (MQ=255)
                                     aaaaCATTTGTCACTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAAc                                                        >  1:2825204/2‑100 (MQ=255)
                                       aaCATTTGTCACTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAg                                                      <  1:3489156/100‑1 (MQ=255)
                                           gtggTCACTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAg                                                  >  1:1784390/4‑100 (MQ=255)
                                              cccACTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGa                                               >  1:3856625/3‑100 (MQ=255)
                                              cacccTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGa                                               >  1:323773/5‑100 (MQ=255)
                                                  cTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGACAGcagatcg                                           >  1:3531938/1‑94 (MQ=255)
                                                    aCGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAgatcggaaga                                         >  1:4197065/1‑91 (MQ=255)
                                                    aCGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGa                                         >  1:25474/1‑100 (MQ=255)
                                                          tCTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCAt                                   >  1:5679759/1‑100 (MQ=255)
                                                           cTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCAtt                                  >  1:1547082/1‑100 (MQ=255)
                                                           cTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCAtt                                  >  1:1919/1‑100 (MQ=255)
                                                           cTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCAtt                                  >  1:1950743/1‑100 (MQ=255)
                                                           cTGTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCAtt                                  >  1:2130039/1‑100 (MQ=255)
                                                             gTAGATGGTATTGATCATGAAATCATGAATAACGCTAGATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCATTGt                                >  1:3922657/1‑100 (MQ=255)
                                                             gTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTCTCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCATTGt                                >  1:5094404/1‑100 (MQ=255)
                                                             gTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAACGCGATGGTGACTGCATTGt                                >  1:4538294/1‑100 (MQ=255)
                                                             gTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCATTGt                                >  1:2979605/1‑100 (MQ=255)
                                                             gTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCATTGt                                >  1:2261646/1‑100 (MQ=255)
                                                             gTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCATTGt                                >  1:167561/1‑100 (MQ=255)
                                                             gTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCATTGt                                >  1:6417744/1‑100 (MQ=255)
                                                             gTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCATTGt                                >  1:5621124/1‑100 (MQ=255)
                                                             gTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCATTGt                                >  1:562397/1‑100 (MQ=255)
                                                             gTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCATTGt                                >  1:6389689/1‑100 (MQ=255)
                                                             gTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCATTGt                                >  1:1847955/1‑100 (MQ=255)
                                                             gTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCATTGt                                >  1:3663109/1‑100 (MQ=255)
                                                             gTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCATTGt                                >  1:2871431/1‑100 (MQ=255)
                                                             gTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCATTGt                                >  1:318853/1‑100 (MQ=255)
                                                             gTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCATTGt                                >  1:3419735/1‑100 (MQ=255)
                                                             gTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCATTGt                                >  1:4444599/1‑100 (MQ=255)
                                                             gTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCATTGt                                >  1:3316586/1‑100 (MQ=255)
                                                             gTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCATTGt                                >  1:6517216/1‑100 (MQ=255)
                                                             gTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCATTGt                                >  1:2327561/1‑100 (MQ=255)
                                                             gTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCATTGt                                >  1:4296162/1‑100 (MQ=255)
                                                             gTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCATTGt                                >  1:4218098/1‑100 (MQ=255)
                                                             gTAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGCGAAATTATAAAGCAACAGAAAGCGATGGTGACTGCATTGt                                >  1:4919858/1‑100 (MQ=255)
                                                             gTAGATGGCATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCATTGt                                >  1:2567788/1‑100 (MQ=255)
                                                              tAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTGTCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCATTGTc                               <  1:1123037/100‑1 (MQ=255)
                                                              tAGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCATTGTc                               <  1:3272444/100‑1 (MQ=255)
                                                               aGATGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCATTGTCa                              >  1:5920638/1‑100 (MQ=255)
                                                                 aTGGTATTGATCATGAGATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCATTGTCAGa                            >  1:895640/1‑100 (MQ=255)
                                                                 aTGGTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGAAACAGAAAGCGATGGTGACTGCATTGTCAGa                            >  1:106728/1‑100 (MQ=255)
                                                                   ggTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCATTGTCAGATg                          >  1:3946314/1‑100 (MQ=255)
                                                                   ggTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCATTGTCAGATg                          >  1:2066671/1‑100 (MQ=255)
                                                                   ggTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCATTGTCAGATg                          >  1:2190413/1‑100 (MQ=255)
                                                                    gTATTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCATTGTCAGATGt                         <  1:5287877/100‑1 (MQ=255)
                                                                      aTTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCATTGTCAGATGTaa                       >  1:5888444/1‑100 (MQ=255)
                                                                      aTTGATCATGAAATCATGAATAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCATTGTCAGATGTaa                       >  1:4489155/1‑100 (MQ=255)
                                                                           tcctgaaatcatgaaTAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCATTGTCAGATGTAaacaac                  >  1:487082/4‑100 (MQ=255)
                                                                           acatgaaatcatgaaTAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCATTGTCAGATGTAaacaac                  >  1:3060625/2‑100 (MQ=255)
                                                                             atgaaatcatgaaTAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCATTGTCAGATGTAAACAacac                <  1:1297283/100‑1 (MQ=255)
                                                                              tgaaatcatgaaTAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCATTGTCAGATGTAAACAACACg               >  1:3578734/1‑100 (MQ=255)
                                                                               gaaatcatgaaTAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCATTGTCAGATGTAAACAACACgg              >  1:4125287/1‑100 (MQ=255)
                                                                                gaatcatgaaTAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCATTGTCAGATGTAAACAACGCGGa             >  1:4394724/2‑100 (MQ=255)
                                                                                    catgaaTAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCATTGTCAGATGTAAACAACACGGATAAt         >  1:6178047/1‑100 (MQ=255)
                                                                                    catgaaTAACGCTAAATTGCATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCATTGTCAGATGTAAacaacaagatcgga         >  1:3097303/1‑92 (MQ=255)
                                                                                         aTAACGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCATTGTCAGATGTAAACAACACGGATAATTATgg    >  1:853164/1‑100 (MQ=255)
                                                                                           aaCGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTTACTGCATTGTCAGATGTAAACAACACGGATAATTATGGTg  >  1:5514110/1‑100 (MQ=255)
                                                                                           aaCGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCATTGTCAGATGTAAACAACACGGATAATTATGGTg  >  1:2024479/1‑100 (MQ=255)
                                                                                           aaCGCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCATTGTCAGATGTAAACAACACGGATAATTATGGTg  >  1:2933986/1‑100 (MQ=255)
                                                                                              |                                                                                                
TCAAACGCATTTATTAAGTACGATTCAGCATAAAGTACAAACATTTGTCACTACGACATCTGTAGATGGTATTGATCATGAAATCATGAATAACCCTAAATTGTATCGTATTAATCAAGGTGAAATTATAAAGTAACAGAAAGCGATGGTGACTGCATTGTCAGATGTAAACAACACGGATAATTATGGTG  >  NC_007793/4916‑5106

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑