Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,867,186 T→C D68D (GAT→GAC yidS → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,867,1860TC86.0% 618.7 / 92.2 265D68D (GAT→GACyidShypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (18/19);  new base C (118/110);  total (136/129)
Fisher's exact test for biased strand distribution p-value = 8.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.65e-01

Only 100 of 265 total aligned reads displayed.

TGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTAC  >  NC_000913/3867095‑3867276
                                                                                           |                                                                                          
 ggATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACgg                                                                                          >  1:1536559/1‑93 (MQ=255)
 ggATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACgg                                                                                          >  2:2946778/1‑93 (MQ=255)
    tAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCt                                                                                       >  1:4072980/1‑93 (MQ=255)
    tAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCt                                                                                       >  2:2567193/1‑93 (MQ=255)
      aaaaGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGa                                                                                     >  1:3480638/1‑93 (MQ=255)
       aaaGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGAc                                                                                    >  1:7487156/1‑93 (MQ=255)
            aCCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTc                                                                               <  2:3698489/93‑1 (MQ=255)
            aCCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTc                                                                               >  2:3406366/1‑93 (MQ=255)
              cAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCt                                                                             <  2:3170303/93‑1 (MQ=255)
               aGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTg                                                                            <  1:3712930/93‑1 (MQ=255)
               aGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTg                                                                            >  1:5779980/1‑93 (MQ=255)
               aGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTg                                                                            <  2:1778796/93‑1 (MQ=255)
                gtgtGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGt                                                                           <  2:4096316/93‑1 (MQ=255)
                gtgtGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGt                                                                           >  1:6845659/1‑93 (MQ=255)
                 tgtgGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTc                                                                          <  2:3276430/93‑1 (MQ=255)
                 tgtgGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTc                                                                          >  2:6117423/1‑93 (MQ=255)
                   tgGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGa                                                                        <  2:4682313/93‑1 (MQ=255)
                     gTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGAtg                                                                      <  2:2366964/93‑1 (MQ=255)
                     gTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGAtg                                                                      >  2:1299641/1‑93 (MQ=255)
                     gTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGAtg                                                                      <  1:6518306/93‑1 (MQ=255)
                        cTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGa                                                                   >  1:2832255/1‑93 (MQ=255)
                        cTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGa                                                                   >  1:1421465/1‑93 (MQ=255)
                         tGAAGGTTTCAGCAACCCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGAt                                                                  <  1:3654964/93‑1 (MQ=255)
                         tGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGAt                                                                  <  1:6644842/93‑1 (MQ=255)
                         tGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGAt                                                                  >  1:822365/1‑93 (MQ=255)
                         tGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGAt                                                                  <  1:5563900/93‑1 (MQ=255)
                         tGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGAt                                                                  >  2:878762/1‑93 (MQ=255)
                          gAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATc                                                                 <  2:1423499/93‑1 (MQ=255)
                          gAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATc                                                                 >  1:3799570/1‑93 (MQ=255)
                           aaGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCg                                                                >  2:1802878/1‑93 (MQ=255)
                           aaGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCg                                                                >  2:4003033/1‑93 (MQ=255)
                            aGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGc                                                               <  1:4972998/93‑1 (MQ=255)
                            aGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGc                                                               >  1:3332324/1‑93 (MQ=255)
                             ggTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCt                                                              >  1:5368915/1‑93 (MQ=255)
                               tttCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCGGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTaa                                                            <  1:2972890/93‑1 (MQ=255)
                                ttCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAAt                                                           <  1:2497352/93‑1 (MQ=255)
                                 tCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATc                                                          <  1:5194652/93‑1 (MQ=255)
                                 tCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATc                                                          <  2:4294867/93‑1 (MQ=255)
                                   aGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCg                                                        >  2:5777418/1‑93 (MQ=255)
                                    gCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGc                                                       >  1:2563674/1‑93 (MQ=255)
                                     cAACCCCTGCGGCGGCCGGCTGCCCCCGGACGCGCAGCGATCTTTTATTCGCGCCGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCa                                                      <  1:2366567/93‑1 (MQ=255)
                                     cAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCa                                                      >  1:1792475/1‑93 (MQ=255)
                                     cAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCa                                                      <  1:6173128/93‑1 (MQ=255)
                                       aaCCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGa                                                    >  2:4181353/1‑93 (MQ=255)
                                        aCCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGAt                                                   >  1:4703280/1‑93 (MQ=255)
                                         cccTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGAtt                                                  >  1:3085610/1‑93 (MQ=255)
                                          ccTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGAttt                                                 >  1:1699766/1‑93 (MQ=255)
                                          ccTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGAttt                                                 >  1:1097712/1‑93 (MQ=255)
                                          ccTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGAttt                                                 >  1:4651106/1‑93 (MQ=255)
                                          ccTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGAttt                                                 <  1:4741485/93‑1 (MQ=255)
                                           cTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGAtttt                                                <  2:6721380/93‑1 (MQ=255)
                                           cTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGAtttt                                                >  2:506984/1‑93 (MQ=255)
                                            tGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTc                                               <  2:2010428/93‑1 (MQ=255)
                                               ggcggTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGc                                            <  2:2510836/93‑1 (MQ=255)
                                                  ggTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCgcc                                         <  1:7233095/93‑3 (MQ=255)
                                                   gTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGCCa                                        <  2:621840/93‑1 (MQ=255)
                                                       gctgGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGGCGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGCCAAAAc                                    <  2:499598/93‑1 (MQ=255)
                                                       gctgGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGCCAAAAc                                    >  2:7263164/1‑93 (MQ=255)
                                                       gctgGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGCCAAAAc                                    >  1:2690703/1‑93 (MQ=255)
                                                        ctgGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGCCAAAAcc                                   <  1:5603708/93‑1 (MQ=255)
                                                         tgGCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGCCAAAACCg                                  >  1:739366/1‑93 (MQ=255)
                                                          ggCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGt                                 <  2:4620654/93‑1 (MQ=255)
                                                          ggCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTTAGCGCCAAAACCGt                                 <  2:6446407/93‑1 (MQ=255)
                                                          ggCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGCCAAAACCGt                                 <  2:2311735/93‑1 (MQ=255)
                                                          ggCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGCCAAAACCGt                                 <  2:3110769/93‑1 (MQ=255)
                                                           gCACCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGCCAAAACCGTc                                >  2:6203157/1‑93 (MQ=255)
                                                             aCCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGCCAAAACCGTCGa                              >  2:5475905/1‑93 (MQ=255)
                                                             aCCGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGCCAAAACCGTCGa                              <  2:380313/93‑1 (MQ=255)
                                                              ccGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGCCAAAACCGTCGAc                             >  2:156953/1‑93 (MQ=255)
                                                              ccGGACGCGCAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGCCAAAACCGTCGAc                             <  1:6914654/93‑1 (MQ=255)
                                                                    gcgcAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGCCAAAACCGTCGACGTcgcc                       <  2:1338577/93‑1 (MQ=255)
                                                                    gcgcAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGCCAAAACCGTCGACGTcgcc                       <  1:6982369/93‑1 (MQ=255)
                                                                      gcAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTcgccgc                     >  1:1090071/1‑93 (MQ=255)
                                                                      gcAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTcgccgc                     >  2:1859334/1‑93 (MQ=255)
                                                                      gcAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTcgccgc                     >  1:2055760/1‑93 (MQ=255)
                                                                      gcAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGCCAAAACCGTCGACGTcgccgc                     >  2:6363492/1‑93 (MQ=255)
                                                                      gcAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGCCAAAACCGTCGACGTcgccgc                     >  1:5964321/1‑93 (MQ=255)
                                                                      gcAGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGCCAAAACCGTCGACGTcgccgc                     >  2:288235/1‑93 (MQ=255)
                                                                        aGCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGCCAAAACCGTCGACGTCGCCGCAt                   <  1:3411977/93‑1 (MQ=255)
                                                                         gCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATc                  >  2:6922619/1‑93 (MQ=255)
                                                                         gCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGCCAAAACCGTCGACGTCGCCGCATc                  >  1:5151761/1‑93 (MQ=255)
                                                                         gCGATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGCCAAAACCGTCGACGTCGCCGCATc                  >  2:5995330/1‑93 (MQ=255)
                                                                           gATCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGc                <  2:3354051/93‑1 (MQ=255)
                                                                            aTCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGCCAAAACCGTCGACGTCGCCGCATCGCt               >  1:3569756/1‑93 (MQ=255)
                                                                            aTCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGCCAAAACCGTCGACGTCGCCGCATCGCt               >  2:7119810/1‑93 (MQ=255)
                                                                             tCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTc              <  1:3258233/93‑1 (MQ=255)
                                                                             tCTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGCCAAAACCGTCGACATCGCCGCATCGCTc              >  2:3466500/1‑93 (MQ=255)
                                                                              cTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGCCAAAACCGTCGACGTCGCCGCATCGCTCa             <  2:112796/93‑1 (MQ=255)
                                                                              cTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGCCAAAACCGTCGACGTCGCCGCATCGCTCa             >  1:2219630/1‑93 (MQ=255)
                                                                              cTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGCCAAAACCGTCGACGTCGCCGCATCGCTCa             <  1:1889272/93‑1 (MQ=255)
                                                                              cTTTTATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGCCAAAACCGTCGACGTCGCCGCATCGCTCa             <  1:3171735/93‑1 (MQ=255)
                                                                               ttttATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTCAc            >  1:6070518/1‑93 (MQ=255)
                                                                               ttttATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGCCAAAACCGTCGACGTCGCCGCATCGCTCAc            <  1:3525909/93‑1 (MQ=255)
                                                                               ttttATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGCCAAAACCGTCGACGTCGCCGCATCGCTCAc            >  2:5641200/1‑93 (MQ=255)
                                                                                  tATTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGCCAAAACCGTCGACGTCGCCGCATCGCTCAcgcg         >  2:6665208/1‑93 (MQ=255)
                                                                                   aTTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGt        <  1:6430253/93‑1 (MQ=255)
                                                                                   aTTCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAAACCGCAGATTTTCAGCGCCAAAACCGTCGACGTCGCCGCATCGCTCACGCGt        >  2:6544186/1‑93 (MQ=255)
                                                                                    ttCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTa       >  1:611511/1‑93 (MQ=255)
                                                                                     tCGCGACGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGCCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTaa      >  2:4393902/1‑93 (MQ=255)
                                                                                         gATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTAc  <  2:2688420/93‑1 (MQ=255)
                                                                                           |                                                                                          
TGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTAC  >  NC_000913/3867095‑3867276

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