Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 457,654 T→A 10.7% Y77N (TAC→AAC)  clpX → ATPase and specificity subunit of ClpX‑ClpP ATP‑dependent serine protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913457,6540TA10.7% 358.8 / 29.3 112Y77N (TAC→AAC) clpXATPase and specificity subunit of ClpX‑ClpP ATP‑dependent serine protease
Reads supporting (aligned to +/- strand):  ref base T (50/50);  new base A (5/7);  total (55/57)
Fisher's exact test for biased strand distribution p-value = 7.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

Only 100 of 117 total aligned reads displayed.

GCTGATTGCCGGTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAAGTAACATTCTGCTGATCGGTCCGACCGGTTCCGGTAAAA  >  NC_000913/457506‑457801
                                                                                                                                                    |                                                                                                                                                   
gCTGATTGCCGGTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGAtt                                                                                                                                                     <  2:731401/149‑1 (MQ=255)
gCTGATTGCCGGTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGAtt                                                                                                                                                     <  2:80420/149‑1 (MQ=255)
  tGATTGCCGGTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTAc                                                                                                                                                   <  3:988078/149‑1 (MQ=255)
  tGATTGCCGGTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTAc                                                                                                                                                   <  2:584945/149‑1 (MQ=255)
    aTTGCCGGTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGt                                                                                                                                                 >  3:678230/1‑149 (MQ=255)
    aTTGCCGGTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGt                                                                                                                                                 <  4:678230/149‑1 (MQ=255)
            tCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGCCGATAACGTAATCGTcc                                                                                                                                         <  1:611431/149‑1 (MQ=255)
            tCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGCCGATAACGTAATCGTcc                                                                                                                                         >  1:396751/1‑149 (MQ=255)
            tCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATc                                                                                                                                             >  3:1383549/1‑145 (MQ=255)
            tCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATc                                                                                                                                             <  4:1383549/145‑1 (MQ=255)
             ccATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGCCGATAACGTAATCGTCCa                                                                                                                                        <  2:396751/149‑1 (MQ=255)
             ccATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGCCGATAACGTAATCGTCCa                                                                                                                                        <  1:1369541/149‑1 (MQ=255)
               aTCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGCCGATAACGTAATCGTCCAg                                                                                                                                       >  4:725982/1‑148 (MQ=255)
               aTCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAgg                                                                                                                                      <  1:2437/149‑1 (MQ=255)
                tCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGa                                                                                                                                     >  2:170069/1‑149 (MQ=255)
                tCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGa                                                                                                                                     >  3:551292/1‑149 (MQ=255)
                       atatCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCg                                                                                                                              <  1:170069/149‑1 (MQ=255)
                               gacgaATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGt                                                                                                                      >  2:708558/1‑149 (MQ=255)
                                   aaTGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTg                                                                                                                  <  4:551292/149‑1 (MQ=255)
                                   aaTGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGCACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTg                                                                                                                  <  3:952969/149‑1 (MQ=255)
                                       tgtTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCgg                                                                                                              >  1:770337/1‑149 (MQ=255)
                                         ttGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGTAACCACCTGGACGATTACGCTATc                                                                                                                                             <  1:1241994/116‑1 (MQ=255)
                                         ttGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTc                                                                                                            >  2:646738/1‑149 (MQ=255)
                                            aTTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGCCGATAACGTAATCGTCCAgg                                                                                                                                      <  2:296777/120‑1 (MQ=255)
                                               tATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGta                                                                                                      <  2:1286482/149‑1 (MQ=255)
                                                 tGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGtata                                                                                                    <  2:770337/149‑1 (MQ=255)
                                                   tAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACa                                                                                                  >  4:919645/1‑149 (MQ=255)
                                                    aaCGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGAtt                                                                                                                                                     >  4:90540/1‑97 (MQ=255)
                                                    aaCGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGAtt                                                                                                                                                     <  3:90540/97‑1 (MQ=255)
                                                    aaCGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATc                                                                                                                                             <  1:243520/105‑1 (MQ=255)
                                                    aaCGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATc                                                                                                                                             >  2:243520/1‑105 (MQ=255)
                                                    aaCGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCa                                                                                                                                        <  2:308875/110‑1 (MQ=255)
                                                    aaCGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCa                                                                                                                                        >  1:308875/1‑110 (MQ=255)
                                                       gACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGCCGATAACGTAATCGTCCAg                                                                                                                                       <  3:873948/108‑1 (MQ=255)
                                                       gACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGCCGATAACGTAATCGTCCAg                                                                                                                                       >  4:873948/1‑108 (MQ=255)
                                                       gACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTa                                                                                                                                               >  4:1009404/1‑100 (MQ=255)
                                                       gACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTa                                                                                                                                               <  3:1009404/100‑1 (MQ=255)
                                                       gACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCg                                                                                                                                            >  4:981344/1‑103 (MQ=255)
                                                       gACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCg                                                                                                                                            <  3:981344/103‑1 (MQ=255)
                                                       gACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCa                                                                                              >  4:601546/1‑149 (MQ=255)
                                                        aCATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTAc                                                                                                                                                   >  4:153042/1‑95 (MQ=255)
                                                        aCATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTAc                                                                                                                                                   <  3:153042/95‑1 (MQ=255)
                                                            catTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCgg                                                                                                                                           <  1:204487/99‑1 (MQ=255)
                                                            catTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACa                                                                                         >  2:1408979/1‑149 (MQ=255)
                                                                 gcgAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGt                                                                                    <  1:708558/149‑1 (MQ=255)
                                                                     agagATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCg                                                                                                               >  1:749548/1‑118 (MQ=255)
                                                                          ttAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAAc                                                                           >  2:39998/1‑149 (MQ=255)
                                                                            aaaGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTAc                                                                                                                                                   <  1:860140/75‑1 (MQ=255)
                                                                            aaaGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTAc                                                                                                                                                   >  2:860140/1‑75 (MQ=255)
                                                                              agaagTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGCCGATAACg                                                                                                                                                  >  2:1278255/1‑74 (MQ=255)
                                                                              agaagTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGCCGATAACg                                                                                                                                                  <  1:1278255/74‑1 (MQ=255)
                                                                              agaagTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCAt                                                                                             >  2:851009/1‑127 (MQ=255)
                                                                              agaagTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCAt                                                                                             <  1:851009/127‑1 (MQ=255)
                                                                              agaagTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCg                                                                       >  3:1334633/1‑149 (MQ=255)
                                                                               gaagTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACa                                                                                                                                  >  4:143581/1‑89 (MQ=255)
                                                                               gaagTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACa                                                                                                                                  <  3:143581/89‑1 (MQ=255)
                                                                               gaagTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGa                                                                      >  2:849102/1‑149 (MQ=255)
                                                                                  gTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATAc                                                                   <  1:475850/149‑1 (MQ=255)
                                                                                   ttGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGc                                                                                                                               >  4:1095065/1‑88 (MQ=255)
                                                                                   ttGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGc                                                                                                                               <  3:1095065/88‑1 (MQ=255)
                                                                                       aCCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCa                                                              >  2:1204365/1‑149 (MQ=255)
                                                                                       aCCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCa                                                              >  3:730296/1‑149 (MQ=255)
                                                                                           cATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATgg                                                          >  4:979754/1‑149 (MQ=255)
                                                                                            aTCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGc                                                         <  1:1204365/149‑1 (MQ=255)
                                                                                                     gCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACaa                                                                                        >  1:905154/1‑109 (MQ=255)
                                                                                                      cAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTgg                                               <  3:919645/149‑1 (MQ=255)
                                                                                                       aGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTggg                                              >  1:986170/1‑149 (MQ=255)
                                                                                                            gcTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCaaaa                                         <  1:39998/149‑1 (MQ=255)
                                                                                                            gcTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCaaaa                                         <  1:646738/149‑1 (MQ=255)
                                                                                                                ccGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGc                                                                                                                <  1:41153/74‑1 (MQ=255)
                                                                                                                ccGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGc                                                                                                                >  2:41153/1‑74 (MQ=255)
                                                                                                                ccGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAAGTaa                                     <  3:601546/149‑1 (MQ=255)
                                                                                                                  gACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCaa                                                             <  3:150598/123‑1 (MQ=255)
                                                                                                                  gACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCaa                                                             >  4:150598/1‑123 (MQ=255)
                                                                                                                  gACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAAGTAACa                                   >  3:27855/1‑149 (MQ=255)
                                                                                                                    cgccgcATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAAGTAACAtt                                 <  4:730296/149‑1 (MQ=255)
                                                                                                                     gccgcATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAAGTAACATTc                                <  2:986170/149‑1 (MQ=255)
                                                                                                                         cATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCa                                                                 >  2:1119488/1‑112 (MQ=255)
                                                                                                                         cATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCa                                                                 <  1:1119488/112‑1 (MQ=255)
                                                                                                                         cATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTgg                                               <  2:681770/130‑1 (MQ=255)
                                                                                                                         cATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTgg                                               >  1:681770/1‑130 (MQ=255)
                                                                                                                          aTGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAAGTAACATTctgctg                           >  2:617232/1‑149 (MQ=255)
                                                                                                                            gAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTc                                                      >  3:523384/1‑120 (MQ=255)
                                                                                                                            gAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTc                                                      <  4:523384/120‑1 (MQ=255)
                                                                                                                             aaaTTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAAGTAACATTCTGCTGATc                        >  4:1342517/1‑149 (MQ=255)
                                                                                                                               aTTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTgg                                               >  3:1311299/1‑124 (MQ=255)
                                                                                                                               aTTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTgg                                               <  4:1311299/124‑1 (MQ=255)
                                                                                                                                ttCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAAGTAACATTCTGCTGATCGGt                     >  4:1134504/1‑149 (MQ=255)
                                                                                                                                     aaCCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACTAGCAATGGCGTCGAGTTGGGCAAAAGTAACATTCTGCTGATCGGTcc                   <  3:80320/146‑1 (MQ=255)
                                                                                                                                     aaCCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAAGTAACATTCTGCTGATCGGTCCGAc                <  3:1134504/149‑1 (MQ=255)
                                                                                                                                           cTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCaaaa                                         >  1:1419063/1‑118 (MQ=255)
                                                                                                                                           cTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCaaaa                                         <  2:1419063/118‑1 (MQ=255)
                                                                                                                                            tGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAAGTAACATTCTGCTGATCGGTCCGACCGGTTcc         >  3:691412/1‑149 (MQ=255)
                                                                                                                                              gACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTgg                                               >  4:308786/1‑109 (MQ=255)
                                                                                                                                              gACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTgg                                               <  3:308786/109‑1 (MQ=255)
                                                                                                                                              gACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAAGTAACATTCTGCTGATCGGTCCGACCGGTTCCgg       >  3:449532/1‑149 (MQ=255)
                                                                                                                                                 gATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAAGTAACATTCTGCTGa                          <  1:1359885/127‑1 (MQ=255)
                                                                                                                                                 gATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAAGTAACATTCTGCTGa                          >  2:1359885/1‑127 (MQ=255)
                                                                                                                                                 gATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCCCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAAGTAACATTCTGCTGATCGGTCCGACCGGTTCCGGTaa    >  1:797337/1‑149 (MQ=255)
                                                                                                                                                   ttACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAAGTAACATTCTGCTGATCGGTCCGACCGGTTCCGGTaaaa  >  1:271147/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCTGATTGCCGGTCCATCCGTGTATATCTGCGACGAATGTGTTGATTTATGTAACGACATCATTCGCGAAGAGATTAAAGAAGTTGCACCGCATCGTGAACGCAGTGCGCTACCGACGCCGCATGAAATTCGCAACCACCTGGACGATTACGTTATCGGCCAGGAACAGGCGAAAAAAGTGCTGGCGGTCGCGGTATACAACCATTACAAACGTCTGCGCAACGGCGATACCAGCAATGGCGTCGAGTTGGGCAAAAGTAACATTCTGCTGATCGGTCCGACCGGTTCCGGTAAAA  >  NC_000913/457506‑457801

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