Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,958,647 A→G 6.5% L33P (CTG→CCG)  fliG ← flagellar motor switch protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,958,6470AG 6.5% 218.0 / 2.6 72L33P (CTG→CCG) fliGflagellar motor switch protein
Reads supporting (aligned to +/- strand):  ref base A (32/35);  new base G (3/2);  total (35/37)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCATCAGAACCGACGCCCAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCACGACCTGAGCCACACGGCCCGGGTCTTC  >  NC_002947/4958498‑4958793
                                                                                                                                                     |                                                                                                                                                  
aCGCATCAGAACCGACGCCCAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCgca                                                                                                                                                    >  1:1094420/1‑150 (MQ=255)
  gCATCAGAACCGACGCCCAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCgcagc                                                                                                                                                  >  1:10406/1‑150 (MQ=255)
      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                                                                                                                                               <  1:217428/96‑1 (MQ=255)
                                                           atcaCCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCAcc                                                                                                                                               >  2:217428/1‑96 (MQ=255)
                                                           atcaCCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGAcg                                                                                         >  2:1452082/1‑150 (MQ=255)
                                                            tcaCCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGAcgc                                                                                        >  2:454538/1‑150 (MQ=255)
                                                               ccTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGt                                                                                              >  2:1014329/1‑141 (MQ=255)
                                                               ccTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGt                                                                                              <  1:1014329/141‑1 (MQ=255)
                                                                cTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCt                                                                                    >  1:427859/1‑150 (MQ=255)
                                                                  gTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGa                                                                                           <  1:941992/141‑1 (MQ=255)
                                                                  gTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGa                                                                                           >  2:941992/1‑141 (MQ=255)
                                                                   ttCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAg                                                                                 <  2:416998/150‑1 (MQ=255)
                                                                      gACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCtt                                                                              >  2:1040629/1‑150 (MQ=255)
                                                                       aCCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATcgcc                                                                                                     <  2:1321891/126‑1 (MQ=255)
                                                                       aCCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATcgcc                                                                                                     >  1:1321891/1‑126 (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                                                         >  1:1169414/1‑132 (MQ=255)
                                                                                                             cacaCCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATc                                                         <  2:1169414/132‑1 (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                                                >  2:35045/1‑120 (MQ=255)
                                                                                                                                  cttcGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCgg                                                <  1:35045/120‑1 (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                                                          <  2:477768/97‑1 (MQ=255)
                                                                                                                                               tgCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTAt                                                          >  1:477768/1‑97 (MQ=255)
                                                                                                                                                  cGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCACGACCTGAGCCACACGGCCCGGGTCTTc  <  1:1347257/150‑1 (MQ=255)
                                                                                                                                                   gcagcaCCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATCCa                                      >  1:1211028/1‑113 (MQ=255)
                                                                                                                                                   gcagcaCCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATCCa                                      <  2:1211028/113‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
ACGCATCAGAACCGACGCCCAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCACGACCTGAGCCACACGGCCCGGGTCTTC  >  NC_002947/4958498‑4958793

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: