Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,958,652 C→A 6.8% 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,6520CA 6.8% 230.2 / 5.4 73Q31H (CAG→CATfliGflagellar motor switch protein
Reads supporting (aligned to +/- strand):  ref base C (33/35);  new base A (3/2);  total (36/37)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.36e-01

CAGAACCGACGCCCAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCACGACCTGAGCCACACGGCCCGGGTCTTCGGCCACC  >  NC_002947/4958504‑4958800
                                                                                                                                                    |                                                                                                                                                    
cAGAACCGACGCCCAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCt                                                                                                                                                     <  2:10406/150‑1 (MQ=255)
 aGAACCGACGCCCAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTg                                                                                                                                                    >  2:457690/1‑150 (MQ=255)
  gAACCGACGCCCAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTgg                                                                                                                                                   >  2:1239974/1‑150 (MQ=255)
     ccGACGCCCAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCAGGTGCTGCGGCACATGGGc                                                                                                                                                 <  2:600390/149‑1 (MQ=255)
     ccGACGCCCAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCAGGTGCTGCGGCACATGGGc                                                                                                                                                 >  1:600390/1‑149 (MQ=255)
      cGACGCCCAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCg                                                                                                                                               <  1:1239974/150‑1 (MQ=255)
             cAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCCCCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGt                                                                                                                                        <  1:177726/150‑1 (MQ=255)
                     tctggtcGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTTGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCa                                                                                                                                <  1:1048287/150‑1 (MQ=255)
                       tggtcGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGc                                                                                                                              <  2:1094420/150‑1 (MQ=255)
                           cGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCgaga                                                                                                                          <  1:1031210/150‑1 (MQ=255)
                           cGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCgaga                                                                                                                          <  1:457690/150‑1 (MQ=255)
                              cgacgaTGTCGACGAACTCGCTCATCATCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCa                                                                                                                       >  2:114839/1‑150 (MQ=255)
                                      tCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATc                                                                                                               >  2:1191793/1‑150 (MQ=255)
                                              cTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCtt                                                                                                       <  1:567634/150‑1 (MQ=255)
                                                     atcaCCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCAcc                                                                                                                                                      >  2:217428/1‑96 (MQ=255)
                                                     atcaCCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCAcc                                                                                                                                                      <  1:217428/96‑1 (MQ=255)
                                                     atcaCCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGAcg                                                                                                >  2:1452082/1‑150 (MQ=255)
                                                      tcaCCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGAcgc                                                                                               >  2:454538/1‑150 (MQ=255)
                                                         ccTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGt                                                                                                     <  1:1014329/141‑1 (MQ=255)
                                                         ccTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGt                                                                                                     >  2:1014329/1‑141 (MQ=255)
                                                          cTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCt                                                                                           >  1:427859/1‑150 (MQ=255)
                                                            gTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGa                                                                                                  >  2:941992/1‑141 (MQ=255)
                                                            gTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGa                                                                                                  <  1:941992/141‑1 (MQ=255)
                                                             ttCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAg                                                                                        <  2:416998/150‑1 (MQ=255)
                                                                gACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCtt                                                                                     >  2:1040629/1‑150 (MQ=255)
                                                                 aCCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATcgcc                                                                                                            >  1:1321891/1‑126 (MQ=255)
                                                                 aCCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATcgcc                                                                                                            <  2:1321891/126‑1 (MQ=255)
                                                                 aCCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTg                                                                                    <  1:454538/150‑1 (MQ=255)
                                                                  ccTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTgg                                                                                   <  1:1234891/150‑1 (MQ=255)
                                                                  ccTGGTCACGCTGCACATTGCCCATCTGCCCCATGGCCAAACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTgg                                                                                   <  1:554665/150‑1 (MQ=255)
                                                                      gTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGt                                                                               >  1:1463841/1‑150 (MQ=255)
                                                                       tCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCgaga                                                                                                                          >  1:946980/1‑106 (MQ=255)
                                                                       tCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTCCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCgaga                                                                                                                          <  2:946980/106‑1 (MQ=255)
                                                                                  aTTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCGTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGc                                                                                                                              <  1:623296/91‑1 (MQ=255)
                                                                                  aTTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGc                                                                                                                              >  2:623296/1‑91 (MQ=255)
                                                                                  aTTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGtt                                                                   <  2:427859/150‑1 (MQ=255)
                                                                                    tGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTAt                                                                 <  1:1452082/150‑1 (MQ=255)
                                                                                            tGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATTTGCCGCAGCACCTGGGCCGCATCGGTCTCGACCAACGAGAGCAGGAGGAGCGCCGCCTGGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATc                                                         <  2:1548913/150‑1 (MQ=255)
                                                                                             gcgcCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCg                                                        >  2:1347257/1‑150 (MQ=255)
                                                                                             gcgcCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCg                                                        <  2:602970/150‑1 (MQ=255)
                                                                                                ccATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTAt                                                                 >  2:1440058/1‑138 (MQ=255)
                                                                                                ccATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTAt                                                                 <  1:1440058/138‑1 (MQ=255)
                                                                                                ccATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCg                                                     <  2:1463841/150‑1 (MQ=255)
                                                                                                 cATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGCGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGt                                                    <  1:114839/150‑1 (MQ=255)
                                                                                                       cacaCCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATc                                                                <  2:1169414/132‑1 (MQ=255)
                                                                                                       cacaCCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATc                                                                >  1:1169414/1‑132 (MQ=255)
                                                                                                        acacCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCACGGAACTCCCGGCATTGTCGACGCGGTTCACCTTCGCAGCCAGGGCTAGGTTGCCACCCATCGGCGTTGATCCa                                             >  1:28719/1‑150 (MQ=255)
                                                                                                          acccacccGCTGCACTTCCTTCGGGCCTTTGAGGCCTAGCAGCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTCACGGCTCGGTTATCACTCATCGGCGTTGATCCAGt                                           <  2:130622/150‑1 (MQ=255)
                                                                                                            ccacccGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCCCTCATCGGCGTTGCTCCCGTCt                                         >  2:424688/1‑150 (MQ=255)
                                                                                                              acccGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATCCAGTCttt                                       >  1:898963/1‑150 (MQ=255)
                                                                                                               cccGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGAc                                                                                                 >  2:768096/1‑91 (MQ=255)
                                                                                                               cccGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGAc                                                                                                 <  1:768096/91‑1 (MQ=255)
                                                                                                               cccGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTc                                      >  2:1644694/1‑150 (MQ=255)
                                                                                                                       ccttccttcGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGAGGTCACGGATCGGCTATCACTGGTCGGCGGTGAGTCAGTCTTGCACGACCTg                              <  2:16756/149‑1 (MQ=255)
                                                                                                                         ttccttcGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCACGACCTGAg                            >  2:410018/1‑150 (MQ=255)
                                                                                                                            cttcGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCgg                                                       <  1:35045/120‑1 (MQ=255)
                                                                                                                            cttcGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCgg                                                       >  2:35045/1‑120 (MQ=255)
                                                                                                                            cttcGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCAACGCCGGCCAGCTTGGCGGGAACGGCTCGGGTATCACTTATCGCCGGTGATCCAGTCTTTCACGACCTGTGCca                         >  1:1263983/1‑150 (MQ=255)
                                                                                                                             ttcGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCa                                     <  2:270424/137‑1 (MQ=255)
                                                                                                                             ttcGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCa                                     >  1:270424/1‑137 (MQ=255)
                                                                                                                             ttcGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCACGACCTGAGCcac                        >  1:803435/1‑150 (MQ=255)
                                                                                                                                  gCCCAGGTGCTGCGGCACATGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCg                                                                                 <  1:673968/88‑1 (MQ=38)
                                                                                                                                  gCCCAGGTGCTGCGGCACATGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCg                                                                                 >  2:673968/1‑88 (MQ=38)
                                                                                                                                    ccAGGTGCTGCGGCACATGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCACGACCTGAGCCACACGGccc                 >  2:241781/1‑150 (MQ=255)
                                                                                                                                     cATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCACGACCTGAGCCACACGGCCCg                >  2:75047/1‑150 (MQ=255)
                                                                                                                                       tgtgCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATc                                               <  1:538321/117‑1 (MQ=255)
                                                                                                                                       tgtgCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATc                                               >  2:538321/1‑117 (MQ=255)
                                                                                                                                        gtgCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCACGACCTGAGCCACACGGCCCgg               <  2:1280036/148‑1 (MQ=255)
                                                                                                                                        gtgCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCACGACCTGAGCCACACGGCCCgg               >  1:1280036/1‑148 (MQ=255)
                                                                                                                                         tgCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTAt                                                                 >  1:477768/1‑97 (MQ=255)
                                                                                                                                         tgCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTAt                                                                 <  2:477768/97‑1 (MQ=255)
                                                                                                                                            cGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCACGACCTGAGCCACACGGCCCGGGTCTTc         <  1:1347257/150‑1 (MQ=255)
                                                                                                                                             gcagcaCCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATCCa                                             >  1:1211028/1‑113 (MQ=255)
                                                                                                                                             gcagcaCCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATCCa                                             <  2:1211028/113‑1 (MQ=255)
                                                                                                                                                  aCCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCACGACCTGAGCCACACGGCCCGGGTCTTCGGccac   >  2:131565/1‑150 (MQ=255)
                                                                                                                                                   ccTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGGCTTGTCGACTCGACTGTGCTTAGCGGTACCGGCTAGGTTATCACCCTTCGCAGTTGCTCCATTCTTCTCCGTCCCGAGCCACACGTGCGGGGTCTTCGTccacc  >  2:283872/1‑150 (MQ=255)
                                                                                                                                                    cTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCgg                                                                                >  2:535506/1‑71 (MQ=255)
                                                                                                                                                    cTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCgg                                                                                <  1:535506/71‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAGAACCGACGCCCAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCACGACCTGAGCCACACGGCCCGGGTCTTCGGCCACC  >  NC_002947/4958504‑4958800

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