Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,764,701 G→A 18.2% intergenic (‑315/+233) 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,7010GA18.2% 266.1 / 49.5 121intergenic (‑315/+233)sufA/ydiHFe‑S cluster assembly protein/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (57/42);  new base A (7/15);  total (64/57)
Fisher's exact test for biased strand distribution p-value = 3.48e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.21e-01

Only 100 of 124 total aligned reads displayed.

ATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCG  >  NC_000913/1764602‑1764800
                                                                                                    |                                                                                                    
 tgaAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGc                                                                                                     >  2:1364973/1‑100 (MQ=255)
  gaaaaagaaaaaGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCt                                                                                                    >  1:1480521/1‑100 (MQ=255)
     aaagaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTaa                                                                                                 <  2:2109487/100‑1 (MQ=255)
       agaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAgg                                                                                               <  1:1032571/100‑1 (MQ=255)
       agaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAgg                                                                                               >  2:1969671/1‑100 (MQ=255)
       agaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTGTCG‑CACAGACTTAAgg                                                                                               <  2:1219737/100‑1 (MQ=37)
       agaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTGTCG‑CACAGACTTAAgg                                                                                               <  1:2418578/100‑1 (MQ=37)
       agaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTGTCG‑CACAGACTTAAgg                                                                                               <  2:2168332/100‑1 (MQ=37)
       agaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTGTCG‑CACAGACTTAAgg                                                                                               <  1:1195683/100‑1 (MQ=37)
        gaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTGTCG‑CACAGACTTAAggg                                                                                              <  1:668406/100‑1 (MQ=37)
        gaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTGTCG‑CACAGACTTAAggg                                                                                              <  2:2772566/100‑1 (MQ=37)
        gaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTG‑TCGCACAGACTTAAggg                                                                                              <  2:1051797/100‑1 (MQ=37)
        gaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTG‑TCGCACAGACTTAAggg                                                                                              <  2:2385261/100‑1 (MQ=37)
        gaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGCTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTaa                                                                                                 >  1:1760031/1‑97 (MQ=255)
         aaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑t                                                                                            >  1:129715/1‑100 (MQ=255)
          aaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTGTCG‑CACAGACTTAAGGG‑tt                                                                                           <  2:800579/100‑1 (MQ=37)
           aaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTGTCG‑CACAGACTTAAGGG‑ttt                                                                                          <  2:681117/100‑1 (MQ=25)
               ttGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTAc                                                                                      >  2:795163/1‑100 (MQ=255)
                tGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTAcc                                                                                     <  2:1163826/100‑1 (MQ=255)
                  tATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTAcccc                                                                                   >  1:436299/1‑100 (MQ=255)
                  tATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTAcccc                                                                                   >  1:1227202/1‑100 (MQ=255)
                  tATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTGTCG‑CACAGACTTAAGGG‑ttttcttatt                                                                                   <  1:2275342/100‑8 (MQ=21)
                   aTTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCa                                                                                  >  1:1567763/1‑100 (MQ=255)
                   aTTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCa                                                                                  >  1:581334/1‑100 (MQ=255)
                   aTTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCa                                                                                  >  1:975718/1‑100 (MQ=255)
                   aTTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCa                                                                                  >  1:2736947/1‑100 (MQ=255)
                   aTTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCa                                                                                  >  1:2219874/1‑100 (MQ=255)
                    ttGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCAt                                                                                 >  2:2782594/1‑100 (MQ=255)
                    tagAAAATGACTAATTAAGAGATAGGTAAAAAAGTGCAGCGTACAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCAt                                                                                 >  2:1976392/3‑100 (MQ=255)
                     tGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTGTCG‑CACAGACTTAAGGG‑ttttcttatttct                                                                                <  2:676990/100‑11 (MQ=18)
                         aaTGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGgc                                                                            <  2:581337/100‑1 (MQ=255)
                          aTGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGgcg                                                                           <  2:1689449/100‑1 (MQ=255)
                           tGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGgcgc                                                                          >  1:307704/1‑100 (MQ=255)
                            gACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCt                                                                         <  2:1554978/100‑1 (MQ=255)
                             aCTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACACTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCtt                                                                        >  1:2553567/1‑100 (MQ=255)
                              cTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTa                                                                       <  2:1923181/100‑1 (MQ=255)
                              cTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTa                                                                       <  2:1300216/100‑1 (MQ=255)
                              cTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTa                                                                       <  2:1790275/100‑1 (MQ=255)
                              cTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTa                                                                       <  1:1105637/100‑1 (MQ=255)
                                aTTTGAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCTCGTAtc                                                                     >  2:2188081/1‑100 (MQ=255)
                                aTTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTAtc                                                                     >  1:539626/1‑100 (MQ=255)
                                aTTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTAtc                                                                     >  2:39846/1‑100 (MQ=255)
                                   tAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTcc                                                                  >  1:1994092/1‑100 (MQ=255)
                                   tAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTcc                                                                  >  1:2493418/1‑100 (MQ=255)
                                    aaGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTATCCCATCCGGCGCTTATCTCCg                                                                 <  1:265064/100‑1 (MQ=255)
                                     agagATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCgg                                                                <  1:628033/100‑1 (MQ=255)
                                        gaTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCAc                                                             >  2:1652584/1‑100 (MQ=255)
                                           gatctaaaaaGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCActcn                                                          <  2:1320636/95‑2 (MQ=255)
                                            ggTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACtctc                                                         >  1:2577993/1‑100 (MQ=255)
                                             gTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCa                                                        <  2:106859/100‑1 (MQ=255)
                                             gTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCa                                                        >  2:1843674/1‑100 (MQ=255)
                                             gTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCa                                                        >  2:708711/1‑100 (MQ=255)
                                             gTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑ATTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCa                                                        >  2:834526/1‑100 (MQ=255)
                                              tAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAg                                                       <  1:834523/100‑1 (MQ=255)
                                              tAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAg                                                       >  1:1881937/1‑100 (MQ=255)
                                                naaaaGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCActctagatc                                                     >  1:1320631/2‑95 (MQ=255)
                                                 aaaaGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑ATTCTACGCCATCCGGCGCTTATCTCCGGCACTCTCAGTgg                                                    >  1:466780/1‑100 (MQ=255)
                                                          ccgatctgAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGctct                                           <  2:2899454/93‑1 (MQ=255)
                                                            gTTCATAAATAAGAAAACCCTTAAGTCAGTG‑CGACAAAGGCTTAAGGT‑TTTCTACGCCATCCGGCGCTTATCTCCGTCACTCTCAGTTACTTAGCTCTTg                                         <  1:2315138/100‑1 (MQ=255)
                                                            attCAGAAATAAGAAAACACTTAAGTCTGTG‑CGACACAGGATTAAGGG‑ATTCTACACCAAAGGGCGCTTAACTACGGCACTGTCAGTGGCTTAActgt                                           >  1:1408935/2‑96 (MQ=255)
                                                              tCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGaa                                       <  2:1953559/100‑1 (MQ=255)
                                                               cAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAg                                      >  2:462891/1‑100 (MQ=255)
                                                               cAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAg                                      >  1:2058803/1‑100 (MQ=255)
                                                               cAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAg                                      >  2:1557753/1‑100 (MQ=255)
                                                                 gAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCTTCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAggg                                    <  1:708688/100‑1 (MQ=255)
                                                                 gAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGctctagatcgg                                    >  1:2899435/1‑93 (MQ=255)
                                                                 gAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAggg                                    >  2:194894/1‑100 (MQ=255)
                                                                   aaTAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGc                                  <  1:2382659/100‑1 (MQ=255)
                                                                   aaTAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGc                                  >  2:2356893/1‑100 (MQ=255)
                                                                    ttccgatcccccTTAAGCCTGTG‑TCGCACAGACTTAAGGGTTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGc                                  <  2:694502/91‑1 (MQ=21)
                                                                       aGAAAACCCTGAAGTCTGTG‑CGACACAGGCTGAAGGG‑TTTCTACCCCATCCGGGGCATATGTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATa                              <  1:2801541/100‑1 (MQ=255)
                                                                         aaaaCCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAg                            <  2:1853816/100‑1 (MQ=255)
                                                                          aaaCCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGa                           <  1:462888/100‑1 (MQ=255)
                                                                           aaCCCTTGAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCGACCCCAGCCGGCGCTTATCGCCGGCACTCGCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaa                          <  1:386191/100‑1 (MQ=255)
                                                                           aaCCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaa                          <  2:1490620/100‑1 (MQ=255)
                                                                           aaCCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaa                          >  2:1163599/1‑100 (MQ=255)
                                                                           aaCCCTTAAGCCTGTGTCG‑CACAGACTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaa                          >  2:1635889/1‑100 (MQ=37)
                                                                            aCCCTTAAGCCTGTGTCG‑CACAGACTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCAGTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaat                         >  2:2778896/1‑100 (MQ=25)
                                                                            aCCCTTAAGCCTGTG‑TCGCACAGACTTAAGGGTTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaa                          >  2:573577/1‑100 (MQ=25)
                                                                             cccTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaata                        >  2:2160650/1‑100 (MQ=255)
                                                                             cccTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaata                        <  2:1352898/100‑1 (MQ=255)
                                                                             cccTTAAGCCTGTGTCG‑CACAGACTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaata                        >  2:1801245/1‑100 (MQ=37)
                                                                             cccTTAAGCCTGTGTCG‑CACAGACTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaata                        >  1:2008406/1‑100 (MQ=37)
                                                                             cccTTAAGCCTGTG‑TCGCACAGACTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaat                         >  1:824994/1‑99 (MQ=37)
                                                                              ccTCAAGTCTGTG‑CGACACAGGCTTAAGGG‑TACATACCCCATCCGGCGCTTATCTCCGGCACCCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaat                         <  1:994326/98‑1 (MQ=255)
                                                                               cTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaataa                       <  1:1955624/99‑1 (MQ=255)
                                                                                 tAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAAtct                    <  2:1898033/100‑1 (MQ=255)
                                                                                   aGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCa                  <  2:2429940/100‑1 (MQ=255)
                                                                                   aGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCa                  <  2:359337/100‑1 (MQ=255)
                                                                                       tgtg‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATaaa              <  1:1668217/100‑1 (MQ=255)
                                                                                        gtg‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAg             <  2:1412572/100‑1 (MQ=255)
                                                                                         tg‑CGACACGGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAgn            <  1:2298884/100‑2 (MQ=255)
                                                                                         tg‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGc            >  1:2144052/1‑100 (MQ=255)
                                                                                             gACACAGGCTTAAGGG‑TTTCTACACCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGc            <  1:1007364/97‑1 (MQ=255)
                                                                                              acacAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAAc        >  2:2089350/1‑100 (MQ=255)
                                                                                               cacaGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAAcc       >  1:1612940/1‑100 (MQ=255)
                                                                                                acaGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAAccc      >  2:2580210/1‑100 (MQ=255)
                                                                                                 caGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACccg     >  1:732595/1‑100 (MQ=255)
                                                                                                   ggCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACccgcc   >  2:1765002/1‑100 (MQ=255)
                                                                                                    gCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACccgccg  >  1:588560/1‑100 (MQ=255)
                                                                                                    |                                                                                                    
ATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCG  >  NC_000913/1764602‑1764800

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