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,4690GA99.6% 784.3 / NA 223T354T (ACG→ACAfliKflagellar hook‑length control protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (92/130);  minor base C (0/1);  total (92/131)
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 228 total aligned reads displayed.

TTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATC  >  NC_000913/2017380‑2017557
                                                                                         |                                                                                        
tttAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACtt                                                                                       <  2:3182193/93‑1 (MQ=255)
   aGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGACCCTCTGGCGGGAGAAGACGACGATACCCTTCCg                                                                                    <  2:1738381/93‑1 (MQ=255)
     tGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGt                                                                                  <  2:4267870/93‑1 (MQ=255)
     tGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGt                                                                                  <  2:335746/93‑1 (MQ=255)
      ggTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGtt                                                                                 <  2:4185367/93‑1 (MQ=255)
       gTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTc                                                                                <  1:4263794/93‑1 (MQ=255)
       gTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGGTc                                                                                >  2:2663497/1‑93 (MQ=255)
        tCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTcc                                                                               <  1:4449092/93‑1 (MQ=255)
        tCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTcc                                                                               <  1:5439710/93‑1 (MQ=255)
         cagcagcagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTccc                                                                              <  1:3740284/93‑1 (MQ=255)
         cagcagcagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTccc                                                                              <  1:588779/93‑1 (MQ=255)
          agcagcagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCg                                                                             >  2:4481290/1‑93 (MQ=255)
          agcagcagGCCGCTCCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGACCCTCGGGCGGGGGAAGACGACGATACACTTCCGGTTCCCg                                                                             <  2:3547349/93‑1 (MQ=255)
           gcagcagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGt                                                                            >  1:2753193/1‑93 (MQ=255)
           gcagcagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGt                                                                            <  2:4927899/93‑1 (MQ=255)
            cagcagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGtc                                                                           <  2:4561155/93‑1 (MQ=255)
             agcagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGATGATACACTTCCGGTTCCCGtct                                                                          >  1:4244933/1‑93 (MQ=255)
              gcagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGtctc                                                                         >  2:1211884/1‑93 (MQ=255)
              gcagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGtctc                                                                         >  1:5060903/1‑93 (MQ=255)
              gcagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGtctc                                                                         >  1:1142669/1‑93 (MQ=255)
               cagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCGCGTctct                                                                        >  1:2204701/1‑93 (MQ=255)
                agGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCtt                                                                       <  2:1147469/93‑1 (MQ=255)
                agGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCtt                                                                       >  1:1523876/1‑93 (MQ=255)
                 ggCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCttt                                                                      >  2:5431658/1‑93 (MQ=255)
                      cTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAg                                                                 >  2:1767356/1‑93 (MQ=255)
                      cTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAg                                                                 <  1:1894587/93‑1 (MQ=255)
                      cTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAg                                                                 >  1:4124686/1‑93 (MQ=255)
                        tCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAggg                                                               <  2:3187857/93‑1 (MQ=255)
                        tCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAggg                                                               <  1:944789/93‑1 (MQ=255)
                        tCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAggg                                                               >  1:5532461/1‑93 (MQ=255)
                           cagcaacagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCgt                                                            >  1:2318687/1‑93 (MQ=255)
                           cagcaacagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCgt                                                            <  1:3092926/93‑1 (MQ=255)
                            agcaacagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCgtg                                                           <  1:5059751/93‑1 (MQ=255)
                            agcaacagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCgtg                                                           <  1:5258399/93‑1 (MQ=255)
                             gcaacagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCgtgt                                                          <  1:5404176/93‑1 (MQ=255)
                              caacagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTa                                                         >  2:4279340/1‑93 (MQ=255)
                                acagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAAc                                                       >  2:2508261/1‑93 (MQ=255)
                                acagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAAc                                                       >  2:1911635/1‑93 (MQ=255)
                                   gcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAgg                                                    >  2:5200914/1‑93 (MQ=255)
                                   gcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAgg                                                    <  1:5814400/93‑1 (MQ=255)
                                    caaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGc                                                   >  2:3854232/1‑93 (MQ=255)
                                     aaaGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCa                                                  >  2:2276503/1‑93 (MQ=255)
                                      aaGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTCCAAGGGCGTGTAACAGGCaa                                                 >  2:3927281/1‑93 (MQ=255)
                                       aGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAAc                                                <  1:3938162/93‑1 (MQ=255)
                                       aGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAAc                                                <  1:1969239/93‑1 (MQ=255)
                                       aGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAAc                                                <  1:719346/93‑1 (MQ=255)
                                         ccAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAg                                              <  2:1043760/93‑1 (MQ=255)
                                          cAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAgc                                             <  2:4773209/93‑1 (MQ=255)
                                          cAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAgc                                             <  1:4658610/93‑1 (MQ=255)
                                           aaCGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAgcg                                            <  2:1677982/93‑1 (MQ=255)
                                            aCGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAgcgg                                           <  2:545699/93‑1 (MQ=255)
                                            aCGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGGGTAATAGGCAACAgcgg                                           >  1:2471048/1‑93 (MQ=255)
                                              gCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACCAGGGCGTGTAACAGGCAACAgcggcg                                         >  1:3801246/1‑93 (MQ=255)
                                                acaGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTTACAGGCAACAgcggcgct                                       >  2:2291469/1‑91 (MQ=255)
                                                 caGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTg                                      <  2:247124/93‑1 (MQ=255)
                                                  aGCAAACCATGATCCTCTGGCGGGGGAAGACGACGATACACTGCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGa                                     <  1:4352578/93‑1 (MQ=255)
                                                  aGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGa                                     >  2:5787761/1‑93 (MQ=255)
                                                   gCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGat                                    <  2:2739607/93‑1 (MQ=255)
                                                   gCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGat                                    <  1:4491778/93‑1 (MQ=255)
                                                    cAAACCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGata                                   <  1:4758079/93‑1 (MQ=255)
                                                    cAAACCATGAACCTCTGGCGGGGGAAGAAGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGata                                   <  2:4603733/93‑1 (MQ=255)
                                                      aaCCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATAtt                                 >  2:3115318/1‑93 (MQ=255)
                                                      aaCCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATAtt                                 >  1:2164837/1‑93 (MQ=255)
                                                      aaCCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATAtt                                 >  1:1075530/1‑93 (MQ=255)
                                                      aaCCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATAtt                                 >  1:4037639/1‑93 (MQ=255)
                                                       aCCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATAttt                                >  2:3286706/1‑93 (MQ=255)
                                                       aCCATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTATCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATAttt                                >  1:3613655/1‑93 (MQ=255)
                                                        ccATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATAtttt                               <  1:5978675/93‑1 (MQ=255)
                                                        ccATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATAtttt                               >  1:5649784/1‑93 (MQ=255)
                                                         cATGAGCCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTc                              <  1:3025796/93‑1 (MQ=255)
                                                         cATGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTc                              <  1:1405268/93‑1 (MQ=255)
                                                          aTGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCg                             <  2:6119537/93‑1 (MQ=255)
                                                           tGAACCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGc                            <  1:5192677/93‑1 (MQ=255)
                                                            gAACCTCTGGCGGGGTAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGcc                           <  2:2782895/93‑1 (MQ=255)
                                                              aCCTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTa                         >  2:3050366/1‑93 (MQ=255)
                                                               ccTCTGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTaa                        >  1:3704376/1‑93 (MQ=255)
                                                                 cctGGCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACg                      <  1:17348/92‑1 (MQ=255)
                                                                     gCGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCag                  <  1:2575356/93‑1 (MQ=255)
                                                                      cGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCaga                 >  2:4040586/1‑93 (MQ=255)
                                                                      cGGGGGAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCaga                 <  1:1696533/93‑1 (MQ=255)
                                                                       gggggAAGCCGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCagag                <  2:1574544/93‑1 (MQ=255)
                                                                       gggggAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCagag                <  2:6185934/93‑1 (MQ=255)
                                                                        ggggAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAgg               >  2:90222/1‑93 (MQ=255)
                                                                        ggggAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAgg               <  1:2643294/93‑1 (MQ=255)
                                                                        ggggAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAgg               <  1:3220718/93‑1 (MQ=255)
                                                                        ggggAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAgg               <  2:5002799/93‑1 (MQ=255)
                                                                        ggggAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAgg               >  2:3916012/1‑93 (MQ=255)
                                                                          ggAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTa             <  1:5934344/93‑1 (MQ=255)
                                                                          ggAAGACGACGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTa             <  2:1727817/93‑1 (MQ=255)
                                                                              gacgacGATACACTTCCGGTTCTCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCAc         <  1:5656977/93‑1 (MQ=255)
                                                                              gacgacGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCAc         <  1:6138146/93‑1 (MQ=255)
                                                                              gacgacGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCAc         <  2:3688154/93‑1 (MQ=255)
                                                                              gacgacGATACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCAc         >  2:6294683/1‑93 (MQ=255)
                                                                                cgacgaTACACTGCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCg       >  1:3405517/1‑93 (MQ=255)
                                                                                   cgaTACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTaa    <  2:4429310/93‑1 (MQ=255)
                                                                                   cgaTACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTaa    <  1:2098507/93‑1 (MQ=255)
                                                                                    gaTACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAAt   >  2:5289424/1‑93 (MQ=255)
                                                                                     aTACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATc  <  1:664618/93‑1 (MQ=255)
                                                                                     aTACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATc  <  2:5152872/93‑1 (MQ=255)
                                                                                     aTACACTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATc  <  2:1874943/93‑1 (MQ=255)
                                                                                         |                                                                                        
TTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATC  >  NC_000913/2017380‑2017557

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