Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 391,703 (G)6→7 pseudogene (1453/1458 nt) yaiT → pseudogene, autotransporter family,putative structure, Not classified, interrupted by IS3, putative flagellin structural protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913391,6971.G100.0% 258.5 / NA 76pseudogene (1447/1458 nt)yaiTpseudogene, autotransporter family,putative structure, Not classified, interrupted by IS3, putative flagellin structural protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (41/35);  total (41/35)

GAATCTAAATTAACTTCAGATACTATCGATTTAACTAACGGTTCTTCGCTGAACATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAAC  >  NC_000913/391562‑391835
                                                                                                                                        |                                                                                                                                          
gAATCTAAATTAACTTCAGATACTATCGATTTAACTAACGGTTCTTCGCTGAACATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACtgggtggg                                                                                                                                         <  1:522346/139‑1 (MQ=255)
    cTAAATTAACTTCAGATACTATCGATTTAACTAACGGTTCTTCGCTGAACATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTggggggg                                                                                                                                     <  1:148912/139‑2 (MQ=255)
         ttAACTTCAGATACTATCGATTTAACTAACGGTTCTTCGCTGAACATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCtgat                                                                                                                                 <  2:2442/139‑1 (MQ=255)
         ttAACTTCAGATACTATCGATTTAACTAACGGTTCTTCGCTGAACATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCtgat                                                                                                                                 <  2:360363/139‑1 (MQ=255)
                  gATACTATCGATTTAACTAACGGTTCTTCGCTGAACATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTAc                                                                                                                        <  2:114126/139‑1 (MQ=255)
                   aTACTATCGATTTAACTAACGGTTCTTCGCTGAACATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACATTACTGTTGCCTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTAcc                                                                                                                       <  1:904985/139‑1 (MQ=255)
                       tATCGATTTAACTAACGGTTCTTCGCTGAACATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACg                                                                                                                   >  2:1425907/1‑139 (MQ=255)
                        aTCGATTTAACTAACGGTTCTTCGCTGAACATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGt                                                                                                                  >  2:324605/1‑139 (MQ=255)
                          cGATTTAACTAACGGTTCTTCGCTGAACATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTaa                                                                                                                >  2:1009332/1‑139 (MQ=255)
                             tttAACTAACGGTTCTTCGCTGAACATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTgggg                                                                                                                                        >  2:446384/1‑111 (MQ=255)
                             tttAACTAACGGTTCTTCGCTGAACATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTgggg                                                                                                                                        <  1:446384/111‑1 (MQ=255)
                                aactaacGGTTCTTCGCTGAACATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGa                                                                                                          >  1:793283/1‑139 (MQ=255)
                                 actaacGGTTCTTCGCTGAACATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTTACCGAATCTTCTGGGTGGGGGGGCTGATTTATCCTACCCCCGTAATATGGac                                                                                                         >  1:465422/1‑139 (MQ=255)
                                   taacGGTTCTTCGCTGAACATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGacac                                                                                                       <  2:653857/139‑1 (MQ=255)
                                   taacGGTTCTTCGCTGAACATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGacac                                                                                                       <  1:324605/139‑1 (MQ=255)
                                      cGGTTCTTCGCTGAACATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACtgggtg                                                                                                                                           <  2:309338/99‑1 (MQ=255)
                                      cGGTTCTTCGCTGAACATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACtgggtg                                                                                                                                           >  1:309338/1‑99 (MQ=255)
                                       ggTTCTTCGCTGAACATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGc                                                                                                   >  2:202743/1‑139 (MQ=255)
                                        gTTCTTCGCTGAACATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGct                                                                                                  <  2:793283/139‑2 (MQ=255)
                                           cttcGCTGAACATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTaa                                                                                               >  2:152447/1‑139 (MQ=255)
                                                   aaCATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGtt                                                                                       >  2:241828/1‑139 (MQ=255)
                                                             aaGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCa                                                                             >  1:1411145/1‑139 (MQ=255)
                                                               gATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCaaa                                                                           >  2:258345/1‑139 (MQ=255)
                                                               gATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCaaa                                                                           <  1:202743/139‑1 (MQ=255)
                                                                            aTACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTAcc                                                                                                                       >  2:328117/1‑82 (MQ=255)
                                                                            aTACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTAcc                                                                                                                       <  1:328117/82‑1 (MQ=255)
                                                                            aTACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTACTATGGCAACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTg                                                              >  1:1012728/1‑139 (MQ=25)
                                                                              aCCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAg                                                            >  2:3205/1‑139 (MQ=255)
                                                                              aCCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAg                                                            >  2:717691/1‑139 (MQ=255)
                                                                                 gATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGgcc                                                         <  1:1061859/139‑1 (MQ=255)
                                                                                    catcTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGacac                                                                                                       >  2:1403943/1‑90 (MQ=255)
                                                                                    catcTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGacac                                                                                                       <  1:1403943/90‑1 (MQ=255)
                                                                                          aCTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACtgggt                                                                                                                                             >  2:308040/1‑46 (MQ=255)
                                                                                          aCTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACtgggt                                                                                                                                             <  1:308040/46‑1 (MQ=255)
                                                                                              ttAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGGTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACtg                                            >  1:302826/1‑139 (MQ=255)
                                                                                              ttAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACtg                                            >  1:53654/1‑139 (MQ=255)
                                                                                              ttAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACtg                                            >  1:806829/1‑139 (MQ=255)
                                                                                              ttAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACtg                                            >  2:729574/1‑139 (MQ=255)
                                                                                               tAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACtgt                                           >  1:274236/1‑139 (MQ=255)
                                                                                                aaCTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCtgat                                                                                                                                 >  2:1129701/1‑52 (MQ=255)
                                                                                                aaCTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCtgat                                                                                                                                 <  1:1129701/52‑1 (MQ=255)
                                                                                                 aCTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGt                                                                       <  1:359891/109‑1 (MQ=255)
                                                                                                 aCTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGt                                                                       >  2:359891/1‑109 (MQ=255)
                                                                                                    ccTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAg                                                            >  2:1047182/1‑117 (MQ=255)
                                                                                                    ccTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAg                                                            <  1:1047182/117‑1 (MQ=255)
                                                                                                     cTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTgccgcc                                     <  1:729574/139‑1 (MQ=255)
                                                                                                         aGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCa                                 <  1:258345/139‑1 (MQ=255)
                                                                                                          gTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCCTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCAc                                <  2:410437/139‑1 (MQ=255)
                                                                                                           tACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGgccgc                                                       >  1:1179363/1‑115 (MQ=255)
                                                                                                           tACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGgccgc                                                       <  2:1179363/115‑1 (MQ=255)
                                                                                                             cTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTAccc                                                                                                                      <  1:86398/50‑1 (MQ=255)
                                                                                                             cTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTAccc                                                                                                                      >  2:86398/1‑50 (MQ=255)
                                                                                                             cTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACg                                                                                                                   >  1:1279050/1‑53 (MQ=39)
                                                                                                             cTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACg                                                                                                                   <  2:1279050/53‑1 (MQ=39)
                                                                                                             cTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTGCGGACTGAGGCCGCCACACCAACTGTgccg                                       >  2:890428/1‑129 (MQ=255)
                                                                                                             cTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTgccg                                       <  1:890428/129‑1 (MQ=255)
                                                                                                             cTGTTGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGa                             >  2:528366/1‑139 (MQ=255)
                                                                                                              ggTTGCGTTGACCGAAGCGACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGGGACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAATTGTGCCGCCGCCACCGAt                            <  1:15653/138‑1 (MQ=17)
                                                                                                               gTTGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGt                                                                                                                  >  2:446526/1‑52 (MQ=39)
                                                                                                               gTTGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGt                                                                                                                  <  1:446526/52‑1 (MQ=39)
                                                                                                                ttGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGacac                                                                                                       >  1:597250/1‑62 (MQ=255)
                                                                                                                ttGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGacac                                                                                                       <  2:597250/62‑1 (MQ=255)
                                                                                                                ttGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTg                          >  2:1452067/1‑139 (MQ=255)
                                                                                                                ttGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTg                          >  1:591691/1‑139 (MQ=255)
                                                                                                                ttGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTg                          >  1:320404/1‑139 (MQ=255)
                                                                                                                ttGCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTg                          >  1:1380815/1‑139 (MQ=255)
                                                                                                                  gCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCaaa                                                                           >  2:1026544/1‑88 (MQ=18)
                                                                                                                  gCGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCaaa                                                                           <  1:1026544/88‑1 (MQ=18)
                                                                                                                   cGTTGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTaa                       >  2:711777/1‑139 (MQ=255)
                                                                                                                     ttGACCGAATCTACTGGGTGGGGGGGCtgattgat                                                                                                                             <  2:434843/35‑1 (MQ=255)
                                                                                                                     ttGACCGAATCTACTGGGTGGGGGGGCtgattgat                                                                                                                             >  1:434843/1‑35 (MQ=255)
                                                                                                                     ttGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACTCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATc                     >  2:809240/1‑139 (MQ=255)
                                                                                                                     ttGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCaaa                                                                           >  2:1348103/1‑85 (MQ=18)
                                                                                                                     ttGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCaaa                                                                           <  1:1348103/85‑1 (MQ=18)
                                                                                                                     ttGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGgccgc                                                       >  2:1042594/1‑105 (MQ=255)
                                                                                                                     ttGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATc                     <  1:241828/139‑1 (MQ=255)
                                                                                                                     ttGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGACGCCACACCAACTGTGCCGCCGCCACCGATTGTAAt                      <  1:153045/138‑1 (MQ=255)
                                                                                                                     ttGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGACGCCACACCAACTGTGCCGCCGCCACCGATTGTAAt                      >  2:153045/1‑138 (MQ=255)
                                                                                                                     ttGACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGATTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGgccgc                                                       <  1:1042594/105‑1 (MQ=17)
                                                                                                                       gACCGAATCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGGGCCGCCGCCACCGATTGTAATcac                   >  2:821366/1‑139 (MQ=255)
                                                                                                                             aTCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGa             <  2:53654/139‑1 (MQ=255)
                                                                                                                             aTCTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGa             <  2:806829/139‑1 (MQ=255)
                                                                                                                               cTACTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGata           <  1:717691/139‑1 (MQ=255)
                                                                                                                                 aCTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAAGTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGAtata         <  1:1271912/139‑1 (MQ=255)
                                                                                                                                  cTGGGTGGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATaa        <  2:1411145/139‑1 (MQ=255)
                                                                                                                                      gtgggGGGCTGATTGATCCTACCCACGTAATATGGACAAAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAac  >  2:292888/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                                          
GAATCTAAATTAACTTCAGATACTATCGATTTAACTAACGGTTCTTCGCTGAACATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAAC  >  NC_000913/391562‑391835

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