Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,362,693 T→C L751P (CTC→CCC)  yhcD → putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,362,6930TC99.8% 1290.3 / ‑5.3 404L751P (CTC→CCC) yhcDputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (197/205);  total (198/206)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.59e-01

Only 100 of 408 total aligned reads displayed.

TCCGACGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCTCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCAGAAGATGCACCAA  >  NC_000913/3362602‑3362782
                                                                                           |                                                                                         
 ccGACGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCa                                                                                         >  1:3245750/1‑93 (MQ=255)
   gACGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAAt                                                                                       <  2:1619338/93‑1 (MQ=255)
    aCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATa                                                                                      <  1:5190703/93‑1 (MQ=255)
    aCGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATa                                                                                      >  2:5269618/1‑93 (MQ=255)
       gCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACa                                                                                   >  1:584733/1‑93 (MQ=255)
       gCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACa                                                                                   >  2:1529068/1‑93 (MQ=255)
         aaaaaaCTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATc                                                                                 <  1:2581610/93‑1 (MQ=255)
          aaaaaCTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCa                                                                                >  2:3230045/1‑93 (MQ=255)
             aaCTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAAc                                                                             <  2:1163760/93‑1 (MQ=255)
               cTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCATCTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCg                                                                           <  2:5907669/93‑1 (MQ=255)
                tGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCgg                                                                          <  1:5102919/93‑1 (MQ=255)
                 gCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGa                                                                         <  2:7088316/93‑1 (MQ=255)
                   ccAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACa                                                                       <  1:1429079/93‑1 (MQ=255)
                   ccAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACa                                                                       >  1:1358613/1‑93 (MQ=255)
                   ccAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACa                                                                       <  1:4040822/93‑1 (MQ=255)
                    cAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACAt                                                                      >  2:2890067/1‑93 (MQ=255)
                    cAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACAt                                                                      >  1:1869051/1‑93 (MQ=255)
                     agagGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATg                                                                     <  1:27318/93‑1 (MQ=255)
                      gagGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATgg                                                                    <  2:5170774/93‑1 (MQ=255)
                      gagGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATgg                                                                    <  1:93826/93‑1 (MQ=255)
                      gagGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATgg                                                                    >  1:4134689/1‑93 (MQ=255)
                      gagGCGTTGCGATTGAAGATAGTGAAGGCAACGATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATgg                                                                    >  2:6584926/1‑93 (MQ=255)
                        ggCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTg                                                                  >  2:5359967/1‑93 (MQ=255)
                         gCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGc                                                                 <  1:1095794/93‑1 (MQ=255)
                           gTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTa                                                               <  2:4303882/93‑1 (MQ=255)
                           gTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTa                                                               <  1:2451595/93‑1 (MQ=255)
                           gTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTa                                                               <  1:4807032/93‑1 (MQ=255)
                           gTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTa                                                               <  1:5293711/93‑1 (MQ=255)
                           gTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTa                                                               >  1:7115476/1‑93 (MQ=255)
                           gTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTa                                                               >  1:5720649/1‑93 (MQ=255)
                            ttGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAg                                                              <  1:1641311/93‑1 (MQ=255)
                            ttGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAg                                                              >  2:3821147/1‑93 (MQ=255)
                             tGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGa                                                             <  1:1785562/93‑1 (MQ=255)
                             tGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGa                                                             >  2:7546139/1‑93 (MQ=255)
                              gCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGat                                                            <  1:2298739/93‑1 (MQ=255)
                              gCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGat                                                            <  2:2629983/93‑1 (MQ=255)
                                gATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGatat                                                          >  1:4304775/1‑93 (MQ=255)
                                 aTTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGAtata                                                         <  1:6410928/93‑1 (MQ=255)
                                 aTTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGAtata                                                         >  1:7837212/1‑93 (MQ=255)
                                  ttGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATaa                                                        >  1:4534513/1‑93 (MQ=255)
                                   tGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATaaa                                                       <  1:7474780/93‑1 (MQ=255)
                                    gAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATaaaa                                                      >  2:2203292/1‑93 (MQ=255)
                                    gAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATaaaa                                                      >  1:104339/1‑93 (MQ=255)
                                     aaGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAg                                                     >  2:4230085/1‑93 (MQ=255)
                                      aGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGa                                                    <  1:6087757/93‑1 (MQ=255)
                                       gATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGAt                                                   <  2:7451674/93‑1 (MQ=255)
                                       gATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGAt                                                   <  1:6486404/93‑1 (MQ=255)
                                       gATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGAt                                                   >  2:6717987/1‑93 (MQ=255)
                                        aTAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATg                                                  >  2:1336032/1‑93 (MQ=255)
                                         tAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGa                                                 <  2:1741301/93‑1 (MQ=255)
                                         tAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGa                                                 >  1:4185988/1‑93 (MQ=255)
                                         tAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGa                                                 >  1:3943321/1‑93 (MQ=255)
                                         tAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGa                                                 <  1:7339854/93‑1 (MQ=255)
                                           gTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGagc                                               >  2:7754809/1‑93 (MQ=255)
                                            tGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGagca                                              <  1:7727532/93‑1 (MQ=255)
                                             gAAGGCAACTATCTGACAAGATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGagcag                                             <  2:5888741/93‑1 (MQ=255)
                                             gAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGagcag                                             <  1:7620673/93‑1 (MQ=255)
                                                 gCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAAt                                         >  2:4837199/1‑93 (MQ=255)
                                                   aaCTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCa                                       >  1:2010345/1‑93 (MQ=255)
                                                    aCTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCAt                                      >  2:1438048/1‑93 (MQ=255)
                                                      tATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGc                                    <  2:6950931/93‑1 (MQ=255)
                                                      tATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGc                                    >  1:6798782/1‑93 (MQ=255)
                                                      tATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGc                                    >  1:2175935/1‑93 (MQ=255)
                                                      tATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGc                                    <  2:4142165/93‑1 (MQ=255)
                                                       aTCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGcc                                   <  2:736455/93‑1 (MQ=255)
                                                       aTCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGcc                                   >  2:3926383/1‑93 (MQ=255)
                                                          tGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTa                                >  2:4992940/1‑93 (MQ=255)
                                                           gACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTAt                               <  2:5093102/93‑1 (MQ=255)
                                                           gACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTAt                               >  1:3632662/1‑93 (MQ=255)
                                                            acaacaTCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTAtt                              >  2:6250998/1‑93 (MQ=255)
                                                             caacaTCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTc                             <  1:2843425/93‑1 (MQ=255)
                                                                caTCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCAcc                          >  2:3491826/1‑93 (MQ=255)
                                                                caTCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCAcc                          <  1:4519399/93‑1 (MQ=255)
                                                                caTCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCAcc                          <  1:6357549/93‑1 (MQ=255)
                                                                caTCAGTGGATGACGGTGTTGTATTCCCCAATAACAGCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCAcc                          <  2:7935427/93‑1 (MQ=255)
                                                                 aTCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCt                         <  1:3454954/93‑1 (MQ=255)
                                                                   cAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTa                       >  2:3127945/1‑93 (MQ=255)
                                                                   cAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAAACGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTa                       >  2:1981368/1‑93 (MQ=255)
                                                                    aGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTAc                      >  1:6222143/1‑93 (MQ=255)
                                                                    aGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTAc                      >  2:5851456/1‑93 (MQ=255)
                                                                     gTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACa                     <  2:6882679/93‑1 (MQ=255)
                                                                     gTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACa                     <  1:5796586/93‑1 (MQ=255)
                                                                     gTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACa                     <  2:6496364/93‑1 (MQ=255)
                                                                      tGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACAt                    <  2:5716351/93‑1 (MQ=255)
                                                                         aTGACGGTGTTGTATTCCCCAATAACATCAAACCGGACGTGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTcc                 >  2:5342097/1‑93 (MQ=255)
                                                                         aTGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTcc                 >  2:4510488/1‑93 (MQ=255)
                                                                          tGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTccc                >  2:2743350/1‑93 (MQ=255)
                                                                           gACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCa               <  1:5998132/93‑1 (MQ=255)
                                                                               gtgtTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCagaag           <  2:6081179/93‑1 (MQ=255)
                                                                               gtgtTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCagaag           >  1:1532642/1‑93 (MQ=255)
                                                                                tgttgtATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCagaaga          >  2:4631244/1‑93 (MQ=255)
                                                                                tgttgtATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCagaaga          >  2:2472357/1‑93 (MQ=255)
                                                                                tgttgtATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCagaaga          >  2:680051/1‑93 (MQ=255)
                                                                                  ttgtATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCAGAAGATg        >  2:3684586/1‑93 (MQ=255)
                                                                                   tgtATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCAGAAGATGc       <  1:817075/93‑1 (MQ=255)
                                                                                    gtATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCAGAAGATGCa      >  1:2922533/1‑93 (MQ=255)
                                                                                    gtATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCAGAAGATGCa      <  2:1084686/93‑1 (MQ=255)
                                                                                     tATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCTTATTCACCTTACATTCCCACAAGATGCAc     >  1:1157571/1‑93 (MQ=255)
                                                                                       ttCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCAGAAGATGCACCa   >  1:6945995/1‑93 (MQ=255)
                                                                                        tCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCAGAAGATGCACCaa  >  1:5513485/1‑93 (MQ=255)
                                                                                           |                                                                                         
TCCGACGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCTCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCAGAAGATGCACCAA  >  NC_000913/3362602‑3362782

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