Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,039,960 Δ1 bp intergenic (+42/‑36) yihD → / → yihE hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,039,9590A.99.8% 3019.3 / ‑3.5 630intergenic (+41/‑37)yihD/yihEhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (318/310);  minor base C (0/1);  total (319/311)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 668 total aligned reads displayed.

TGAAGAAGACTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATAATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACCATCATGGACGCTC  >  NC_000913/4039864‑4040053
                                                                                               |                                                                                              
tGAAGAAGACTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATa                                                                                                <  1:6517647/96‑1 (MQ=255)
 gaagaaGACTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATat                                                                                               >  2:4615298/1‑95 (MQ=255)
   agaagaCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGGAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATatcc                                                                                             >  2:4833851/1‑93 (MQ=255)
    gaagaCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATatcct                                                                                            >  2:12383349/1‑92 (MQ=255)
     aagaCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCtt                                                                                          >  1:14357369/1‑96 (MQ=255)
     aagaCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCtt                                                                                          <  1:10915048/96‑1 (MQ=255)
       gaCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCg                                                                                        >  2:53908/1‑96 (MQ=255)
        aCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTcct                                                                                       <  1:10575158/96‑3 (MQ=255)
         cTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTc                                                                                      <  1:9526824/96‑1 (MQ=255)
          ttCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCa                                                                                     <  2:11521771/96‑1 (MQ=255)
          ttCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCa                                                                                     >  2:7010821/1‑96 (MQ=255)
            cAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCAtt                                                                                   <  2:8331683/96‑1 (MQ=255)
               aaCGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGt                                                                                >  1:9937882/1‑96 (MQ=255)
                 cGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTaa                                                                              <  2:3946536/96‑1 (MQ=255)
                  ggCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAAt                                                                             >  2:7204026/1‑96 (MQ=255)
                  ggCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAAt                                                                             <  2:10973609/96‑1 (MQ=255)
                   gcgcTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATg                                                                            <  2:13077821/96‑1 (MQ=255)
                     gcTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGtt                                                                          >  1:13955658/1‑96 (MQ=255)
                       tGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGGGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGtttt                                                                        >  2:13823404/1‑96 (MQ=255)
                       tGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGtttt                                                                        >  1:12427571/1‑96 (MQ=255)
                       agTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGtttt                                                                        >  1:11306838/2‑96 (MQ=255)
                        gTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTc                                                                       >  2:3762487/1‑96 (MQ=255)
                         ttACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGtttttc                                                                      <  1:11516561/96‑3 (MQ=255)
                          tACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCg                                                                     >  2:9396601/1‑96 (MQ=255)
                           aggcgcTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCgg                                                                    >  1:9318120/3‑96 (MQ=255)
                           aCGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCgg                                                                    <  2:163583/96‑1 (MQ=255)
                            cgcgCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGga                                                                   <  2:5588536/96‑1 (MQ=255)
                              cgcTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGgatg                                                                 <  2:7458500/96‑1 (MQ=255)
                               gcTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGgatga                                                                >  1:10186959/1‑96 (MQ=255)
                               gcTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGgatga                                                                >  2:6239728/1‑96 (MQ=255)
                                cTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGgatgat                                                               <  2:6664293/96‑1 (MQ=255)
                                 tCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGatgatg                                                              <  1:8151513/96‑1 (MQ=255)
                                  cgcgGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATgg                                                             >  2:3197485/1‑96 (MQ=255)
                                   gcggcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATggg                                                            <  1:6382145/96‑1 (MQ=255)
                                    cggcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGa                                                           <  2:9209907/96‑1 (MQ=255)
                                     ggcgTAATTAAAGAGTAAAAGCTTGTAAGCGGGGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGAt                                                          >  1:6839645/1‑96 (MQ=255)
                                     ggcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGAt                                                          >  2:6013535/1‑96 (MQ=255)
                                       cgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGa                                                        <  2:505756/96‑1 (MQ=255)
                                        gTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCGCCAAAATCATCGTGAAATGATCTCCTTCGTCCCTCGTAATGTTTTCCGGATGATGGGATGaa                                                       <  1:3143086/96‑1 (MQ=255)
                                        gTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGaa                                                       >  1:9193601/1‑96 (MQ=255)
                                        gTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGaa                                                       <  1:3737864/96‑1 (MQ=255)
                                          aaTTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATa                                                     >  1:4504632/1‑96 (MQ=255)
                                          aaTTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATa                                                     >  1:2872006/1‑96 (MQ=255)
                                          aaTTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATAATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATa                                                     <  2:10427687/96‑1 (MQ=255)
                                            ttAAAGACTACAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAAc                                                   <  1:5483629/96‑1 (MQ=255)
                                             tAAAGAGTAAAAGCTTGTAAGGGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACa                                                  >  1:5629679/1‑96 (MQ=255)
                                               aaGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAgc                                                >  2:10810328/1‑96 (MQ=255)
                                                agagTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTAATTCGTAATGTTTTCCGGATCATGGGATGAATAATAgcg                                               <  1:12621715/96‑1 (MQ=255)
                                                 gagTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTTGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAgcgc                                              <  1:3880287/96‑1 (MQ=255)
                                                  agTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCt                                             >  2:13437520/1‑96 (MQ=255)
                                                   gTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCtt                                            >  1:834831/1‑96 (MQ=255)
                                                    tAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCttt                                           <  1:3517446/96‑1 (MQ=255)
                                                     aaaaGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCtttt                                          >  2:2036810/1‑96 (MQ=255)
                                                      aaaGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTTGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTa                                         <  1:4257887/96‑1 (MQ=255)
                                                       aaGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTAc                                        >  2:2275839/1‑96 (MQ=255)
                                                         gCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACtt                                      <  1:9096088/96‑1 (MQ=255)
                                                          cTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTCATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACttt                                     <  2:1167838/96‑1 (MQ=255)
                                                           ttGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACtctc                                    >  1:900186/1‑96 (MQ=255)
                                                           ttGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTc                                    <  1:40714/96‑1 (MQ=255)
                                                           ttGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTc                                    >  2:2854045/1‑96 (MQ=255)
                                                           ttGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTc                                    >  1:13047216/1‑96 (MQ=255)
                                                           ttGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGACGATGGGATGAATAACAGCGCTTTTACTTTc                                    >  2:14122888/1‑96 (MQ=255)
                                                            tGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTcc                                   <  2:14269512/96‑1 (MQ=255)
                                                            tGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTcc                                   <  2:4991468/96‑1 (MQ=255)
                                                               aaGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGa                                >  2:8107396/1‑96 (MQ=255)
                                                                aGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGac                               <  1:13263628/96‑1 (MQ=255)
                                                                 gcggcgCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGaca                              <  1:5762717/96‑1 (MQ=255)
                                                                  cggcgCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGacac                             <  2:3354439/96‑1 (MQ=255)
                                                                    gcgcCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTa                           >  1:13563611/1‑96 (MQ=255)
                                                                     cgcCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTac                          >  1:1701564/1‑96 (MQ=255)
                                                                     cgcCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTac                          <  2:8601227/96‑1 (MQ=255)
                                                                     cgcCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTac                          <  2:4175203/96‑1 (MQ=255)
                                                                      gcCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTaca                         <  2:12541344/96‑1 (MQ=255)
                                                                      gcCACCAAAATCATCGTGAAATGATATCCTTCGGCATTCGTAATGTTTTCCGGATGCTGGGATGAATAACAGCGCTTTTACTTTCCAGACACTaca                         >  1:8790947/1‑96 (MQ=255)
                                                                       ccaccaAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTacac                        >  1:6093390/1‑96 (MQ=255)
                                                                         accaAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACAccc                      >  1:1252728/1‑96 (MQ=255)
                                                                         accaAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACAccc                      <  1:7401236/96‑1 (MQ=255)
                                                                         accaAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACAccc                      >  2:11147496/1‑96 (MQ=255)
                                                                         accaAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACAccc                      <  1:12174751/96‑1 (MQ=255)
                                                                         accaAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACAccc                      >  2:2679945/1‑96 (MQ=255)
                                                                          ccaAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCg                     <  2:2467166/96‑1 (MQ=255)
                                                                          ccaAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCg                     >  1:257355/1‑96 (MQ=255)
                                                                            aaaaTCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCggc                   >  1:7835180/1‑95 (MQ=255)
                                                                            aaaaTCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGa                   >  2:11964114/1‑96 (MQ=255)
                                                                            aaaaTCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGa                   >  1:4700442/1‑96 (MQ=255)
                                                                                  atcGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACcat             >  2:433265/1‑96 (MQ=255)
                                                                                  atcGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACcat             <  1:7010821/96‑1 (MQ=255)
                                                                                   tcGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACcatc            >  2:6515143/1‑96 (MQ=255)
                                                                                    cGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACcatca           >  1:8611463/1‑96 (MQ=255)
                                                                                    cGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACcatca           >  2:10185828/1‑96 (MQ=255)
                                                                                     gTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACcatcat          >  1:4957807/1‑96 (MQ=255)
                                                                                     gTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACcatcat          <  1:2434146/96‑1 (MQ=255)
                                                                                     gTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACcatcat          <  2:892533/96‑1 (MQ=255)
                                                                                     gTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACcatcat          >  1:3398659/1‑96 (MQ=255)
                                                                                       gAAATGATATCCTTCGTCATTCGTAATGTTTACCGGATGATGGGATGAATAACAGCGCTTTTCCTTTCCAGACACTACACCCGGATACCATCATgg        <  2:7980088/96‑1 (MQ=255)
                                                                                         aaTGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACCATCATggtc      <  1:11823971/96‑3 (MQ=255)
                                                                                         aaTGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACCATCATGGAc      <  2:13004212/96‑1 (MQ=255)
                                                                                         aaTGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACCATCATGGAc      <  1:100023/96‑1 (MQ=255)
                                                                                         aaTGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACCATCATGGAc      <  1:1960391/96‑1 (MQ=255)
                                                                                          aTGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACCATCATGGACg     >  2:9562791/1‑96 (MQ=255)
                                                                                               |                                                                                              
TGAAGAAGACTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATAATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACCATCATGGACGCTC  >  NC_000913/4039864‑4040053

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: