Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,958,644 C→T 6.7% R34Q (CGG→CAG)  fliG ← flagellar motor switch protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,958,6440CT 6.7% 212.1 / 4.0 74R34Q (CGG→CAG) fliGflagellar motor switch protein
Reads supporting (aligned to +/- strand):  ref base C (35/34);  new base T (3/2);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

TACGCATCAGAACCGACGCCCAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCACGACCTGAGCCACACGGCCCGGGTCTTC  >  NC_002947/4958497‑4958793
                                                                                                                                                   |                                                                                                                                                     
tACGCATCAGAACCGACGCCCAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTACCCAGCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCgc                                                                                                                                                     >  1:16756/1‑150 (MQ=255)
 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                                                                                                                                          >  1:600390/1‑149 (MQ=255)
            ccGACGCCCAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCAGGTGCTGCGGCACATGGGc                                                                                                                                          <  2:600390/149‑1 (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:457690/150‑1 (MQ=255)
                                  cGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCgaga                                                                                                                   <  1:1031210/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                                                                                                     <  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:602970/150‑1 (MQ=255)
                                                                                                    gcgcCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCg                                                 >  2:1347257/1‑150 (MQ=255)
                                                                                                       ccATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTAt                                                          <  1:1440058/138‑1 (MQ=255)
                                                                                                       ccATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTAt                                                          >  2:1440058/1‑138 (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                                                                                          <  1:768096/91‑1 (MQ=255)
                                                                                                                      cccGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGAc                                                                                          >  2:768096/1‑91 (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                              >  1:270424/1‑137 (MQ=255)
                                                                                                                                    ttcGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCa                              <  2:270424/137‑1 (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                                        >  2:538321/1‑117 (MQ=255)
                                                                                                                                              tgtgCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATc                                        <  1:538321/117‑1 (MQ=255)
                                                                                                                                               gtgCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCACGACCTGAGCCACACGGCCCgg        >  1:1280036/1‑148 (MQ=255)
                                                                                                                                               gtgCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCACGACCTGAGCCACACGGCCCgg        <  2:1280036/148‑1 (MQ=255)
                                                                                                                                                tgCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTAt                                                          <  2:477768/97‑1 (MQ=255)
                                                                                                                                                tgCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTAt                                                          >  1:477768/1‑97 (MQ=255)
                                                                                                                                                   cGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCACGACCTGAGCCACACGGCCCGGGTCTTc  <  1:1347257/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TACGCATCAGAACCGACGCCCAGGCTGGTCTGGTCGCCGACGATGTCGACGAACTCGCTCATCACCTGTTCGACCTGGTCACGGTGCACATTGCCCATCTGCGCCATGGCCACACCCACCCGCTGCACTTCCTTCGGGCCCATGTGCCGCAGCACCTGGGCCGCATCGGTCTCGCCCAGCGAGAGCAGGAGGATCGCCGCCTTGTCGACGCGGCTCAGCTTGGCGGTAACGGCTCGGTTATCACTCATCGGCGTTGATCCAGTCTTTCACGACCTGAGCCACACGGCCCGGGTCTTC  >  NC_002947/4958497‑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: