Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,130,388 G→A A159V (GCA→GTA)  flgA ← assembly protein for flagellar basal‑body periplasmic P ring

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,130,3880GA98.1% 702.6 / ‑1.4 211A159V (GCA→GTA) flgAassembly protein for flagellar basal‑body periplasmic P ring
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base A (86/121);  total (88/123)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 211 total aligned reads displayed.

GACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAGCACCGTACGC  >  NC_000913/1130291‑1130483
                                                                                                 |                                                                                               
  cGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTg                                                                                             >  1:5495294/1‑100 (MQ=255)
   ggCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTgg                                                                                            <  1:1835611/100‑1 (MQ=255)
    gCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGc                                                                                           <  2:2380324/100‑1 (MQ=255)
    gCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGc                                                                                           <  1:1926333/100‑1 (MQ=255)
    gCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGc                                                                                           >  1:1689936/1‑100 (MQ=255)
     cTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCg                                                                                          <  1:1012121/100‑1 (MQ=255)
       gCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGaa                                                                                        >  1:4933978/1‑100 (MQ=255)
        cATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGaaa                                                                                       >  1:274617/1‑100 (MQ=255)
         aTTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAAc                                                                                      <  1:3173959/100‑1 (MQ=255)
            gTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTgg                                                                                   >  2:936581/1‑100 (MQ=255)
            gTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTgg                                                                                   <  2:1501714/100‑1 (MQ=255)
            gTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTgg                                                                                   <  1:241387/100‑1 (MQ=255)
             ttCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTggg                                                                                  >  1:664785/1‑100 (MQ=255)
             ttCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTggg                                                                                  <  2:1811050/100‑1 (MQ=255)
             ttCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTggg                                                                                  <  2:4244688/100‑1 (MQ=255)
             ttCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTggg                                                                                  <  2:391847/100‑1 (MQ=255)
                  cgcCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAAc                                                                             <  1:1296353/100‑1 (MQ=255)
                  cgcCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAAc                                                                             >  2:4351453/1‑100 (MQ=255)
                    ccTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTg                                                                           <  1:1459532/100‑1 (MQ=255)
                     cTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTgg                                                                          >  2:1632395/1‑100 (MQ=255)
                     cTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTgg                                                                          <  1:2825985/100‑1 (MQ=255)
                       gACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTACCTGGAt                                                                        >  2:3858869/1‑100 (MQ=255)
                       gACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGAt                                                                        >  2:2531450/1‑100 (MQ=255)
                       gACCTTCTGCGTTGGCGCTAAACCCACCACCGCTGGCGATCACATTGACGCGTGGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGAt                                                                        <  2:4861605/100‑1 (MQ=255)
                        aCCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATa                                                                       >  1:1502816/1‑100 (MQ=255)
                         ccTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAg                                                                      >  2:3754533/1‑100 (MQ=255)
                          cTTCTGCGTTGGCGCGATACGCATCACCGCTGGCGATAACATTGACGCGTTGTCCGGCTTATACACGCCATACTTGGCGAAGCTTGGTTAACTGGATAgg                                                                     >  2:442401/1‑100 (MQ=255)
                           ttCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGt                                                                    >  2:4160276/1‑100 (MQ=255)
                           ttCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGt                                                                    <  1:1608514/100‑1 (MQ=255)
                           ttCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGt                                                                    <  2:271978/100‑1 (MQ=255)
                            tcggCGTTGGCGCTAACCCCATCACCCCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGtt                                                                   <  2:2515368/97‑1 (MQ=255)
                             cTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTg                                                                  >  2:3818601/1‑100 (MQ=255)
                              tGCGTTGGCGCTAAACCCATCACCGCTCACGATCACATTCACGCGTTGTCCCACTTTTACCCGCCATACCTGGCGAATCTGGGTTAACTGGATCGGTTGa                                                                 >  1:5421455/1‑100 (MQ=255)
                                tgTTGGCGCTAAACCCATCTCCGCGGGCGATCACATGGACGCTTGGTCCCGCTTGTCCCCGCCATACCTGGCGGAACTGGGTTAACTGGATAGGTTGATc                                                               <  1:4693074/99‑1 (MQ=255)
                                cGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATc                                                               <  2:834943/100‑1 (MQ=255)
                                cGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATc                                                               <  1:3020169/100‑1 (MQ=255)
                                 gTTGGCGCTAAACCTATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCg                                                              <  2:245949/100‑1 (MQ=255)
                                 gTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGGGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCg                                                              <  2:506381/100‑1 (MQ=255)
                                 gTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCg                                                              >  1:2336296/1‑100 (MQ=255)
                                   tGGCGCTAAACCCATCACCGCGGGCGATCACATTGTCGCGTTGGCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCggg                                                            <  2:3304028/100‑1 (MQ=255)
                                    ggCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGt                                                           <  1:3440702/100‑1 (MQ=255)
                                     gcgcTAACCCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTg                                                          <  1:1135658/100‑1 (MQ=255)
                                     gcgcTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTg                                                          <  2:4727378/100‑1 (MQ=255)
                                     gcgcTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTg                                                          >  2:2079806/1‑100 (MQ=255)
                                        cTAAACCCTGCCCCGGGGGCGAGCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATa                                                       <  1:4191762/100‑1 (MQ=255)
                                        cTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATa                                                       >  2:4878729/1‑100 (MQ=255)
                                          aaaCCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGa                                                     >  1:4353986/1‑100 (MQ=255)
                                            ccccATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATc                                                   <  1:5200434/99‑1 (MQ=255)
                                            aCCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATc                                                   >  2:4740170/1‑100 (MQ=255)
                                            aCCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATc                                                   >  1:3586519/1‑100 (MQ=255)
                                             cccATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATcg                                                  <  2:264351/100‑1 (MQ=255)
                                             cccATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATcg                                                  <  1:719030/100‑1 (MQ=255)
                                              ccATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATcgc                                                 >  2:3477452/1‑100 (MQ=255)
                                                aTCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCgcgc                                               >  2:3246904/1‑100 (MQ=255)
                                                   aCCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAgg                                            <  1:1978133/100‑1 (MQ=255)
                                                   aCCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAgg                                            >  1:2231351/1‑100 (MQ=255)
                                                         ggCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATg                                      <  1:2164819/100‑1 (MQ=255)
                                                          gCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATgg                                     <  1:4830595/100‑1 (MQ=255)
                                                           cGATCACATTGACGCGTTGACCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCACGCAGGCTAATGGc                                    >  2:5032485/1‑100 (MQ=255)
                                                            gATCACATTGGCGCGAGGTCCGGCTTTTCCCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCa                                   <  2:1994095/100‑1 (MQ=255)
                                                             aTCACATTGCCGCGTTGTCCCGCTTTTCCCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCAt                                  <  1:3894398/100‑1 (MQ=255)
                                                             aTCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCAt                                  >  2:2778129/1‑100 (MQ=255)
                                                             aTCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCAt                                  >  1:1630965/1‑100 (MQ=255)
                                                             aTCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCAt                                  <  2:2601148/100‑1 (MQ=255)
                                                              tCACATTGACGCGTTGTCCCGCTTTTACCCGTCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATc                                 <  2:3670944/100‑1 (MQ=255)
                                                              tCACATTGACGCGTTGTCCCGCTTTTACCCGCCATACATGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATc                                 <  1:4525763/100‑1 (MQ=255)
                                                               cacaTTGACGGGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATca                                <  2:4615837/100‑1 (MQ=255)
                                                               cacaTTGACGCGTTGTCCCGCGTTTTCCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATca                                <  2:1304042/100‑1 (MQ=255)
                                                                acaTTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATcaa                               >  1:3923336/1‑100 (MQ=255)
                                                                 caTTGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATcaac                              >  2:2332143/1‑100 (MQ=255)
                                                                   ttGACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATcaacaa                            >  1:4263940/1‑100 (MQ=255)
                                                                    ggACGCGTTGTCCCGCTTTTACCCGCCATACCTCGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAg                           <  1:5387079/99‑1 (MQ=255)
                                                                     gACGCGTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGt                          >  2:4109138/1‑100 (MQ=255)
                                                                         cgTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGttgat                      <  1:441437/100‑1 (MQ=255)
                                                                         cgTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGttgat                      <  2:1104703/100‑1 (MQ=255)
                                                                          gTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGTTCAAGTGATCGATCGCTCAGGCTAATGGCATCAACAAGttgatt                     >  2:5245461/1‑100 (MQ=255)
                                                                          gTTGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGttgatt                     >  1:1706975/1‑100 (MQ=255)
                                                                           ttGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGttgattg                    <  2:49858/100‑1 (MQ=255)
                                                                           ttGTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGttgattg                    <  2:3152952/100‑1 (MQ=255)
                                                                             gTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTACTGGCATCAACAAGTtgattgat                  >  2:4337456/1‑100 (MQ=255)
                                                                             gTCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTtgattgat                  >  2:5448483/1‑100 (MQ=255)
                                                                              tCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGata                 <  2:1118714/100‑1 (MQ=255)
                                                                              tCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGata                 >  1:2086590/1‑100 (MQ=255)
                                                                              tCCCGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGata                 <  1:3852654/100‑1 (MQ=255)
                                                                               cccGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGatat                >  1:2724267/1‑100 (MQ=255)
                                                                               cccGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGatat                >  1:2557454/1‑100 (MQ=255)
                                                                               cccGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGAACGCGCAGGCTAATGGCATCAACAAGTTGATTGatat                >  1:5249291/1‑100 (MQ=255)
                                                                                ccGCTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATc               >  2:1854761/1‑100 (MQ=255)
                                                                                   cTTTTACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAg            >  1:324833/1‑100 (MQ=255)
                                                                                       tACCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAGCAcc        >  2:4038695/1‑100 (MQ=255)
                                                                                        aCCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGTATCAACAAGTTGATTGATATCCAGCACCg       <  2:3975888/100‑1 (MQ=255)
                                                                                        aCCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAGCACCg       >  2:2016858/1‑100 (MQ=255)
                                                                                        aCCCGCCATACCTGGCGAAACTGGGTTAACTGGATAGCTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAGCACCg       >  2:5366718/1‑100 (MQ=255)
                                                                                         cccGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGGTCGCGCAGGCTAATGGCATCAGCAAGTTGATTGATATCCAGCACCGt      <  1:4977759/100‑1 (MQ=255)
                                                                                         cccGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAGCACCGt      >  1:3319517/1‑100 (MQ=255)
                                                                                          ccGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAGCACCGTa     >  2:2250250/1‑100 (MQ=255)
                                                                                           cGCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAGCACCGTAc    >  1:2490981/1‑100 (MQ=255)
                                                                                            gCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAGCACCGTAcg   <  2:4476608/100‑1 (MQ=255)
                                                                                            gCCATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAGCACCGTAcg   <  1:447310/100‑1 (MQ=255)
                                                                                             ccATACCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTTATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAGCACCGTAcgc  <  2:804767/100‑1 (MQ=255)
                                                                                                 |                                                                                               
GACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGGCGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAGCACCGTACGC  >  NC_000913/1130291‑1130483

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: