Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,325,494 C→G A261G (GCG→GGG)  yjdA → putative vimentin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,325,4940CG99.6% 922.1 / ‑5.5 267A261G (GCG→GGG) yjdAputative vimentin
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (133/133);  minor base T (0/1);  total (133/134)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 290 total aligned reads displayed.

AGTTGACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAATCGGTGC  >  NC_000913/4325404‑4325583
                                                                                          |                                                                                         
aGTTGACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTcggggg                                                                                         >  1:448993/1‑93 (MQ=255)
 gTTGACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGt                                                                                        >  2:6099829/1‑93 (MQ=255)
 gTTGACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGt                                                                                        >  2:5926007/1‑93 (MQ=255)
   tGACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTAc                                                                                      <  1:7412620/93‑1 (MQ=255)
   tGACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTAc                                                                                      <  2:8257342/93‑1 (MQ=255)
   tGACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTAc                                                                                      <  2:1655249/93‑1 (MQ=255)
    gACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACt                                                                                     >  2:4029350/1‑93 (MQ=255)
    gACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACt                                                                                     <  1:1004356/93‑1 (MQ=255)
     aCGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTg                                                                                    <  2:5112558/93‑1 (MQ=255)
     aCGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTg                                                                                    >  2:1544213/1‑93 (MQ=255)
     aCGTGACTGGATACCCCCGGGCCAAATGAAGCCGTGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTg                                                                                    <  1:8130579/93‑1 (MQ=255)
       gTTACTGGATACCCCCGGGCCAACTGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGc                                                                                  <  1:7523031/93‑1 (MQ=255)
        ttACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCg                                                                                 <  1:1277299/93‑1 (MQ=255)
        ttACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCg                                                                                 <  1:4271035/93‑1 (MQ=255)
         tACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCgg                                                                                <  1:3124087/93‑1 (MQ=255)
          aCTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGt                                                                               <  1:5679202/93‑1 (MQ=255)
          aCTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGt                                                                               <  2:8226305/93‑1 (MQ=255)
           cTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTg                                                                              <  2:2742294/93‑1 (MQ=255)
           cTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTg                                                                              <  1:1104878/93‑1 (MQ=255)
            tGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGc                                                                             >  1:1981280/1‑93 (MQ=255)
             ggATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCt                                                                            <  2:3619157/93‑1 (MQ=255)
              gATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTg                                                                           >  1:3658891/1‑93 (MQ=255)
                   ccccGAGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATta                                                                      <  1:789096/93‑1 (MQ=255)
                    cccGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATtat                                                                     >  1:842735/1‑93 (MQ=255)
                    cccGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATtat                                                                     >  1:3398017/1‑93 (MQ=255)
                     ccGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATtata                                                                    <  1:225491/93‑1 (MQ=255)
                      cGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATAc                                                                   >  1:5858072/1‑93 (MQ=255)
                        ggCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGc                                                                 <  1:2648173/93‑1 (MQ=255)
                         gCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCa                                                                >  1:8005553/1‑93 (MQ=255)
                          caaaaTGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCaa                                                               <  1:6406678/91‑1 (MQ=255)
                           cAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAAc                                                              <  2:6195252/93‑1 (MQ=255)
                            aaaTGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACt                                                             <  1:4794095/93‑1 (MQ=255)
                            aaaTGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACCCAACt                                                             >  1:6741984/1‑93 (MQ=255)
                               tGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATCCGCAACTGaa                                                          >  2:1438689/1‑93 (MQ=255)
                               tGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGaa                                                          <  1:6247871/93‑1 (MQ=255)
                                 aaGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAAt                                                        <  2:5383706/93‑1 (MQ=255)
                                  aGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATc                                                       <  2:7752685/93‑1 (MQ=255)
                                  aGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATc                                                       <  2:7384989/93‑1 (MQ=255)
                                   gCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCg                                                      >  2:6417162/1‑93 (MQ=255)
                                      gggCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGAtc                                                   <  1:6875678/93‑1 (MQ=255)
                                      gggCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGAtc                                                   >  2:6003724/1‑93 (MQ=255)
                                      gggCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGAtc                                                   <  1:973282/93‑1 (MQ=255)
                                      gggCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGAtc                                                   >  1:4021960/1‑93 (MQ=255)
                                       ggCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGAtct                                                  >  1:1845567/1‑93 (MQ=255)
                                       ggCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGAtct                                                  >  2:7210936/1‑93 (MQ=255)
                                       ggCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGAtct                                                  >  2:8106104/1‑93 (MQ=255)
                                        gCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGAtctc                                                 <  2:2077507/93‑1 (MQ=255)
                                         cAACCGCGTCTGCAAAAAATGCTGAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTcc                                                <  1:8629578/93‑1 (MQ=255)
                                         cAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTcc                                                >  2:4553120/1‑93 (MQ=255)
                                         cAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTcc                                                <  1:7194576/93‑1 (MQ=255)
                                         cAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTcc                                                <  2:2334343/93‑1 (MQ=255)
                                           aCCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGa                                              <  1:5332797/93‑1 (MQ=255)
                                             cGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATg                                            <  2:69082/93‑1 (MQ=255)
                                              gCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGa                                           >  2:3079138/1‑93 (MQ=255)
                                               cATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGCTGaa                                          >  1:2704936/1‑93 (MQ=255)
                                               cATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGaa                                          >  2:4363846/1‑93 (MQ=255)
                                               cATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGaa                                          <  1:8303115/93‑1 (MQ=255)
                                                aTCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAag                                         <  2:8752096/93‑1 (MQ=255)
                                                aTCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGAACTTGCGGTGCTCGATTATACGCAACTGAAATCGATGCCCGATAAag                                         >  2:7092738/1‑93 (MQ=255)
                                                 tCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAaga                                        <  1:5575400/93‑1 (MQ=255)
                                                 tCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAaga                                        >  1:1639180/1‑93 (MQ=255)
                                                 tCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAaga                                        <  2:5568889/93‑1 (MQ=255)
                                                   tGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGCAATCGATCTCAGATGAAGAgg                                      >  2:2147958/1‑93 (MQ=255)
                                                   tGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAgg                                      >  2:2808388/1‑93 (MQ=255)
                                                   tGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAgg                                      <  1:3482549/93‑1 (MQ=255)
                                                      aaaaaaTGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTcc                                   >  1:5058070/1‑93 (MQ=255)
                                                      aaaaaaTGCTTAACCAGCAGCTGGCACGAGGCTCGGGGGTACTGGCGGGGCTGGATTATACGCAACTGAAATCGCTCTCCGATGAAGAGGTcc                                   >  2:8459125/1‑93 (MQ=255)
                                                      aaaaaaTGCTTAACCAGCAGCGGGCACGCGCCTCGGGGGTACTGGCGGGGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTcc                                   <  1:376286/93‑1 (MQ=255)
                                                           aTGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAg                              >  2:1382640/1‑93 (MQ=255)
                                                             gCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGc                            >  2:5288264/1‑93 (MQ=255)
                                                              cTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCg                           <  2:4229267/93‑1 (MQ=255)
                                                                 aaCCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGAtt                        >  2:1782788/1‑93 (MQ=255)
                                                                    cagcagCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTg                     <  1:8726860/93‑1 (MQ=255)
                                                                    cagcagCCGGCACGCGCCTCGGGGGTACTGGCGGTGCTTGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTg                     <  1:6595539/93‑1 (MQ=255)
                                                                     agcagcTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTgg                    <  1:7330799/93‑1 (MQ=255)
                                                                     agcagcTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTgg                    <  1:4117282/93‑1 (MQ=255)
                                                                      gcagcTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGc                   <  2:3867895/93‑1 (MQ=255)
                                                                       cagcTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCg                  >  1:5355592/1‑93 (MQ=255)
                                                                        agcTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCgg                 >  2:6734076/1‑93 (MQ=255)
                                                                        agcTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCgg                 >  1:6787890/1‑93 (MQ=255)
                                                                        agcTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCgg                 >  2:3435627/1‑93 (MQ=255)
                                                                         gcTGGCCCGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGGAATCGATCTCCGATTACGAGGTCCGTGAGGGGATTTTCGCggg                >  1:5826837/1‑92 (MQ=255)
                                                                         gcTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCggg                >  1:780445/1‑92 (MQ=255)
                                                                         gcTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCggg                >  2:6774725/1‑92 (MQ=255)
                                                                          gtGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTg               <  1:7650204/92‑1 (MQ=255)
                                                                          cTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTg               <  1:1577871/93‑1 (MQ=255)
                                                                          cTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTACGTGAGGCGATTTTGGCGGTg               <  1:1516685/93‑1 (MQ=255)
                                                                          cTGGCACGAGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATTTACGATGAATAGGTCCGTGAGGCGATTTTGGCGGTg               >  2:5710179/1‑93 (MQ=255)
                                                                               aCGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCa          >  1:6856472/1‑93 (MQ=255)
                                                                                 gcgcCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAAt        <  1:5174667/93‑1 (MQ=255)
                                                                                 gcgcCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAAt        <  1:4701388/93‑1 (MQ=255)
                                                                                 gcgcCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAAt        <  2:7948015/93‑1 (MQ=255)
                                                                                  cgcCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAATc       <  2:3181406/93‑1 (MQ=255)
                                                                                  cgcCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAATc       <  1:2390940/93‑1 (MQ=255)
                                                                                    ccTCGGGGTTGCTGGCGGTGCTGGATTCTACGCAATTGAAATCCATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAATCgg     <  1:2991168/93‑1 (MQ=255)
                                                                                    ccTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAATCgg     <  1:7087826/93‑1 (MQ=255)
                                                                                     cTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAATCGGt    >  2:4904350/1‑93 (MQ=255)
                                                                                       cgggggTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAATCGGTGc  <  1:933196/93‑1 (MQ=255)
                                                                                       cgggggTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAATCGGTGc  <  2:575137/93‑1 (MQ=255)
                                                                                       cgggggTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAATCGGTGc  <  2:4967447/93‑1 (MQ=255)
                                                                                          |                                                                                         
AGTTGACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAATCGGTGC  >  NC_000913/4325404‑4325583

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