Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,022,007 C→T 6.4% G21G (GGC→GGT tatA → TatABCE protein translocation system subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,022,0070CT 6.4% 509.6 / 18.7 156G21G (GGC→GGTtatATatABCE protein translocation system subunit
Reads supporting (aligned to +/- strand):  ref base C (75/71);  new base T (5/5);  total (80/76)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.62e-01

Only 100 of 157 total aligned reads displayed.

GCTGATTTTTTCATCGCTCAAGGCGGGCCGTGTAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTATCGCCGATAAGCAGGCGGATA  >  NC_000913/4021862‑4022155
                                                                                                                                                 |                                                                                                                                                    
     ttttttCATCGCTCAAGGCGGGCCGTGTAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCaaaaa                                                                                                                                              >  2:192082/1‑149 (MQ=255)
         ttCATCGCTCAAGGCGGGCCGTGTAACCTATAATGGGGCTTTGTTTAATCATCCTCTACCACAGAGGAACATTTATGGGTGGTAACAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTc                                                                                                                                          <  1:648598/149‑1 (MQ=255)
            aTCGCTCAAGGCGGGCCGTGTAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAg                                                                                                                                             >  4:546756/1‑143 (MQ=255)
            aTCGCTCAAGGCGGGCCGTGTAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAg                                                                                                                                             <  3:546756/143‑1 (MQ=255)
              cGCTCAAGGCGGGCCGTGTAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGca                                                                                                                                           >  2:898190/1‑142 (MQ=255)
                 tCAAGGCGGGCCGTGTAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCAt                                                                                                                                  >  4:120003/1‑149 (MQ=255)
                        gggCCGTGTAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTcc                                                                                                                           <  1:1310097/149‑1 (MQ=255)
                                tAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGca                                                                                                                                           <  2:1451319/125‑2 (MQ=255)
                                tAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGca                                                                                                                                           >  1:1451319/1‑124 (MQ=255)
                                tAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTgg                                                                                                                   >  2:1029081/1‑149 (MQ=255)
                                       taATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCg                                                                                                            <  4:226990/149‑1 (MQ=255)
                                         aTGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGAt                                                                                                          <  1:1029081/149‑1 (MQ=255)
                                          tGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATc                                                                                                         >  1:643761/1‑149 (MQ=255)
                                               cTTTGTTTAATCATCATCTACCACAGATGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAgg                                                                                                    >  3:975922/1‑149 (MQ=255)
                                                tttGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGc                                                                                                   >  2:976188/1‑149 (MQ=255)
                                                tttGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGc                                                                                                   >  3:977707/1‑149 (MQ=255)
                                                    tttAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATc                                                                                                                       >  1:1062637/1‑125 (MQ=255)
                                                       aaTCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTcc                                                                                                                                    <  1:144938/109‑1 (MQ=255)
                                                       aaTCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTcc                                                                                                                                    >  2:144938/1‑109 (MQ=255)
                                                       aaTCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCtt                                                                                                                     <  2:814553/124‑1 (MQ=255)
                                                       aaTCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCtt                                                                                                                     >  1:814553/1‑124 (MQ=255)
                                                         tcatcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTGTTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGca                                                                                                                                           <  4:850785/100‑2 (MQ=255)
                                                          catcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGc                                                                                                   >  3:36023/1‑139 (MQ=255)
                                                          catcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGc                                                                                                   <  4:36023/139‑1 (MQ=255)
                                                           atcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGCTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAATGc                                                                                        <  4:977707/149‑1 (MQ=255)
                                                             catcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCaa                                                                                      <  1:946838/149‑1 (MQ=255)
                                                                 tACCTCAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCaaa                                                                                                      <  4:1423909/129‑1 (MQ=255)
                                                                 tACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAg                                                                                  >  4:582193/1‑149 (MQ=255)
                                                                 tACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGATCCATCGGTTCCGATCTTGGTGCGTCGATCaaa                                                                                                      >  3:1423909/1‑129 (MQ=255)
                                                                  aCCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGc                                                                                 >  3:1169263/1‑149 (MQ=255)
                                                                       agagGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGa                                                                                   >  1:150036/1‑142 (MQ=255)
                                                                       agagGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGa                                                                                   <  2:150036/142‑1 (MQ=255)
                                                                           gAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTg                                                                                                                    <  1:195654/105‑1 (MQ=255)
                                                                           gAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGaa                                                                        <  3:292417/149‑1 (MQ=255)
                                                                           gAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGaa                                                                        >  4:935419/1‑149 (MQ=255)
                                                                               aTGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTaaa                                                                                             >  3:879885/1‑124 (MQ=255)
                                                                               aTGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCaa                                                                    >  2:1090565/1‑149 (MQ=255)
                                                                               aTGTATGCGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTaaa                                                                                             <  4:879885/124‑1 (MQ=255)
                                                                                  tATGGGTGGTAGCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGc                                                                 >  1:766251/1‑149 (MQ=255)
                                                                                     gggTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGt                                                                                                                  >  1:549077/1‑97 (MQ=255)
                                                                                       gtggtATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGAt                                                            >  3:1122512/1‑149 (MQ=255)
                                                                                          gtATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCaaa                                                                   >  4:1448121/1‑139 (MQ=255)
                                                                                           tATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATaaaa                                                        <  3:120003/149‑1 (MQ=255)
                                                                                            aTCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAgg                                                                                                    <  3:396660/104‑1 (MQ=255)
                                                                                            aTCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAgg                                                                                                    >  4:396660/1‑104 (MQ=255)
                                                                                                gTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGa                                                                                                           <  1:1235227/93‑1 (MQ=255)
                                                                                                gTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGa                                                                                                           >  2:1235227/1‑93 (MQ=255)
                                                                                                 tATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCg                                                                                                                                         >  3:300057/1‑62 (MQ=255)
                                                                                                 tATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCg                                                                                                                                         <  4:300057/62‑1 (MQ=255)
                                                                                                 tATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATaaaa                                                        >  1:898683/1‑143 (MQ=255)
                                                                                                  aTTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGt                                                                                                                  >  3:1217785/1‑84 (MQ=255)
                                                                                                   tttGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAg                                                <  2:643761/149‑1 (MQ=255)
                                                                                                   tttGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGATTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAg                                                >  1:924360/1‑149 (MQ=255)
                                                                                                      ggCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTc                                                  <  3:707479/144‑1 (MQ=255)
                                                                                                      ggCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTc                                                  >  4:707479/1‑144 (MQ=255)
                                                                                                         aGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCt                                          >  3:195480/1‑149 (MQ=255)
                                                                                                         aGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCt                                          >  4:833164/1‑149 (MQ=255)
                                                                                                          gTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCga                                                                               <  1:875277/111‑1 (MQ=255)
                                                                                                          gTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCga                                                                               >  2:875277/1‑111 (MQ=255)
                                                                                                            tATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGAt                                       <  4:1169263/149‑1 (MQ=255)
                                                                                                            tATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGAt                                       <  1:1090565/149‑1 (MQ=255)
                                                                                                              ttGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTcc                                                                                                                           <  2:1045045/63‑1 (MQ=255)
                                                                                                              ttGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGAttt                                     >  3:1367417/1‑149 (MQ=255)
                                                                                                              ttGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGAttt                                     >  3:1405590/1‑149 (MQ=255)
                                                                                                              ttGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGAttt                                     >  1:1357592/1‑149 (MQ=255)
                                                                                                                gATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGaa                                                                        >  3:810613/1‑112 (MQ=255)
                                                                                                                gATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGaa                                                                        <  4:810613/112‑1 (MQ=255)
                                                                                                                gATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTa                                   <  4:195480/149‑1 (MQ=255)
                                                                                                                gATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTa                                   <  4:1367417/149‑1 (MQ=255)
                                                                                                                gATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTa                                   >  2:1177575/1‑149 (MQ=255)
                                                                                                                 attattGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTg                                         <  1:1180038/142‑1 (MQ=255)
                                                                                                                 attattGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTAc                                  <  2:766251/149‑1 (MQ=255)
                                                                                                                  ttattGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACt                                 <  1:1177575/149‑1 (MQ=255)
                                                                                                                     ttGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGa                                                                                   <  4:453766/96‑1 (MQ=255)
                                                                                                                     ttGCCGTCATCGTTGTACTCCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGt                                                                                                                  <  3:841953/65‑1 (MQ=255)
                                                                                                                        ccGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGaaa                           >  2:1156125/1‑149 (MQ=255)
                                                                                                                         cGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTcc                                                                                                                           >  3:930995/1‑52 (MQ=255)
                                                                                                                          gTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGc                                           >  3:759560/1‑131 (MQ=255)
                                                                                                                          gTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAAc                         >  1:1088018/1‑149 (MQ=255)
                                                                                                                           tCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATaaaa                                                        >  4:1256949/1‑117 (MQ=255)
                                                                                                                            cATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTa                       >  1:94053/1‑149 (MQ=255)
                                                                                                                             aTCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGATTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTAt                      <  2:924360/149‑1 (MQ=255)
                                                                                                                              tCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTATc                     >  3:649508/1‑149 (MQ=255)
                                                                                                                                gTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGAtt                                      >  2:1074201/1‑130 (MQ=255)
                                                                                                                                gTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGAtt                                      <  1:1074201/130‑1 (MQ=255)
                                                                                                                                gTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTATCGc                   >  1:190342/1‑149 (MQ=255)
                                                                                                                                 ttGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCgat                                                                              >  3:369299/1‑89 (MQ=255)
                                                                                                                                 ttGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCgat                                                                              <  4:369299/89‑1 (MQ=255)
                                                                                                                                   gTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTATCGCCGa                <  2:1088018/149‑1 (MQ=255)
                                                                                                                                    tACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTATCGCCGAt               >  3:173547/1‑149 (MQ=255)
                                                                                                                                    tACTGCTGTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTAc                                  >  2:938869/1‑130 (MQ=255)
                                                                                                                                     aCTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTATCGCCGATa              >  4:937576/1‑149 (MQ=255)
                                                                                                                                      cTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTATCGCCGATaa             <  2:1357592/149‑1 (MQ=255)
                                                                                                                                       tGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAg                                                                                         <  1:4310/72‑1 (MQ=255)
                                                                                                                                       tGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAg                                                                                         >  2:4310/1‑72 (MQ=255)
                                                                                                                                       tGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATaaa                                                         <  2:394245/104‑1 (MQ=255)
                                                                                                                                       tGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTATCGCCGa                >  3:476957/1‑145 (MQ=255)
                                                                                                                                       tGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTATCGCCGa                <  2:1226809/145‑1 (MQ=255)
                                                                                                                                          tttttGGTGCCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTATCGCCGATAAGCAg         >  1:1006158/1‑149 (MQ=255)
                                                                                                                                               ggCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTATCGCCGATAAGCAGGCGGa    <  4:1122512/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GCTGATTTTTTCATCGCTCAAGGCGGGCCGTGTAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTATCGCCGATAAGCAGGCGGATA  >  NC_000913/4021862‑4022155

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