Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,001,732 Δ1 bp coding (2517/2601 nt) elfC → putative fimbrial usher protein ElfC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,001,7300C.99.6% 1204.9 / ‑2.3 232coding (2515/2601 nt)elfCputative fimbrial usher protein ElfC
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (103/128);  total (104/128)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 240 total aligned reads displayed.

CGGCGTACAGGCTAATAGCGGCATAATTGATGATGATGGCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCCGATCAAATTTGTCATATCAATTACGAGCTTACCGAACAACAAATTAACTCTGCGATTACGCGAATGGATGCCATATGCAGATAATCTTTGGAGA  >  NC_000913/1001630‑1001826
                                                                                                    |                                                                                                
cGGCGTACAGGCTAATAGCGGCATAATTGATGATGATGGCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTc                                                                                                  >  1:1545068/1‑101 (MQ=255)
 ggCGTACAGGCTAATAGCGGCATAATTGATGATGATGGCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTcc                                                                                                 <  2:1544693/101‑1 (MQ=255)
   cGTACAGGCTAATAGCGGCATAATTGATGATGATGGCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTccga                                                                                               <  2:1544731/101‑3 (MQ=255)
    gTACAGGCTAATAGCGGCATAATTGATGATGATGGCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTccgat                                                                                              >  1:1545074/1‑98 (MQ=255)
      aCAGGCTAATAGCGGCATAATTGATGATGATGGCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCa                                                                                           <  2:1544631/101‑1 (MQ=255)
        aGGCTAATAGCGGCATAATTGATGATGATGGCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCaaa                                                                                         <  1:1544317/101‑1 (MQ=255)
         ggCTAATAGCGGCATAATTGATGATGATGGCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAAt                                                                                        <  1:1544475/101‑1 (MQ=255)
           cTAATAGCGGCATAATTGATGATGATGGCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAAttt                                                                                      >  1:1545086/1‑101 (MQ=255)
             aaTAGCGGCATAATTGATGATGATGGCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGt                                                                                    >  1:1545089/1‑101 (MQ=255)
              aTAGCGGCATAATTGATGATGATGGCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTc                                                                                   <  1:1544583/101‑1 (MQ=255)
                 gCGGCATAATTGATGATGATGGCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCata                                                                                <  2:1544766/101‑1 (MQ=255)
                 gCGGCATAATTGATGATGATGGCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCata                                                                                >  1:1545096/1‑101 (MQ=255)
                    gCATAATTGATGATGATGGCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCa                                                                             <  2:1544792/101‑1 (MQ=255)
                     cATAATTGATGATGATGGCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCaa                                                                            <  1:1544505/101‑1 (MQ=255)
                     cATAATTGATGATGATGGCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCaa                                                                            <  1:1544803/101‑1 (MQ=255)
                     cATAATTGATGATGATGGCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCaa                                                                            <  1:1544660/101‑1 (MQ=255)
                     cATAATTGATGATGATGGCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCaa                                                                            <  1:1544864/101‑1 (MQ=255)
                       tAATTGATGATGATGGCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAAtt                                                                          <  1:1544729/101‑1 (MQ=255)
                       tAATTGATGATGATGGCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAAtt                                                                          >  1:1545103/1‑101 (MQ=255)
                        aaTTGATGATGATGGCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTa                                                                         <  1:1544671/101‑1 (MQ=255)
                           tgatgatgaTGGCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACgg                                                                      <  1:1544679/101‑2 (MQ=255)
                           tgatgatgaTGGCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACgg                                                                      <  1:1544591/101‑2 (MQ=255)
                           tgatgatgaTGGCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACgg                                                                      >  1:1545108/1‑100 (MQ=255)
                            gatgatgatGGCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACggg                                                                     <  1:1544888/101‑3 (MQ=255)
                              tgatgatgGCTCGCTCTATATGGCGGGGTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCt                                                                   >  1:1545111/1‑101 (MQ=255)
                                atgatgGCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCata                                                                                >  1:1545113/1‑86 (MQ=255)
                                atgatgGCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTa                                                                 <  2:1544622/101‑1 (MQ=255)
                                atgatgGCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTa                                                                 <  2:1544711/101‑1 (MQ=255)
                                   atgGCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCg                                                              >  1:1545116/1‑101 (MQ=255)
                                     ggCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGaa                                                            <  2:1544611/101‑1 (MQ=255)
                                     ggCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGaa                                                            <  1:1544718/101‑1 (MQ=255)
                                      gctcgctcTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGaac                                                           <  2:1544813/101‑1 (MQ=255)
                                       ctcgctcTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGaaca                                                          <  2:1544851/101‑1 (MQ=255)
                                       ctcgctcTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGaaca                                                          <  2:1544722/101‑1 (MQ=255)
                                        tcgctcTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGaacaa                                                         >  1:1545119/1‑101 (MQ=255)
                                          gctcTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAacaaca                                                       >  1:1545122/1‑101 (MQ=255)
                                            tctATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACaaa                                                     <  1:1544840/101‑1 (MQ=255)
                                            tctATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACaaa                                                     >  1:1545126/1‑101 (MQ=255)
                                             ctATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAAt                                                    >  1:1545128/1‑101 (MQ=255)
                                              tataTGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAAtt                                                   >  1:1545130/1‑101 (MQ=255)
                                               atatGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTa                                                  >  2:1545133/1‑101 (MQ=255)
                                                tatGGCGGTTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTaa                                                 <  2:1544678/101‑1 (MQ=255)
                                                 atGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAAcaacaa                                                      >  1:1545137/1‑95 (MQ=255)
                                                 atGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAAc                                                <  2:1544590/101‑1 (MQ=255)
                                                  tGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAAct                                               <  2:1544741/101‑1 (MQ=255)
                                                  tGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAAct                                               >  2:1545138/1‑101 (MQ=255)
                                                  tGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAAct                                               <  1:1544641/101‑1 (MQ=255)
                                                    gCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAActct                                             <  2:1544703/101‑1 (MQ=255)
                                                     cGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTg                                            >  1:1545140/1‑101 (MQ=255)
                                                     cGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTg                                            <  2:1544574/101‑1 (MQ=255)
                                                     cGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTg                                            <  1:1544634/101‑1 (MQ=255)
                                                     cGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTg                                            <  2:1544981/101‑1 (MQ=255)
                                                       ggTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCg                                          <  2:1544893/101‑1 (MQ=255)
                                                          ttACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGAtt                                       >  1:1545144/1‑101 (MQ=255)
                                                          ttACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGAtt                                       <  1:1544753/101‑1 (MQ=255)
                                                           tACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTa                                      <  1:1544551/101‑1 (MQ=255)
                                                            aCCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTAc                                     <  2:1544858/101‑1 (MQ=255)
                                                             ccGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTAcg                                    >  2:1545147/1‑101 (MQ=255)
                                                             ccGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTAcg                                    >  1:1545146/1‑101 (MQ=255)
                                                              cGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTAcgc                                   >  1:1545149/1‑101 (MQ=255)
                                                               ggCGAAGGGAACAATACGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTAcgcg                                  <  1:1544694/101‑1 (MQ=255)
                                                               ggCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTAcgcg                                  >  1:1545153/1‑101 (MQ=255)
                                                               ggCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTAcgcg                                  <  1:1544757/101‑1 (MQ=255)
                                                                 cGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTACGCGaa                                >  2:1545156/1‑101 (MQ=255)
                                                                 cGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTACGCGaa                                <  1:1544897/101‑1 (MQ=255)
                                                                 cGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTACGCGaa                                <  1:1544736/101‑1 (MQ=255)
                                                                  gAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTACGCGAAt                               >  1:1545159/1‑101 (MQ=255)
                                                                  gAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTACGCGAAt                               <  1:1544603/101‑1 (MQ=255)
                                                                  gAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTACGCGAAt                               >  2:1545158/1‑101 (MQ=255)
                                                                  gAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTACGCGAAt                               <  2:1544758/101‑1 (MQ=255)
                                                                   aaGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTACGCGAATg                              <  1:1544824/101‑1 (MQ=255)
                                                                   aaGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTACGCGAATg                              >  2:1545163/1‑101 (MQ=255)
                                                                      ggAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTACGCGAATGGAt                           <  1:1544568/101‑1 (MQ=255)
                                                                       gAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTACGCGAATGGATg                          >  2:1545167/1‑101 (MQ=255)
                                                                       gAACAATAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTACGCGAATGGATg                          >  1:1545166/1‑101 (MQ=255)
                                                                           aaTAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTACGCGAATGGATGCCat                      >  2:1545169/1‑101 (MQ=255)
                                                                            aTAAGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTACGCGAATGGATGCCata                     >  1:1545172/1‑101 (MQ=255)
                                                                              aaGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTACGCGAATGGATGCCATATg                   <  2:1544822/101‑1 (MQ=255)
                                                                               aGCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTACGCGAATGGATGCCATATGc                  >  2:1545175/1‑101 (MQ=255)
                                                                                gCGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTACGCGAATGGATGCCATATGCa                 <  2:1544734/101‑1 (MQ=255)
                                                                                 cGTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTACGCGAATGGATGCCATATGCAg                <  2:1544616/101‑1 (MQ=255)
                                                                                  gTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTACGCGAATGGATGc                         <  2:1544799/91‑1 (MQ=255)
                                                                                  gTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTACGCGAATGGATGCCATATGCAGa               >  1:1545177/1‑101 (MQ=255)
                                                                                  gTGCGCTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTACGCGAATGGATGCCATATGCAGa               >  2:1545178/1‑101 (MQ=255)
                                                                                    gcgcTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTACGCGAATGGATGCCATATGCAGATa             <  1:1544812/101‑1 (MQ=255)
                                                                                      gcTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTACGCGAATGGATGCCATATGCAGATAAt           <  2:1544774/101‑1 (MQ=255)
                                                                                      gcTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTACGCGAATGGATGCCATATGCAGATAAt           >  2:1545183/1‑101 (MQ=255)
                                                                                       cTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTACGCGAATGGATGCCATATGCAGATAATc          <  1:1544743/101‑1 (MQ=255)
                                                                                       cTGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTACGCGAATGGATGCCATATGCAGATAATc          >  2:1545185/1‑101 (MQ=255)
                                                                                        tGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTACGCGAATGGATGCCATATGCa                 <  2:1544667/93‑1 (MQ=255)
                                                                                        tGGGGCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTACGCGAATGGATGCCATATGCAGATAATCt         >  1:1545187/1‑101 (MQ=255)
                                                                                         ggggCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTACGCGAATGGATGCCATATGCAGATAATCtt        >  2:1545189/1‑101 (MQ=255)
                                                                                          gggCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTACGCGAATGGATGCCATATGCAGATAATCttt       >  1:1545193/1‑101 (MQ=255)
                                                                                          gggCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTACGCGAATGGATGCCATATGCAGATAATCttt       >  2:1545192/1‑101 (MQ=255)
                                                                                           ggCGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTACGCGAATGGATGCCATATGCAGATAATCTTTg      <  2:1544841/101‑1 (MQ=255)
                                                                                             cGAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTACGCGAATGGATGCCATATGCAGATAATCTTTGga    <  1:1544997/101‑1 (MQ=255)
                                                                                              gAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTACGCGAATGGATGCCATATGCAGATAATCTTTGgag   >  1:1545197/1‑101 (MQ=255)
                                                                                              gAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTACGCGAATGGATGCCATATGCAGATAATCTTTGgag   >  1:1545199/1‑101 (MQ=255)
                                                                                              gAAGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTACGCGAATGGATGCCATATGCAGATAATCTTTGgag   >  2:1545200/1‑101 (MQ=255)
                                                                                               aaGCTCCGATCAAATTTGTCATATCAATTACGGGCTTACCGAACAACAAATTAACTCTGCGATTACGCGAATGGATGCCATATGCAGATAATCTTTGgaga  >  1:1545203/1‑101 (MQ=255)
                                                                                                    |                                                                                                
CGGCGTACAGGCTAATAGCGGCATAATTGATGATGATGGCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGAACAATAAGCGTGCGCTGGGGCGAAGCTCCCGATCAAATTTGTCATATCAATTACGAGCTTACCGAACAACAAATTAACTCTGCGATTACGCGAATGGATGCCATATGCAGATAATCTTTGGAGA  >  NC_000913/1001630‑1001826

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