Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,958,607 C→A 7.0% Q31H (CAG→CAT fliG ← flagellar motor switch protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,958,6070CA 7.0% 204.6 / 4.1 65Q31H (CAG→CATfliGflagellar motor switch protein
Reads supporting (aligned to +/- strand):  ref base C (32/28);  new base A (2/3);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

ACCGACGCCCAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAAC‑‑‑‑GGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCACGACCTGAGCCACACGGCCCGGGTCTTCGGCCAC  >  NC_002947/4958463‑4958754
                                                                                                                                                |                                                                                                                                                       
aCCGACGCCCAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTggg                                                                                                                                                     <  7:231622/149‑1 (MQ=255)
aCCGACGCCCAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCAGGTGCTGCGGCACATggg                                                                                                                                                     <  6:134318/149‑1 (MQ=255)
   gACGCCCAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCAGGTGCTGCGGCACATGGGccc                                                                                                                                                  <  8:276440/149‑2 (MQ=255)
    aCGCCCAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGc                                                                                                                                                 <  5:242107/149‑1 (MQ=255)
      gCCCAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCAGGTGCTGCGGCACATGGGcccgaa                                                                                                                                               <  5:109328/149‑5 (MQ=255)
                 tctggtcGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGACCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGcc                                                                                                                                     >  3:156336/1‑148 (MQ=255)
                     gtcGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGc                                                                                                                                 >  3:1324/1‑148 (MQ=255)
                     gtcGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGc                                                                                                                                 <  4:1324/148‑1 (MQ=255)
                        gCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCga                                                                                                                               >  2:254590/1‑147 (MQ=255)
                          cgacgaTGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGagag                                                                                                                            >  3:242043/1‑148 (MQ=255)
                          cgacgaTGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGagag                                                                                                                            >  1:245875/1‑148 (MQ=255)
                             cgaTGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCagg                                                                                                                        <  1:320170/149‑1 (MQ=255)
                             cgaTGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCagg                                                                                                                        <  7:145306/149‑1 (MQ=255)
                             cgaTGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCagg                                                                                                                        <  2:245875/149‑1 (MQ=255)
                                 gTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAggagga                                                                                                                    >  8:136495/1‑149 (MQ=255)
                                    gacgaACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATcg                                                                                                                 >  2:282861/1‑149 (MQ=255)
                                    gacgaACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATcg                                                                                                                 >  3:108727/1‑149 (MQ=255)
                                        aaCTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATcgccgc                                                                                                             <  2:230432/149‑1 (MQ=255)
                                         aCTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCgccgcc                                                                                                            <  3:97806/149‑1 (MQ=255)
                                           tCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCtt                                                                                                          >  6:191679/1‑149 (MQ=255)
                                             gCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGt                                                                                                        <  7:82562/149‑1 (MQ=255)
                                               tcatcaCCTGTTCGACCTGGTCACGGAGCACAATGCCCATCTCCGCCATGGCCACACCCACCAGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACATGGCGCGCAACCGTCTCGCCCAGCGAGAGAAGGAGGATCGCCGCCATGTc                                                                                                       >  8:205012/1‑148 (MQ=255)
                                                 atcaCCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCtt                                                                                                          <  6:242212/143‑1 (MQ=255)
                                                 atcaCCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACACACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCtt                                                                                                          >  5:242212/1‑143 (MQ=255)
                                                  tcaCCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGAcg                                                                                                   >  2:101401/1‑149 (MQ=255)
                                                   caCCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGAccc                                                                                                  >  2:269313/1‑147 (MQ=255)
                                                     ccTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCgg                                                                                                <  3:109437/149‑1 (MQ=255)
                                                     ccTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCgg                                                                                                >  4:225428/1‑149 (MQ=255)
                                                              ccTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTg                                                                                       <  7:65308/149‑1 (MQ=255)
                                                               cTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTgg                                                                                      >  2:191746/1‑149 (MQ=255)
                                                                  gTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGc                                                                                     <  4:144814/147‑1 (MQ=255)
                                                                         tGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAAC‑‑‑‑GGc                                                                        <  4:108727/149‑1 (MQ=255)
                                                                           cacaTTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAAC‑‑‑‑GGCt                                                                       <  5:146756/148‑1 (MQ=255)
                                                                           cacaTTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAAC‑‑‑‑GGCTc                                                                      >  5:256549/1‑149 (MQ=255)
                                                                            acaTTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAAC‑‑‑‑GGCTCg                                                                     >  4:241843/1‑149 (MQ=255)
                                                                             caTTGCCCATCTGCGCCATGGCGACAGCCACGCGCTGCACTTCCTTCGGGCCGATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAAC‑‑‑‑GGCTCgg                                                                    <  1:123975/149‑1 (MQ=255)
                                                                              aTTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAAC‑‑‑‑GGCTCGGt                                                                   >  1:48423/1‑149 (MQ=255)
                                                                                   ccATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAAC‑‑‑‑GGCTCGGTTAt                                                                >  8:204486/1‑147 (MQ=255)
                                                                                   ccATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAAC‑‑‑‑GGCTCGGTTAt                                                                <  7:204486/147‑1 (MQ=255)
                                                                                   ccATCTGCGCCATGGCCACACCCACCCGCTCCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCATGTCGACGCGGCACAGCTTGGCGCTAAC‑‑‑‑GGCTCGGTTATCa                                                              >  2:61684/1‑149 (MQ=255)
                                                                                             cATGGCCACACCCACACGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGAGGTAAC‑‑‑‑GGCTCGGTTATCACTCATCGGCg                                                    >  6:319564/1‑149 (MQ=255)
                                                                                               tGGCCACACCCACCCGCTGCACTTCCTTCGGTCCCATGTGCCGCAGCACCTTGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAAC‑‑‑‑GGCTCGGTTATCACTCATCGGCGtt                                                  <  6:43045/149‑1 (MQ=255)
                                                                                                 gCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAAC‑‑‑‑GGCTCGGTTATCACTCATCGGCGTTGa                                                <  6:256549/149‑1 (MQ=255)
                                                                                                  ccACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAAC‑‑‑‑GGCTCGGTTATCACTCATCGGCGTTGAt                                               <  1:254590/149‑1 (MQ=255)
                                                                                                   cacaCCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAAC‑‑‑‑GGCTCGGTTATCACTCATCGGCGTTGATc                                              <  4:242043/149‑1 (MQ=255)
                                                                                                      acccacccGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAAC‑‑‑‑GGCTCGGTTATCACTCATCGGCGTTGATCCAg                                           >  6:14208/1‑149 (MQ=255)
                                                                                                      acccacccGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAAC‑‑‑‑GGCTCGGTTATCACTCATCGGCCTTGATCCAg                                           >  8:34463/1‑149 (MQ=255)
                                                                                                        ccacccGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAAC‑‑‑‑GGCTCGGt                                                                   >  4:192018/1‑123 (MQ=255)
                                                                                                        ccacccGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAAC‑‑‑‑GGCTCGGt                                                                   <  3:192018/123‑1 (MQ=255)
                                                                                                         cacccGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAAC‑‑‑‑GGCTCGGTTATCACTCCTCGGCGTTGATCCAGTCt                                        >  7:116472/1‑149 (MQ=255)
                                                                                                         cacccGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAAC‑‑‑‑GGCTCGGTTATCACTCATCGGCGTTGATCCAGTCt                                        <  2:48423/149‑1 (MQ=255)
                                                                                                                       ccttcGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAAC‑‑‑‑GGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCACGACCTGa                            <  3:241843/147‑1 (MQ=255)
                                                                                                                        cttcGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAAC‑‑‑‑GGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCACGACCTGAGc                          <  7:136495/148‑1 (MQ=255)
                                                                                                                         ttcGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAAC‑‑‑‑GGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCACGACCTGAGCca                        >  1:63363/1‑149 (MQ=255)
                                                                                                                              gCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAAC‑‑‑‑GGCTCg                                                                     >  1:10314/1‑99 (MQ=255)
                                                                                                                              gCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAAC‑‑‑‑GGCTCg                                                                     <  2:10314/99‑1 (MQ=255)
                                                                                                                               cccAGGTGCTGCGGCACATGGGCCGCATAGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAAC‑‑‑‑GGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCACGACCTGAGCCACACGGc                  >  1:185309/1‑149 (MQ=255)
                                                                                                                                ccATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAAC‑‑‑‑GGCTCGGTTATAACTCATCGGCGTTGATCCAGTCTTTCACGACCTGAGCCACACGGcc                 >  8:248974/1‑149 (MQ=255)
                                                                                                                                     tgCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAAC‑‑‑‑GGCTCGGGTATCACTCTTCGGCGTTGATCCAGTCTTTCACGACCTGAGCCACACGGCCCggg             <  4:156336/148‑1 (MQ=255)
                                                                                                                                      gCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAAC‑‑‑‑GGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCACGACCTGAGCCACACGGCCCGGGTc           >  7:91940/1‑149 (MQ=255)
                                                                                                                                       ccGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAAC‑‑‑‑GGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCACGACCTGAGCCACACGGCCCGGGTCt          >  6:260192/1‑149 (MQ=255)
                                                                                                                                        cGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAAC‑‑‑‑GGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCACGACCTGAGCCACACGGCCCGGGTCtt         >  2:44403/1‑149 (MQ=255)
                                                                                                                                         gcagcaCCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAAC‑‑‑‑GGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCACGACCTGAGCCACACGGCCCGGGTCTTc        >  1:228054/1‑149 (MQ=255)
                                                                                                                                             caCCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAAC‑‑‑‑GGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCACGACCTGAGCCACACGGCCCGGGTCTTCGGcc    >  1:79321/1‑149 (MQ=255)
                                                                                                                                              aCCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTACCTCCTGGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCACGACCTGAGCCACACGGCCCGGGTCTTCg       <  1:101401/149‑1 (MQ=255)
                                                                                                                                               ccTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAAC‑‑‑‑GGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCACGACCTGAGCCACACGGCCCGGGTCTTCGGccac  <  7:34463/149‑1 (MQ=255)
                                                                                                                                                cTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAAC‑‑‑‑GGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCACGACCTGAGCCACACGGCCCGGGTCTTCGGccac  >  5:205377/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                       
ACCGACGCCCAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAAC‑‑‑‑GGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCACGACCTGAGCCACACGGCCCGGGTCTTCGGCCAC  >  NC_002947/4958463‑4958754

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