Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,017,469 G→A T354T (ACG→ACA fliK → flagellar hook‑length control protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,017,4690GA98.4% 436.7 / ‑3.4 128T354T (ACG→ACAfliKflagellar hook‑length control protein
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (59/67);  total (59/69)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.98e-01

Only 100 of 139 total aligned reads displayed.

CTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCC  >  NC_000913/2017379‑2017558
                                                                                          |                                                                                         
cTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACt                                                                                         >  2:4263750/1‑93 (MQ=255)
 tttAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGATGACGATACACtc                                                                                        >  2:4216992/1‑92 (MQ=255)
 tttAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACtt                                                                                        >  2:1212779/1‑93 (MQ=255)
 tttAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACtt                                                                                        >  1:1646615/1‑93 (MQ=255)
  ttAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTc                                                                                       <  2:496101/93‑1 (MQ=255)
   tAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTcc                                                                                      <  1:480498/93‑1 (MQ=255)
   tAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTcc                                                                                      <  1:3789432/93‑1 (MQ=255)
   tAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTcc                                                                                      <  2:4581462/93‑1 (MQ=255)
    aGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATTAACCTCTGGCGGGGGAAGACGACGATACACTTCCg                                                                                     >  2:339113/1‑93 (MQ=255)
       ggTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGtt                                                                                  <  2:1274851/93‑1 (MQ=255)
          cagcagcagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTccc                                                                               <  2:3686021/93‑1 (MQ=255)
            gcagcagGCCGCTTCCCAGCAACAGCAAATCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGt                                                                             >  1:4731432/1‑93 (MQ=255)
            gcagcagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGt                                                                             >  1:2487941/1‑93 (MQ=255)
              agcagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGtct                                                                           >  2:261875/1‑93 (MQ=255)
               gcagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGCACCTCTGGCGGGGGACGACGACGATCCGCTTGCGCTTCCCGtctc                                                                          >  2:2831346/1‑93 (MQ=255)
               gcagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGACCCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGtctc                                                                          <  2:2147342/93‑1 (MQ=255)
                 agGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCtt                                                                        >  2:70336/1‑93 (MQ=255)
                 agGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCtt                                                                        <  2:1599433/93‑1 (MQ=255)
                  ggCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCttt                                                                       >  2:1583023/1‑93 (MQ=255)
                  ggCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCttt                                                                       >  2:3123187/1‑93 (MQ=255)
                     cGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACa                                                                    >  1:1168398/1‑93 (MQ=255)
                     cGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACa                                                                    >  2:4024257/1‑93 (MQ=255)
                      gCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACaa                                                                   >  2:1488740/1‑93 (MQ=255)
                      gCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACaa                                                                   >  2:1488741/1‑93 (MQ=255)
                       cTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAg                                                                  <  2:3201949/93‑1 (MQ=255)
                       cTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAg                                                                  <  2:736911/93‑1 (MQ=255)
                         tCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCGTCTCTTTACAAggg                                                                <  1:3993350/93‑1 (MQ=255)
                         tCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAggg                                                                >  1:557663/1‑93 (MQ=255)
                         tCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAggg                                                                <  1:1271554/93‑1 (MQ=255)
                            cagcaacagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGTGCgt                                                             <  2:4053568/93‑1 (MQ=255)
                              gcaacagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCgtgt                                                           >  2:1819443/1‑93 (MQ=255)
                              gcaacagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCgtgt                                                           <  1:3876315/93‑1 (MQ=255)
                              gcaacagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCgtgt                                                           <  2:2570622/93‑1 (MQ=255)
                                aacagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTaa                                                         <  2:3799703/93‑1 (MQ=255)
                                aacagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTaa                                                         >  1:3203752/1‑93 (MQ=255)
                                aacagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTaa                                                         >  1:4552281/1‑93 (MQ=255)
                                aacagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTaa                                                         <  1:1012518/93‑1 (MQ=255)
                                 acagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAAc                                                        <  1:1381917/93‑1 (MQ=255)
                                    gcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAgg                                                     >  1:316285/1‑93 (MQ=255)
                                     caaAGCTAATGCACAGCAAACCATGAACCTCTAGCGGGGGAAGACGACAATAAACTTCCTGTTCCTGTGTTTTTACAAGGACGTGCAACAGGc                                                    >  2:3945442/1‑93 (MQ=255)
                                     caaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGc                                                    <  2:3892511/93‑1 (MQ=255)
                                       aaGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCaa                                                  <  2:2448629/93‑1 (MQ=255)
                                        aGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAAc                                                 >  2:3893389/1‑93 (MQ=255)
                                        aGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAAc                                                 <  2:785845/93‑1 (MQ=255)
                                         gCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACa                                                <  2:4374196/93‑1 (MQ=255)
                                           cAACGCACAGCAAACCATGAGCCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAgc                                              <  1:1484110/93‑1 (MQ=255)
                                           cAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAgc                                              <  1:28969/93‑1 (MQ=255)
                                            aaCGCACAGCAAACCATGACCCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAgcg                                             <  2:3368014/93‑1 (MQ=255)
                                               gCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAgcggcg                                          >  1:4463536/1‑93 (MQ=255)
                                               gCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAgcggcg                                          <  2:2150161/93‑1 (MQ=255)
                                                cacaGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGt                                         >  1:269521/1‑93 (MQ=255)
                                                 acaGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGtt                                        >  1:1945910/1‑93 (MQ=255)
                                                   aGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGa                                      <  1:4185416/93‑1 (MQ=255)
                                                    gCAAACCATGAACCTCTGGCGGGTGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGat                                     <  1:3528981/93‑1 (MQ=255)
                                                     cAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGata                                    <  2:261350/93‑1 (MQ=255)
                                                      aaaCCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGatat                                   >  1:3950463/1‑93 (MQ=255)
                                                       aaCCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATAtt                                  >  2:723827/1‑93 (MQ=255)
                                                       aaCCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATAtt                                  >  1:1530638/1‑93 (MQ=255)
                                                       aaCCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATAtt                                  >  2:1121396/1‑93 (MQ=255)
                                                        aCCATGAACGTCTGGCGGGGGAAGACGGCGATCCACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATAttt                                 <  1:2229887/93‑1 (MQ=255)
                                                        aCCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATAttt                                 >  2:3782493/1‑93 (MQ=255)
                                                         ccATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATAtttt                                <  1:4128007/93‑1 (MQ=255)
                                                         ccATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATAtttt                                <  1:4319676/93‑1 (MQ=255)
                                                         ccATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATAtttt                                <  1:3061152/93‑1 (MQ=255)
                                                         ccATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATAtttt                                <  2:519907/93‑1 (MQ=255)
                                                          cATGAACCTCTGGTGGGGGAAGACGACGATACACTTCAGGTTCCCGTCGCTTTACAAGGGCGTGTAACAGGCAACAGCGGCTTTGATATTTTc                               <  2:2490189/93‑1 (MQ=255)
                                                          cATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTc                               >  2:1567677/1‑93 (MQ=255)
                                                          cATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTc                               >  1:699873/1‑93 (MQ=255)
                                                           aTGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCg                              >  1:3658885/1‑93 (MQ=255)
                                                           aTGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCg                              >  2:3593853/1‑93 (MQ=255)
                                                            tGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGc                             >  1:3662426/1‑93 (MQ=255)
                                                               aCCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTa                          <  1:2684490/93‑1 (MQ=255)
                                                                 ctctGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAAc                        <  2:3806337/93‑1 (MQ=255)
                                                                 ctctGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAAc                        <  1:4576439/93‑1 (MQ=255)
                                                                 ctctGGCGGGGGAAGACGACGATACACTTCCGGGTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAAc                        <  1:4576437/93‑1 (MQ=255)
                                                                    tGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTc                     <  2:2808689/93‑1 (MQ=255)
                                                                       cGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCaga                  <  2:2442860/93‑1 (MQ=255)
                                                                       cGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCaga                  <  1:2199777/93‑1 (MQ=255)
                                                                        gggggAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCagag                 >  1:4085676/1‑93 (MQ=255)
                                                                        gggggAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCagag                 <  1:137903/93‑1 (MQ=255)
                                                                        gggggAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCagag                 >  1:95908/1‑93 (MQ=255)
                                                                        gggggAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCagag                 >  1:1461509/1‑93 (MQ=255)
                                                                        gggggAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCagag                 >  1:369239/1‑93 (MQ=255)
                                                                        ggggaAAGACGATGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGTCGTTGATATTTTCGCCTAACGTCagag                 <  2:580864/93‑1 (MQ=255)
                                                                         ggggAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAgg                <  1:2779651/93‑1 (MQ=255)
                                                                           ggAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTa              <  2:4449251/93‑1 (MQ=255)
                                                                           ggAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTa              >  1:1867493/1‑93 (MQ=255)
                                                                            gAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAg             >  1:1036865/1‑93 (MQ=255)
                                                                               gacgacGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCAc          <  2:3680656/93‑1 (MQ=255)
                                                                               gacgacGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCAc          <  1:3074780/93‑1 (MQ=255)
                                                                                 cgacgaTACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCg        >  1:2249298/1‑93 (MQ=255)
                                                                                 cgacgaTACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCg        <  2:3150720/93‑1 (MQ=255)
                                                                                 cgacgaTACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCg        <  2:1813011/93‑1 (MQ=255)
                                                                                   acgaTACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTa      >  1:2463078/1‑93 (MQ=255)
                                                                                    cgaTACACTTCCGGCTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTaa     >  2:2659955/1‑93 (MQ=255)
                                                                                     gaTACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAAt    >  2:2311829/1‑93 (MQ=255)
                                                                                     gaTACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAAt    <  1:643608/93‑1 (MQ=255)
                                                                                       tACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATcc  >  1:1557372/1‑93 (MQ=255)
                                                                                       tACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATcc  <  2:1887699/93‑1 (MQ=255)
                                                                                       tACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATcc  <  1:909799/93‑1 (MQ=255)
                                                                                          |                                                                                         
CTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCC  >  NC_000913/2017379‑2017558

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