Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,091,304 T→C D70G (GAC→GGC)  ycdS ← putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,091,3040TC97.9% 453.6 / ‑0.6 145D70G (GAC→GGC) ycdSputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (80/62);  total (81/64)
Fisher's exact test for biased strand distribution p-value = 5.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.71e-01

Only 100 of 158 total aligned reads displayed.

CCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCTTGCGAGCTTCAATAATCAAT  >  NC_000913/1091214‑1091393
                                                                                          |                                                                                         
ccACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCa                                                                                         <  2:3593393/93‑1 (MQ=255)
ccACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCa                                                                                         >  1:1822091/1‑93 (MQ=255)
ccACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCa                                                                                         <  1:3680178/93‑1 (MQ=255)
 cACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAAccagccag                                                                                        <  1:1171463/93‑1 (MQ=255)
 cACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAAccagccag                                                                                        >  2:1978832/1‑93 (MQ=255)
  aCGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACcagccagc                                                                                       >  1:1543252/1‑93 (MQ=255)
     cgcTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAAt                                                                                    <  1:4714350/93‑1 (MQ=255)
      gcTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAAtt                                                                                   >  2:3512904/1‑93 (MQ=255)
      gcTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAAtt                                                                                   >  2:1126616/1‑93 (MQ=255)
      gcTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAAtt                                                                                   >  1:369158/1‑93 (MQ=255)
      gcTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAAtt                                                                                   <  1:2676909/93‑1 (MQ=255)
       cTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAAttt                                                                                  <  2:116576/93‑1 (MQ=255)
       cTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAAttt                                                                                  <  1:4688190/93‑1 (MQ=255)
        tGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTg                                                                                 >  1:2661628/1‑93 (MQ=255)
        tGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTg                                                                                 >  1:4407295/1‑93 (MQ=255)
        tGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTg                                                                                 >  1:4385370/1‑93 (MQ=255)
        tGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGGTTATCTTGCCCGGCCCATAAGGCAATCTGCAACCCGCCAGCAATTTg                                                                                 >  1:2661627/1‑93 (MQ=255)
         ggTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTgg                                                                                >  2:1012858/1‑93 (MQ=255)
         ggTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTgg                                                                                >  2:704782/1‑93 (MQ=255)
         ggTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTgg                                                                                >  1:3777579/1‑93 (MQ=255)
         ggTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTgg                                                                                <  2:2646696/93‑1 (MQ=255)
          gTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGt                                                                               >  1:2977559/1‑93 (MQ=255)
             aTTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTAt                                                                            <  2:692955/93‑1 (MQ=255)
             aTTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTAt                                                                            <  2:720307/93‑1 (MQ=255)
             aTTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTAt                                                                            <  1:3530635/93‑1 (MQ=255)
                gCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTGGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGc                                                                         <  2:3182677/93‑1 (MQ=255)
                 cTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCt                                                                        >  1:1055067/1‑93 (MQ=255)
                 cTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCt                                                                        >  2:4706396/1‑93 (MQ=255)
                  tgatgaCGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTg                                                                       >  1:3537619/1‑93 (MQ=255)
                   gatgaCGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGa                                                                      >  1:4081803/1‑93 (MQ=255)
                     tgaCGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGt                                                                    <  2:1232898/93‑1 (MQ=255)
                        cGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCt                                                                 <  2:1280084/93‑1 (MQ=255)
                        cGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCt                                                                 <  2:1237157/93‑1 (MQ=255)
                         ggTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTg                                                                >  1:343395/1‑93 (MQ=255)
                         ggTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTg                                                                <  2:3480315/93‑1 (MQ=255)
                         ggTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGGGCTg                                                                >  2:848921/1‑93 (MQ=255)
                          gTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGa                                                               <  2:3444204/93‑1 (MQ=255)
                          gTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGa                                                               >  1:4693212/1‑93 (MQ=255)
                                gTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAg                                                         >  2:3854827/1‑93 (MQ=255)
                                gTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAg                                                         >  2:636878/1‑93 (MQ=255)
                                 ttGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGt                                                        >  2:1063835/1‑93 (MQ=255)
                                 ttGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGt                                                        >  1:3388764/1‑93 (MQ=255)
                                  tGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTg                                                       >  1:1314018/1‑93 (MQ=255)
                                    tAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCa                                                     <  1:2200854/93‑1 (MQ=255)
                                      aaCGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCaaa                                                   >  1:3110982/1‑93 (MQ=255)
                                       aCGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTACAAAc                                                  >  1:3232107/1‑93 (MQ=255)
                                         ggTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCa                                                >  1:666490/1‑93 (MQ=255)
                                         ggTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCa                                                <  1:236497/93‑1 (MQ=255)
                                          gTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCAt                                               <  1:2014825/93‑1 (MQ=255)
                                          gTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCAt                                               <  1:907024/93‑1 (MQ=255)
                                           taataaCCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATg                                              >  2:3724918/1‑93 (MQ=255)
                                            aataaCCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGa                                             >  2:4013839/1‑93 (MQ=255)
                                              taaCCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACa                                           >  2:229030/1‑93 (MQ=255)
                                                 ccTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAg                                        <  2:4094190/93‑1 (MQ=255)
                                                  cTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAgc                                       >  2:1448051/1‑93 (MQ=255)
                                                   tGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAgct                                      >  1:1297234/1‑93 (MQ=255)
                                                   tGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAgct                                      >  1:2137543/1‑93 (MQ=255)
                                                     tttATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAgctgg                                    <  2:2715447/93‑1 (MQ=255)
                                                      ttATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAgctggc                                   <  1:4035935/93‑1 (MQ=255)
                                                       tATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAgctggct                                  <  2:1676681/93‑1 (MQ=255)
                                                         tCTTGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGa                                >  1:4546808/1‑93 (MQ=255)
                                                           ttGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGt                              <  2:387654/93‑1 (MQ=255)
                                                            tGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTa                             >  2:1632740/1‑93 (MQ=255)
                                                            tGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTa                             >  2:171082/1‑93 (MQ=255)
                                                            tGCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTa                             <  1:3087239/93‑1 (MQ=255)
                                                             gCCCGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTAt                            <  1:4482127/93‑1 (MQ=255)
                                                              cccGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTAtt                           <  2:2874773/93‑1 (MQ=255)
                                                               ccGGCCCATAAGGCAATCTGTAACTAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTa                          >  2:3160223/1‑93 (MQ=255)
                                                               ccGGCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTa                          >  1:604062/1‑93 (MQ=255)
                                                                 ggCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTAcc                        <  2:2368605/93‑1 (MQ=255)
                                                                  gCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTAccc                       <  2:4164696/93‑1 (MQ=255)
                                                                  gCCCATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTTCCCt                      <  2:71645/93‑1 (MQ=255)
                                                                   gccATAAGGCAATCTGTACCCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCAGGGCAAATCTGGCTGAGTATTACCCt                      <  1:3114650/92‑1 (MQ=255)
                                                                   cccATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCt                      >  2:4493178/1‑93 (MQ=255)
                                                                     cATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCTTg                    >  1:3020393/1‑93 (MQ=255)
                                                                     cATAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCTTg                    >  2:594540/1‑93 (MQ=255)
                                                                      aTAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCTTGc                   >  1:2420336/1‑93 (MQ=255)
                                                                      aTAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCTTGc                   >  2:1553761/1‑93 (MQ=255)
                                                                       tAAGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCTTGCg                  >  1:4631937/1‑93 (MQ=255)
                                                                        aaGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCTTGCGa                 >  1:252737/1‑93 (MQ=255)
                                                                         aGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCTTGCGAg                <  2:2058166/93‑1 (MQ=255)
                                                                         aGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCTTGCGAg                <  2:4190015/93‑1 (MQ=255)
                                                                         aGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCTTGCGAg                >  2:4213993/1‑93 (MQ=255)
                                                                         aGGCAATCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCTTGCGAg                >  2:4317244/1‑93 (MQ=255)
                                                                             aaTCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCTTGCGAGCtta            >  2:1916647/1‑92 (MQ=255)
                                                                             aaTCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCTTGCGAGCTTc            <  2:946623/93‑1 (MQ=255)
                                                                             aaTCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCTTGCGAGCTTc            <  1:1898262/93‑1 (MQ=255)
                                                                             aaTCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCTTGCGAGCTTc            >  2:4365564/1‑93 (MQ=255)
                                                                             aaTCTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCTTGCGAGCTTc            >  2:4327419/1‑93 (MQ=255)
                                                                                cTGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCTTGCGAGCTTCaat         <  2:870648/93‑1 (MQ=255)
                                                                                 tGTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCTTGCGAGCTTCaata        >  2:796913/1‑93 (MQ=255)
                                                                                  gTAACCAGCCAGCAATTTGGTTATTGCTTAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCTTGCGAGCTTCaataa       <  1:1934491/93‑1 (MQ=255)
                                                                                  gTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCTTGCGAGCTTCaataa       <  2:4426276/93‑1 (MQ=255)
                                                                                  gTAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCTTGCGAGCTTCaataa       >  2:544213/1‑93 (MQ=255)
                                                                                   tAACCAGCCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCTTGCGAGCTTCaataat      >  1:2850860/1‑93 (MQ=255)
                                                                                      ccagccagCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCTTGCGAGCTTCAATAATCaa   <  1:47809/93‑1 (MQ=255)
                                                                                      ccagccagCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCTTGCGAGCTTCAATAATCaa   <  1:1968072/93‑1 (MQ=255)
                                                                                      ccagccagCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCTTGCGAGCTTCAATAATCaa   <  1:1498977/93‑1 (MQ=255)
                                                                                      ccagccagCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCTTGCGAGCTTCAATAATCaa   <  2:2690861/93‑1 (MQ=255)
                                                                                      ccagccagCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCTTGCGAGCTTCAATAATCaa   <  2:3697055/93‑1 (MQ=255)
                                                                                          |                                                                                         
CCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCTTGCGAGCTTCAATAATCAAT  >  NC_000913/1091214‑1091393

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