Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,764,687 T→C 15.4% intergenic (‑301/+247) sufA ← / ← ydiH Fe‑S cluster assembly protein/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,764,6870TC15.4% 365.6 / 57.1 165intergenic (‑301/+247)sufA/ydiHFe‑S cluster assembly protein/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (82/57);  new base C (12/13);  total (95/70)
Fisher's exact test for biased strand distribution p-value = 3.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.05e-01

Only 100 of 166 total aligned reads displayed.

GCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCAT  >  NC_000913/1764588‑1764785
                                                                                                   |                                                                                                   
 cATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTc                                                                                                    <  2:1993985/100‑1 (MQ=255)
   ttAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCtg                                                                                                  >  1:275490/1‑100 (MQ=255)
   ttAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCt                                                                                                   >  2:1436743/1‑99 (MQ=255)
     aaagctgaaatgaAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCtg                                                                                                  <  2:2943438/89‑1 (MQ=255)
      aGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑c                                                                                              >  1:1196556/1‑100 (MQ=255)
       ggCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGATTATTTAAGAAATATGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑Cg                                                                                             <  2:1594554/100‑1 (MQ=255)
       ggCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑Cg                                                                                             >  1:1925192/1‑100 (MQ=255)
         cTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGac                                                                                           >  2:3411718/1‑100 (MQ=255)
         cTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGac                                                                                           <  1:1097473/100‑1 (MQ=255)
         cTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGac                                                                                           >  2:1302230/1‑100 (MQ=255)
           gtaatgaAAAAGAAAAAGTTGTATTTATAATGACTATTTAATAGATAGGTAAAAAATTTCAGCGTTCAGAAATATTAAAACCGTTAAGTCtgtttg                                                                                              <  2:1340257/93‑4 (MQ=255)
             gatgaAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAg                                                                                       >  1:2608684/1‑100 (MQ=255)
               tgaAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGc                                                                                     >  1:1257163/1‑100 (MQ=255)
                gaaaaagaaaaaGTTGTATTGAAAATGACTATTTAAGAGATAGGT‑AAAAAGTGCAGCGTTAAGAAATAAAAAAACCCTTAAGTCTGTG‑CGACACAGGCt                                                                                    >  2:1192644/1‑99 (MQ=255)
                 aaaaagaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCt                                                                                    >  1:727216/1‑99 (MQ=255)
                 aaaaagaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCt                                                                                    <  2:364492/99‑1 (MQ=255)
                 aaaaagaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGATCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAg                                                                                       >  2:2453831/1‑96 (MQ=255)
                    aagaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAg                                                                                <  2:1500786/100‑1 (MQ=255)
                    aagaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAg                                                                                <  2:1736096/100‑1 (MQ=255)
                     agaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAgg                                                                               >  2:1574501/1‑100 (MQ=255)
                     agaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAgg                                                                               >  1:3450296/1‑100 (MQ=255)
                     agaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTGTCG‑CACAGACTTAAgg                                                                               <  1:513161/100‑1 (MQ=37)
                     agaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTG‑TCGCACAGACTTAAgg                                                                               <  2:1676621/100‑1 (MQ=37)
                     agaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTG‑TCGCACAGACTTAAgg                                                                               <  1:2461452/100‑1 (MQ=37)
                      gaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAggg                                                                              <  2:3279640/100‑1 (MQ=255)
                      gaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAggg                                                                              >  2:3226763/1‑100 (MQ=255)
                      gaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTGTCG‑CACAGACTTAAggg                                                                              <  1:3657292/100‑1 (MQ=37)
                      gaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTGTCG‑CACAGACTTAAggg                                                                              <  1:1507293/100‑1 (MQ=37)
                      gaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTGTCG‑CACAGACTTAAggg                                                                              <  1:3182908/100‑1 (MQ=37)
                      gaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTG‑TCGCACAGACTTAAggg                                                                              <  2:3449390/100‑1 (MQ=37)
                      gaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTAAGAAATAAGAAAACGCTTAAGTCAGTG‑CGACACAGGCTTAAgg                                                                               >  2:1260010/1‑99 (MQ=255)
                        aaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGtt                                                                            >  1:1945664/1‑100 (MQ=255)
                         aaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGttt                                                                           <  1:3834691/100‑1 (MQ=255)
                           agTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCt                                                                         <  1:3318781/100‑1 (MQ=255)
                           agTTGTATTGAAAAAGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CCACACAGGCTTAAGGGttcc                                                                          <  2:2708038/99‑3 (MQ=255)
                              tGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTAcc                                                                      <  2:462514/100‑1 (MQ=255)
                              tGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTAcc                                                                      <  1:290782/100‑1 (MQ=255)
                               gTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTAccc                                                                     <  2:3868225/100‑1 (MQ=255)
                               gTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTAccc                                                                     >  2:3387954/1‑100 (MQ=255)
                                 aTTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCa                                                                   >  1:3107960/1‑100 (MQ=255)
                                 aTTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCa                                                                   >  1:3047698/1‑100 (MQ=255)
                                    gAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATcc                                                                >  1:331736/1‑100 (MQ=255)
                                      aaaTGACTATTTAAGAGATAGGTAAATAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCgg                                                              >  2:1366626/1‑100 (MQ=255)
                                       aaTGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAATCCCTTATGTTTGAG‑CAACAGAGGCCTACCGGTTTCTACCCCATCCGgc                                                             <  2:697064/100‑1 (MQ=255)
                                        aTGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGgcg                                                            >  1:1100947/1‑100 (MQ=255)
                                         tGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGgcgc                                                           >  1:1592289/1‑100 (MQ=255)
                                           aCTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCtt                                                         <  1:759192/100‑1 (MQ=255)
                                           aCTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCtt                                                         <  1:2609265/100‑1 (MQ=255)
                                            cTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTa                                                        <  2:751924/100‑1 (MQ=255)
                                             tATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTAt                                                       >  2:2610338/1‑100 (MQ=255)
                                              aTTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTAtc                                                      <  2:190381/100‑1 (MQ=255)
                                               tttAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTAtct                                                     >  2:2847791/1‑100 (MQ=255)
                                               tttAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTAtct                                                     >  1:7254/1‑100 (MQ=255)
                                                ttAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTAtctc                                                    >  2:3226673/1‑100 (MQ=255)
                                                 tAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTcc                                                   >  1:1042298/1‑100 (MQ=255)
                                                 tAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTcc                                                   <  1:1574510/100‑1 (MQ=255)
                                                 tAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTATCTACCCCATCCGGCGCTTATCTcc                                                   >  2:484742/1‑100 (MQ=255)
                                                   agagATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCgg                                                 <  1:3853402/100‑1 (MQ=255)
                                                    gagaTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGc                                                >  1:1144725/1‑100 (MQ=255)
                                                      gaTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTGAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACTCCATCCGGCGCTTATCTCCGGCAc                                              >  1:2672301/1‑100 (MQ=255)
                                                      gaTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCAc                                              >  2:1417884/1‑100 (MQ=255)
                                                        tAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCActc                                            <  2:2061944/100‑1 (MQ=255)
                                                        tAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCActc                                            >  1:2914367/1‑100 (MQ=255)
                                                         agatAACAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGGG‑CGACACAGGCTTAAGGGTTGCGACCCCATCCGGCGCTTATCTCCGGCActc                                            <  2:847967/96‑1 (MQ=255)
                                                         aGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCActct                                           >  1:2018221/1‑100 (MQ=255)
                                                          ggTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACtctc                                          >  2:1211407/1‑100 (MQ=255)
                                                           gTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCa                                         >  1:3076088/1‑100 (MQ=255)
                                                             aacaaaGTGCAGCGTTCAGCAATGAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGt                                       <  1:2556300/97‑1 (MQ=255)
                                                             aaaaaaGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAATGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGt                                       >  2:1291934/1‑100 (MQ=255)
                                                              aaaaaGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTg                                      >  1:2798918/1‑100 (MQ=255)
                                                               aaaaGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTgg                                     <  2:3922458/100‑1 (MQ=255)
                                                                aaaGTGCAACGTTCAAAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTgg                                     <  2:2745297/99‑1 (MQ=255)
                                                                      cAGCGTTCATAAATAAGTAACCCCTTAAGTCCGTG‑CCACACCCGCTTAAGGGTTTATACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTa                                 <  2:3833850/97‑1 (MQ=255)
                                                                       aGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGctc                             <  1:1567387/100‑1 (MQ=255)
                                                                        gCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGctct                            >  1:3816364/1‑100 (MQ=255)
                                                                        gCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGctct                            >  1:1758989/1‑100 (MQ=255)
                                                                         cGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCtt                           <  2:559548/100‑1 (MQ=255)
                                                                          gTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑GGACACAGGCTTAAGGGTTTCGACCCCGTCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTg                          <  2:418632/100‑1 (MQ=255)
                                                                               gAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCATAAGGGTTTCTACCCCCTCCGGCGCTTATCTCCGGCAGTCTCAGTGGCTTAGCTCTTCAAggg                     >  2:1133792/1‑100 (MQ=255)
                                                                                aaaTAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAgggg                    >  1:87690/1‑100 (MQ=255)
                                                                                   tAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGc                   <  1:240967/98‑1 (MQ=255)
                                                                                   tAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGa                 >  2:56439/1‑100 (MQ=255)
                                                                                     aGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGTCTTAGCTCTTGAAGGGGCGATa               >  1:1819728/1‑100 (MQ=255)
                                                                                     aGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATa               <  2:500342/100‑1 (MQ=255)
                                                                                     aGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATa               >  2:1364056/1‑100 (MQ=255)
                                                                                     aGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATa               <  1:1018829/100‑1 (MQ=255)
                                                                                          aCCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaat          >  2:3912245/1‑100 (MQ=255)
                                                                                          aCCCTTAAGCCTGTGTCG‑CACAGACTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaat          >  2:1613662/1‑100 (MQ=37)
                                                                                           cccTTAAGCCTGTGTCG‑CACAGACTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaata         >  1:2734103/1‑100 (MQ=37)
                                                                                           cccTTAAGCCTGTGTCG‑CACAGACTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaa           >  2:3709274/1‑98 (MQ=37)
                                                                                           accTTAAGACTGTG‑TTGCACAGAATTAAGGTTTTCTTCCCTAGCTGGCGCTTGTCGCCGGCACTCTCAGTAGAATAGCTCTTGAAGGGGCGATAAGaata         >  2:3633752/2‑100 (MQ=255)
                                                                                            ccTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaataa        >  1:3604576/1‑100 (MQ=255)
                                                                                            ccTTAAGCCTGTGTCG‑CACAGACTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaataa        >  1:3450678/1‑100 (MQ=37)
                                                                                            ccTTAAGCCTGTG‑TCTCACAGACTTAAGGGTTTCTACACCATCAGGCGCTTATCTCCGGCACTGTCAGTGGCTTAGCTCTTGAAGGGTCGATAAGtataa        >  2:1753101/1‑100 (MQ=25)
                                                                                             cTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaataat       >  1:3402040/1‑100 (MQ=255)
                                                                                               tAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAAtct     >  2:3844424/1‑100 (MQ=255)
                                                                                               tAAGCCTGTG‑TCGCACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGCAGGGGCGATAAGAATAAtct     >  1:214171/1‑100 (MQ=38)
                                                                                                aaGTCGGTG‑CGACACAGGCTTAAGGGTTTCTACCACATTCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAAtctc    <  2:722693/100‑1 (MQ=255)
                                                                                                  gTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCAt  >  2:3168949/1‑100 (MQ=255)
                                                                                                  gTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCAt  <  1:3752928/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCAT  >  NC_000913/1764588‑1764785

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: ‑