Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 391,703 (G)6→5 pseudogene (1453/1458 nt) yaiT → putative autotransporter YaiT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913391,6980G.99.5% 1080.5 / ‑1.8 209pseudogene (1448/1458 nt)yaiTputative autotransporter YaiT
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base . (107/101);  minor base A (1/0);  total (108/101)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 233 total aligned reads displayed.

AACGGTTCTTCGCTGAACATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGT  >  NC_000913/391598‑391794
                                                                                                    |                                                                                                
aacGGTTCTTCGCTGAACATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACtgggtg                                                                                                  >  1:586346/1‑101 (MQ=255)
   ggTTCTTCGCTGAACATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTgggg                                                                                               <  1:585891/101‑1 (MQ=255)
    gTTCTTCGCTGAACATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTggggg                                                                                              >  1:586349/1‑101 (MQ=255)
     ttcttcGCTGAACATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTgggggc                                                                                             <  2:586228/101‑2 (MQ=255)
      tcttcGCTGAACATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTgggggct                                                                                            <  1:585834/101‑3 (MQ=255)
      tcttcGCTGAACATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTgggggct                                                                                            >  2:586353/1‑99 (MQ=255)
       cttcGCTGAACATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTgggggctg                                                                                           <  2:586127/101‑4 (MQ=255)
       cttcGCTGAACATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTgggggctg                                                                                           <  2:585867/101‑4 (MQ=255)
        ttcGCTGAACATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTgggggctga                                                                                          >  2:586358/1‑97 (MQ=255)
         tcGCTGAACATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCtgat                                                                                        >  1:586359/1‑101 (MQ=255)
         tcGCTGAACATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCtgat                                                                                        <  1:585924/101‑1 (MQ=255)
             tGAACATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACtgggtg                                                                                                  >  1:586364/1‑88 (MQ=255)
             tGAACATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCtgattgat                                                                                    <  1:585595/101‑1 (MQ=255)
             tGAACATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCtgattgat                                                                                    >  2:586363/1‑101 (MQ=255)
             tGAACATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTACTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCtgattgat                                                                                    <  2:585904/101‑1 (MQ=255)
                aCATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCt                                                                                 >  1:586368/1‑101 (MQ=255)
                  aTTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTAc                                                                               <  2:585878/101‑1 (MQ=255)
                     ggTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCa                                                                            <  2:585841/101‑1 (MQ=255)
                      gTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCAc                                                                           <  2:585827/101‑1 (MQ=255)
                      gTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCAc                                                                           >  1:586372/1‑101 (MQ=255)
                       tGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACg                                                                          >  2:586374/1‑101 (MQ=255)
                         aaGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTa                                                                        >  2:586376/1‑101 (MQ=255)
                          aGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTaa                                                                       >  2:586378/1‑101 (MQ=255)
                           gATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAat                                                                      <  1:586088/101‑1 (MQ=255)
                            aTGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAata                                                                     <  1:585954/101‑1 (MQ=255)
                            aTGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAata                                                                     >  1:586380/1‑101 (MQ=255)
                             tGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAatat                                                                    <  2:585959/101‑1 (MQ=255)
                             tGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAatat                                                                    <  2:585699/101‑1 (MQ=255)
                              ggCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATg                                                                   >  1:586383/1‑101 (MQ=255)
                               gCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATgg                                                                  <  2:586026/101‑1 (MQ=255)
                                cTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACtgggtggg                                                                                                <  1:585957/71‑1 (MQ=255)
                                cTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGa                                                                 >  2:586386/1‑101 (MQ=255)
                                  aCGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGaca                                                               >  2:586389/1‑101 (MQ=255)
                                   cGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGacac                                                              <  1:586028/101‑1 (MQ=255)
                                    gTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGAcaca                                                             <  1:585718/101‑1 (MQ=255)
                                    gTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGAcaca                                                             <  2:585932/101‑1 (MQ=255)
                                    gTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGAcaca                                                             >  1:586391/1‑101 (MQ=255)
                                      tGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAgg                                                           <  2:585923/101‑1 (MQ=255)
                                       gATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGc                                                          <  1:585792/101‑1 (MQ=255)
                                       gATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGc                                                          <  2:586032/101‑1 (MQ=255)
                                       gATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGc                                                          >  2:586392/1‑101 (MQ=255)
                                         tACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCt                                                        >  1:586395/1‑101 (MQ=255)
                                          aCCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTa                                                       <  1:585983/101‑1 (MQ=255)
                                          aCCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTa                                                       <  2:586055/101‑1 (MQ=255)
                                           ccGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTaa                                                      >  1:586399/1‑101 (MQ=255)
                                            cGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAg                                                     >  2:586402/1‑101 (MQ=255)
                                             gATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGc                                                    >  1:586403/1‑101 (MQ=255)
                                             gATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGc                                                    <  1:585944/101‑1 (MQ=255)
                                              atcatcTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCg                                                   >  2:586405/1‑101 (MQ=255)
                                                 atcTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAgg                                                >  2:586407/1‑101 (MQ=255)
                                                    tGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTAGCCACGTAATATGGACACAGGCCTAAGCGAGGTTc                                             >  1:586409/1‑101 (MQ=255)
                                                    tGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAgg                                                >  1:586412/1‑98 (MQ=255)
                                                    tGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTc                                             <  2:585997/101‑1 (MQ=255)
                                                     gACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCt                                            <  1:585479/101‑1 (MQ=255)
                                                      aCTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCtt                                           <  1:586014/101‑1 (MQ=255)
                                                        tATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGt                                         <  2:585967/101‑1 (MQ=255)
                                                         aTTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGtt                                        >  1:586416/1‑101 (MQ=255)
                                                          ttAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGttt                                       >  2:586417/1‑101 (MQ=255)
                                                          ttAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGttt                                       <  1:586067/101‑1 (MQ=255)
                                                           tAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGGCTAAGCGAGGTTCTTGtttt                                      <  2:586063/101‑1 (MQ=255)
                                                           tAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGtttt                                      >  2:586419/1‑101 (MQ=255)
                                                             aCTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCa                                    >  1:586421/1‑101 (MQ=255)
                                                             aCTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCa                                    >  1:586423/1‑101 (MQ=255)
                                                             aCTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCa                                    <  1:586201/101‑1 (MQ=255)
                                                             aCTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCa                                    >  2:586428/1‑101 (MQ=255)
                                                             aCTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCa                                    >  1:586426/1‑101 (MQ=255)
                                                              cTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCaa                                   <  1:585883/101‑1 (MQ=255)
                                                              cTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCaa                                   >  1:586430/1‑101 (MQ=255)
                                                              cTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCaa                                   <  1:586072/101‑1 (MQ=255)
                                                               tCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCaaa                                  >  1:586432/1‑101 (MQ=255)
                                                                ccTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAAt                                 <  2:585915/101‑1 (MQ=255)
                                                                ccTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAAt                                 <  2:586009/101‑1 (MQ=255)
                                                                 cTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAAtt                                <  2:586259/101‑1 (MQ=255)
                                                                  tACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTg                               <  2:586075/101‑1 (MQ=255)
                                                                    cAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGtt                             <  2:586044/101‑1 (MQ=255)
                                                                    cAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGtt                             >  1:586435/1‑101 (MQ=255)
                                                                     aGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTaa                                                      >  1:586437/1‑75 (MQ=255)
                                                                      gTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTcc                           <  2:585951/101‑1 (MQ=255)
                                                                       tACTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCg                          >  1:586440/1‑101 (MQ=255)
                                                                         cTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTc                            >  2:586445/1‑97 (MQ=255)
                                                                         cTGTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGa                        >  1:586444/1‑101 (MQ=255)
                                                                           gTTGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACt                      <  2:586168/101‑1 (MQ=255)
                                                                             tGCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGa                    >  1:586447/1‑101 (MQ=255)
                                                                              gCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAg                   >  2:586450/1‑101 (MQ=255)
                                                                              gCGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAg                   <  2:586001/101‑1 (MQ=255)
                                                                               cGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAAtt                                >  1:586451/1‑87 (MQ=255)
                                                                               cGTTGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAgg                  <  2:586121/101‑1 (MQ=255)
                                                                                 ttGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGgcc                <  2:585658/101‑1 (MQ=255)
                                                                                  tGACCGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGgccg               >  1:586454/1‑101 (MQ=255)
                                                                                     ccGAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCca            >  2:586457/1‑101 (MQ=255)
                                                                                       gAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCcaca          >  1:586460/1‑101 (MQ=255)
                                                                                       gAATCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCcaca          <  1:585992/101‑1 (MQ=255)
                                                                                        aaTCTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCacac         >  2:586461/1‑101 (MQ=255)
                                                                                           cTACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCaa      >  1:586464/1‑101 (MQ=255)
                                                                                            tACTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTATTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAAc     >  2:586465/1‑101 (MQ=255)
                                                                                             aCTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACt    >  2:586470/1‑101 (MQ=255)
                                                                                             aCTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACt    <  2:586094/101‑1 (MQ=255)
                                                                                             aCTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACt    >  1:586469/1‑101 (MQ=255)
                                                                                             aCTGGGTGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACt    <  1:585912/101‑1 (MQ=255)
                                                                                               tgggtgggGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGa                    >  1:586474/1‑83 (MQ=255)
                                                                                                    |                                                                                                
AACGGTTCTTCGCTGAACATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGT  >  NC_000913/391598‑391794

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