Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,022,008 A→G 6.4% T22A (ACC→GCC)  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,0080AG 6.4% 529.4 / 17.9 157T22A (ACC→GCC) tatATatABCE protein translocation system subunit
Reads supporting (aligned to +/- strand):  ref base A (75/72);  new base G (5/5);  total (80/77)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

Only 100 of 159 total aligned reads displayed.

GCTGATTTTTTCATCGCTCAAGGCGGGCCGTGTAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTATCGCCGATAAGCAGGCGGATA  >  NC_000913/4021862‑4022155
                                                                                                                                                  |                                                                                                                                                   
     ttttttCATCGCTCAAGGCGGGCCGTGTAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCaaaaa                                                                                                                                              >  2:192082/1‑149 (MQ=255)
      tttttCATCGCTCAAGGCGGGCCGTGTAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAg                                                                                                                                             <  1:192082/149‑1 (MQ=255)
       ttttCATCGCTCAAGGCGGGCCGTGTAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGc                                                                                                                                            <  1:551809/149‑1 (MQ=255)
            aTCGCTCAAGGCGGGCCGTGTAACGTATAATGCGGCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGc                                                                                                                                            >  2:1081652/1‑144 (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)
                                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)
                                             ggCTTTGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCaaa                                                                                                      <  3:1411291/149‑1 (MQ=255)
                                               cTTTGTTTAATCATCATCTACCACAGATGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAgg                                                                                                    >  3:975922/1‑149 (MQ=255)
                                                tttGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGc                                                                                                   <  3:717976/149‑1 (MQ=255)
                                                  tGTTTAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCtt                                                                                                 <  1:976188/149‑1 (MQ=255)
                                                    tttAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATc                                                                                                                       <  2:1062637/125‑1 (MQ=255)
                                                    tttAATCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATc                                                                                                                       >  1:1062637/1‑125 (MQ=255)
                                                       aaTCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTcc                                                                                                                                    <  1:144938/109‑1 (MQ=255)
                                                       aaTCATCATCTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTcc                                                                                                                                    >  2:144938/1‑109 (MQ=255)
                                                         tcatcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTGTTGCCGTCATCGTTGTACTGCTTTTTGGTGCCAAAAAGca                                                                                                                                           <  4:850785/100‑2 (MQ=255)
                                                          catcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGc                                                                                                   >  3:36023/1‑139 (MQ=255)
                                                           atcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGCTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAATGc                                                                                        <  4:977707/149‑1 (MQ=255)
                                                            tcatcTACCACAGAGGAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCa                                                                                       >  2:946838/1‑149 (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)
                                                                           gAACATGTATGGGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTg                                                                                                                    >  2:195654/1‑105 (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)
                                                                               aTGTATGCGTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTaaa                                                                                             <  4:879885/124‑1 (MQ=255)
                                                                                  tATGGGTGGTAGCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGc                                                                 >  1:766251/1‑149 (MQ=255)
                                                                                     gggTGGTATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGt                                                                                                                  <  2:549077/97‑1 (MQ=255)
                                                                                       gtggtATCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGAt                                                            >  3:38062/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)
                                                                                             tCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGc                                                                                                   <  2:1457875/104‑1 (MQ=255)
                                                                                             tCAGTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGc                                                                                                   >  1:1457875/1‑104 (MQ=255)
                                                                                                gTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGa                                                                                                           >  2:1235227/1‑93 (MQ=255)
                                                                                                gTATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGa                                                                                                           <  1:1235227/93‑1 (MQ=255)
                                                                                                 tATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCg                                                                                                                                         <  4:300057/62‑1 (MQ=255)
                                                                                                 tATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCg                                                                                                                                         >  3:300057/1‑62 (MQ=255)
                                                                                                 tATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATaaaa                                                        >  1:898683/1‑143 (MQ=255)
                                                                                                 tATTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATaaaa                                                        <  2:898683/143‑1 (MQ=255)
                                                                                                  aTTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGt                                                                                                                  >  3:1217785/1‑84 (MQ=255)
                                                                                                  aTTTGGCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGt                                                                                                                  <  4:1217785/84‑1 (MQ=255)
                                                                                                      ggCAGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTc                                                  >  4:707479/1‑144 (MQ=255)
                                                                                                         aGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCt                                          >  4:833164/1‑149 (MQ=255)
                                                                                                         aGTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCt                                          >  3:195480/1‑149 (MQ=255)
                                                                                                          gTTATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCga                                                                               <  1:875277/111‑1 (MQ=255)
                                                                                                           ttATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCgat                                                                              <  2:1433177/111‑1 (MQ=255)
                                                                                                            tATTGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGAt                                       <  1:1090565/149‑1 (MQ=255)
                                                                                                              ttGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGAttt                                     >  1:1357592/1‑149 (MQ=255)
                                                                                                              ttGATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGAttt                                     >  3:1367417/1‑149 (MQ=255)
                                                                                                                gATTATTGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTa                                   <  4:1367417/149‑1 (MQ=255)
                                                                                                                 attattGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTg                                         >  2:1180038/1‑142 (MQ=255)
                                                                                                                 attattGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTAc                                  <  2:766251/149‑1 (MQ=255)
                                                                                                                  ttattGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACt                                 <  1:1177575/149‑1 (MQ=255)
                                                                                                                  ttattGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACt                                 <  3:935419/149‑1 (MQ=255)
                                                                                                                     ttGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGa                                                                                   <  4:1102026/96‑1 (MQ=255)
                                                                                                                     ttGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGa                                                                                   >  3:1102026/1‑96 (MQ=255)
                                                                                                                     ttGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGa                                                                                   <  4:453766/96‑1 (MQ=255)
                                                                                                                     ttGCCGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGAt                                       >  2:1248323/1‑140 (MQ=255)
                                                                                                                     ttGCCGTCATCGTTGTACTCCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGt                                                                                                                  <  3:841953/65‑1 (MQ=255)
                                                                                                                        ccGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGaaa                           >  2:1156125/1‑149 (MQ=255)
                                                                                                                         cGTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTcc                                                                                                                           <  2:1095343/52‑1 (MQ=255)
                                                                                                                          gTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGc                                           <  4:759560/131‑1 (MQ=255)
                                                                                                                          gTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGc                                           >  3:759560/1‑131 (MQ=255)
                                                                                                                          gTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAAc                         >  1:1088018/1‑149 (MQ=255)
                                                                                                                          gTCATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAAc                         <  3:833164/149‑1 (MQ=255)
                                                                                                                            cATCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTa                       >  1:94053/1‑149 (MQ=255)
                                                                                                                             aTCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTaaaaaa                                                                                          >  4:354365/1‑81 (MQ=255)
                                                                                                                              tCGTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTATc                     >  3:649508/1‑149 (MQ=255)
                                                                                                                                gTTGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGAtt                                      <  1:1074201/130‑1 (MQ=255)
                                                                                                                                 ttGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCgat                                                                              <  4:369299/89‑1 (MQ=255)
                                                                                                                                 ttGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGAt                                       >  1:870948/1‑128 (MQ=255)
                                                                                                                                 ttGTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGAt                                       <  2:870948/128‑1 (MQ=255)
                                                                                                                                   gTACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTATCGCCGa                <  2:1088018/149‑1 (MQ=255)
                                                                                                                                    tACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTAc                                  <  1:938869/130‑1 (MQ=255)
                                                                                                                                    tACTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTATCGCCGAt               >  3:173547/1‑149 (MQ=255)
                                                                                                                                    tACTGCTGTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTAc                                  >  2:938869/1‑130 (MQ=255)
                                                                                                                                     aCTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTATCGCCGATa              >  4:937576/1‑149 (MQ=255)
                                                                                                                                     aCTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTATCGCCGATa              <  4:1405590/149‑1 (MQ=255)
                                                                                                                                      cTGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTATCGCCGATaa             <  4:173547/149‑1 (MQ=255)
                                                                                                                                       tGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAg                                                                                         >  2:4310/1‑72 (MQ=255)
                                                                                                                                       tGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAg                                                                                         <  1:4310/72‑1 (MQ=255)
                                                                                                                                       tGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGaa                                                                        >  2:693881/1‑89 (MQ=255)
                                                                                                                                       tGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGaa                                                                        <  1:693881/89‑1 (MQ=255)
                                                                                                                                       tGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTATCGCCGa                >  3:476957/1‑145 (MQ=255)
                                                                                                                                       tGCTTTTTGGCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTATCGCCGa                >  1:1226809/1‑145 (MQ=255)
                                                                                                                                          tttttGGTGCCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTATCGCCGATAAGCAg         >  1:1006158/1‑149 (MQ=255)
                                                                                                                                               ggCACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTATCGCCGATAAGCAGGCGGa    <  4:1122512/149‑1 (MQ=255)
                                                                                                                                                 cACCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAGCAGGATAAAACCAGTCAGGATGCTGATTTTACTGCGAAAACTATCGCCGATAAGCAGGCGGATa  <  4:38062/149‑1 (MQ=255)
                                                                                                                                                  aCCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAg                                                                  <  3:1192735/84‑1 (MQ=255)
                                                                                                                                                  aCCAAAAAGCTCGGCTCCATCGGTTCCGATCTTGGTGCGTCGATCAAAGGCTTTAAAAAAGCAATGAGCGATGATGAACCAAAg                                                                  >  4:1192735/1‑84 (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: