Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,253,781 T→C E171G (GAG→GGG)  yhaNM ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,253,7810TC99.7% 1011.9 / ‑4.9 316E171G (GAG→GGG) yhaNMhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (158/157);  minor base A (1/0);  total (159/157)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01

Only 100 of 335 total aligned reads displayed.

ATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACAC  >  NC_000913/3253690‑3253870
                                                                                            |                                                                                         
aTCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCccct                                                                                         <  1:3260315/93‑1 (MQ=255)
aTCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGCAGACTCTTGCTCCCCccct                                                                                         >  1:1076719/1‑93 (MQ=255)
  cgccgcAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGc                                                                                       >  1:7601714/1‑93 (MQ=255)
   gccgcAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAATGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGcc                                                                                      >  1:6084229/1‑93 (MQ=255)
    ccgcAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCca                                                                                     <  1:2622995/93‑1 (MQ=255)
     cgcAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCcac                                                                                    >  1:6890879/1‑93 (MQ=255)
      gcAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCcaca                                                                                   >  2:6485437/1‑93 (MQ=255)
      gcAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCcaca                                                                                   <  1:2197234/93‑1 (MQ=255)
      gcAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCcaca                                                                                   >  1:2025493/1‑93 (MQ=255)
        aaaCGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCcacaca                                                                                 >  2:3905480/1‑93 (MQ=255)
        aaaCGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCcacaca                                                                                 >  2:5921280/1‑93 (MQ=255)
         aaCGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCacacac                                                                                >  2:4350842/1‑93 (MQ=255)
         aaCGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGCCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGATATCCCCCTGTCgcacac                                                                                >  2:4166290/1‑93 (MQ=255)
           cGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGc                                                                              >  1:6017715/1‑93 (MQ=255)
            gggACTTTGTTGACGAACTTCAGTATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGcc                                                                             <  2:4016365/93‑1 (MQ=255)
            gggACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGcc                                                                             <  1:6943677/93‑1 (MQ=255)
             ggACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCct                                                                            >  2:3942182/1‑93 (MQ=255)
             ggACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCct                                                                            >  2:3317611/1‑93 (MQ=255)
              gACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCctg                                                                           >  2:1246664/1‑93 (MQ=255)
              gACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCctg                                                                           >  2:1937071/1‑93 (MQ=255)
               aCTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCctgc                                                                          <  1:10005629/93‑1 (MQ=255)
                 ttCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCctgctg                                                                        <  1:4493902/93‑1 (MQ=255)
                  tCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTgg                                                                       >  2:8710719/1‑93 (MQ=255)
                  tCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTgg                                                                       <  2:7412908/93‑1 (MQ=255)
                  tCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTgg                                                                       <  1:5247106/93‑1 (MQ=255)
                    aTTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGt                                                                     >  2:1550316/1‑93 (MQ=255)
                     ttGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTa                                                                    <  1:338961/93‑1 (MQ=255)
                       gACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTaaa                                                                  <  2:2663420/93‑1 (MQ=255)
                        aCGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAac                                                                 >  1:705869/1‑93 (MQ=255)
                        aCGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAac                                                                 >  2:1999817/1‑93 (MQ=255)
                        aCGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAac                                                                 <  1:3525815/93‑1 (MQ=255)
                          gAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAacac                                                               <  1:8408671/93‑1 (MQ=255)
                           aaCTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAacacc                                                              >  2:6118994/1‑93 (MQ=255)
                            aCTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAacacca                                                             <  1:1873478/93‑1 (MQ=255)
                            aCTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAacacca                                                             >  1:6703184/1‑93 (MQ=255)
                             cTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAacaccac                                                            <  1:1546598/93‑1 (MQ=255)
                              ttCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAacaccaca                                                           <  1:4973244/93‑1 (MQ=255)
                               tCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAacaccacac                                                          >  1:7357494/1‑93 (MQ=255)
                                 aGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAAcaccacacca                                                        >  2:9537664/1‑93 (MQ=255)
                                 aGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAAcaccacacca                                                        >  1:1662319/1‑93 (MQ=255)
                                  ggATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCCTGCCACACACGCCTGCTGGGTAAAcaccacacca                                                        >  2:8479936/1‑93 (MQ=255)
                                  ggATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCAt                                                       <  2:5841496/93‑1 (MQ=255)
                                  ggATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCAt                                                       >  2:7989584/1‑93 (MQ=255)
                                  ggATCTCAACCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCAt                                                       <  1:8012118/93‑1 (MQ=255)
                                    aTCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCg                                                     <  1:471582/93‑1 (MQ=255)
                                    aTCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCg                                                     <  2:2997492/93‑1 (MQ=255)
                                    aTCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCg                                                     >  2:2924466/1‑93 (MQ=255)
                                     tctcAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCACCTGCCACACACGCCTGCTGGGTCAACACCACAACAtaga                                                    >  1:7836788/1‑90 (MQ=255)
                                      ctcAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTg                                                   >  2:5187336/1‑93 (MQ=255)
                                        cagccagcGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCg                                                 >  2:9217285/1‑93 (MQ=255)
                                          gccagcGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGtc                                               <  1:9273964/93‑1 (MQ=255)
                                            cagcGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGtctc                                             >  1:5109106/1‑93 (MQ=255)
                                             agcGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCg                                            >  2:10024096/1‑93 (MQ=255)
                                             agcGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCg                                            >  2:1457402/1‑93 (MQ=255)
                                             agcGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCg                                            >  2:6368813/1‑93 (MQ=255)
                                             agcGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCg                                            >  1:5450199/1‑93 (MQ=255)
                                              gcGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGa                                           >  1:2520567/1‑93 (MQ=255)
                                              gcGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGa                                           >  2:7749993/1‑93 (MQ=255)
                                                 gtcgTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGata                                         <  1:9729321/90‑1 (MQ=255)
                                                 gtcgTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGata                                         <  1:3808071/90‑1 (MQ=255)
                                                 gtcgTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGata                                         <  2:8401802/90‑1 (MQ=255)
                                                  tcgTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGatat                                        <  2:4943019/91‑1 (MQ=255)
                                                  tcgTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGatat                                        >  1:5626185/3‑93 (MQ=255)
                                                     ttCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGca                                     <  1:8211009/93‑1 (MQ=255)
                                                     ttCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGca                                     >  2:6040680/1‑93 (MQ=255)
                                                       cTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGcaca                                   >  2:612978/1‑93 (MQ=255)
                                                       cTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGcaca                                   <  2:826518/93‑1 (MQ=255)
                                                        tGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACaa                                  >  1:8781779/1‑93 (MQ=255)
                                                        tGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACaa                                  >  2:1670878/1‑93 (MQ=255)
                                                         ggAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAAt                                 >  2:7354948/1‑93 (MQ=255)
                                                            aaGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGtt                              <  1:9528317/93‑1 (MQ=255)
                                                            aaGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGtt                              <  2:907636/93‑1 (MQ=255)
                                                             aGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTg                             <  1:5649093/93‑1 (MQ=255)
                                                              gCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTgg                            <  1:9443885/93‑1 (MQ=255)
                                                               gaCCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGt                           <  2:9330936/92‑1 (MQ=255)
                                                               cACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGt                           <  2:6221905/93‑1 (MQ=255)
                                                                cccGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTa                          <  2:4837742/92‑1 (MQ=255)
                                                                cccGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTa                          <  1:3706473/92‑1 (MQ=255)
                                                                aCCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTa                          <  2:4507681/93‑1 (MQ=255)
                                                                 ccGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTAt                         >  1:8633368/1‑93 (MQ=255)
                                                                  cGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATg                        <  1:4343430/93‑1 (MQ=255)
                                                                    tCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGcc                      <  2:4775002/93‑1 (MQ=255)
                                                                    tCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGcc                      >  2:7124438/1‑93 (MQ=255)
                                                                    tCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGcc                      <  2:3502376/93‑1 (MQ=255)
                                                                    ccAGCGGAGACCCTTTCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGcc                      <  1:2987716/92‑1 (MQ=255)
                                                                     cAGCGGAGCCTCTTCCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGccc                     <  2:7167217/93‑1 (MQ=255)
                                                                     cAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGccc                     <  1:1408459/93‑1 (MQ=255)
                                                                           agcctctTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCcgac               <  1:6399446/90‑1 (MQ=255)
                                                                           agaCTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCcgac               >  2:6808970/1‑93 (MQ=255)
                                                                            gaCTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCcgacg              >  1:4065185/1‑93 (MQ=255)
                                                                             aCTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCcgacga             <  1:6574146/93‑1 (MQ=255)
                                                                              ccctTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAACCACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGAt            <  1:4030639/91‑1 (MQ=255)
                                                                                ctTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATgg          >  1:588982/1‑93 (MQ=255)
                                                                                 ttGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGt         <  2:7601713/93‑1 (MQ=255)
                                                                                  tGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTg        >  1:1098386/1‑93 (MQ=255)
                                                                                  tGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTg        <  2:7479320/93‑1 (MQ=255)
                                                                                  tGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCACCGACGATGGTg        >  2:7084572/1‑93 (MQ=255)
                                                                                      cGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGacac    <  1:3881140/93‑1 (MQ=255)
                                                                                        cccccCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGacaca   <  2:5660601/92‑1 (MQ=255)
                                                                                         ccccCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGacacac  >  2:2204506/1‑93 (MQ=255)
                                                                                            |                                                                                         
ATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACAC  >  NC_000913/3253690‑3253870

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