Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,133,157 A→T 7.0% D195V (GAC→GTC)  flgE → flagellar hook protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,133,1570AT 7.0% 313.9 / 11.3 86D195V (GAC→GTC) flgEflagellar hook protein
Reads supporting (aligned to +/- strand):  ref base A (40/40);  new base T (3/3);  total (43/43)
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.10e-01

CCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGATGG  >  NC_000913/1133019‑1133297
                                                                                                                                          |                                                                                                                                            
ccaccaCGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTcatgaca                                                                                                                                            <  1:196591/141‑1 (MQ=255)
ccaccaCGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTcatgaca                                                                                                                                            >  2:196591/1‑141 (MQ=255)
              aTGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGa                                                                                                                       >  3:888764/1‑148 (MQ=255)
               tGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCt                                                                                                                                   <  1:1014726/135‑1 (MQ=255)
               tGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCt                                                                                                                                   >  2:1014726/1‑135 (MQ=255)
                 cAGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCg                                                                                                                    <  4:888764/148‑1 (MQ=255)
                   gATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCg                                                                                                                                      >  1:1077826/1‑128 (MQ=255)
                   gATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCg                                                                                                                                      <  2:1077826/128‑1 (MQ=255)
                           tGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTa                                                                                                                                  >  4:287184/1‑124 (MQ=255)
                           tGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTa                                                                                                                                  <  3:287184/124‑1 (MQ=255)
                             aaTTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTtg                                                                                                                             <  3:219334/127‑1 (MQ=255)
                             aaTTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGATTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTtg                                                                                                                             >  4:219334/1‑127 (MQ=255)
                                tCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCTGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGt                                                                                                                                     >  1:731800/1‑116 (MQ=255)
                                tCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGTCATGAGcatt                                                                                                                                    <  3:100097/117‑4 (MQ=255)
                                tCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAg                                                                                                                                        >  3:472433/1‑113 (MQ=255)
                                tCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAg                                                                                                                                        <  4:472433/113‑1 (MQ=255)
                                tCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGt                                                                                                                                     <  2:731800/116‑1 (MQ=255)
                                tCCAGTGATCAGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATTCTCATGTCATGAGcatt                                                                                                                                    >  4:100097/1‑114 (MQ=255)
                                 ccAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCa                                                                                                    >  2:762462/1‑148 (MQ=255)
                                   aGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAgg                                                                                                  >  4:753237/1‑148 (MQ=255)
                                    gTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGt                                                                                                 >  2:700886/1‑148 (MQ=255)
                                         ccGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTtg                                                                                                                             >  3:855187/1‑115 (MQ=255)
                                         ccGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTtg                                                                                                                             <  4:855187/115‑1 (MQ=255)
                                          cGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGTCATGAGca                                                                                                                                      >  1:58578/1‑104 (MQ=255)
                                          cGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGTCATGAGca                                                                                                                                      <  2:58578/105‑2 (MQ=255)
                                                 tACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTtg                                                                                                                             <  1:180672/107‑1 (MQ=255)
                                                 tACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTtg                                                                                                                             >  2:180672/1‑107 (MQ=255)
                                                  aCTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATa                                                                                   <  3:753237/148‑1 (MQ=255)
                                                    tGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGc                                                                                 <  1:700886/148‑1 (MQ=255)
                                                      ttACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATa                                                                                   >  2:451502/1‑144 (MQ=255)
                                                      ttACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATa                                                                                   <  1:451502/144‑1 (MQ=255)
                                                              ttCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCaa                                                                       <  1:611111/148‑1 (MQ=255)
                                                                      cAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGc                                                                                 <  3:422995/130‑1 (MQ=255)
                                                                      cAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGc                                                                                 >  4:422995/1‑130 (MQ=255)
                                                                          aaTGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGc                                                                                 >  4:688437/1‑126 (MQ=255)
                                                                          aaTGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGc                                                                                 <  3:688437/126‑1 (MQ=255)
                                                                            tGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAg                                                         >  1:16270/1‑148 (MQ=255)
                                                                                gATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAg                                                     <  1:762462/148‑1 (MQ=255)
                                                                                      tatgACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTTTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAg                                                         >  4:837976/1‑138 (MQ=255)
                                                                                      tataACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCg                                                                                                                    >  3:975355/1‑79 (MQ=255)
                                                                                      tataACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCg                                                                                                                    <  4:975355/79‑1 (MQ=255)
                                                                                      tataACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAg                                                         <  3:837976/138‑1 (MQ=255)
                                                                                       ataACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATGACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGacaac                                              >  3:352213/1‑148 (MQ=255)
                                                                                        taACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGacaaca                                             <  2:16270/148‑1 (MQ=255)
                                                                                         aaCAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGa                                                           <  2:155186/133‑1 (MQ=255)
                                                                                         aaCAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGa                                                           >  1:155186/1‑133 (MQ=255)
                                                                                             aaaaaGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCa                                                                                        <  4:805906/100‑1 (MQ=255)
                                                                                             aaaaaGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCa                                                                                        >  3:805906/1‑100 (MQ=255)
                                                                                              aaaaGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGATATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCg                                                   <  1:813057/136‑1 (MQ=255)
                                                                                              aaaaGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGATATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCg                                                   >  2:813057/1‑136 (MQ=255)
                                                                                                aaGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAgg                                                                                      <  1:190755/99‑1 (MQ=255)
                                                                                                aaGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAgg                                                                                      >  2:190755/1‑99 (MQ=255)
                                                                                                 aGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGc                                                                  >  1:159281/1‑118 (MQ=255)
                                                                                                 aGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGc                                                                  <  2:159281/118‑1 (MQ=255)
                                                                                                    ttCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGa                                        <  4:677415/141‑1 (MQ=255)
                                                                                                    ttCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGa                                        >  3:677415/1‑141 (MQ=255)
                                                                                                     tCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGa                                                                            >  3:1031255/1‑104 (MQ=255)
                                                                                                     tCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGa                                                                            <  4:1031255/104‑1 (MQ=255)
                                                                                                      cGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCa                                                                                                    >  4:684235/1‑79 (MQ=255)
                                                                                                      cGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCa                                                                                                    <  3:684235/79‑1 (MQ=255)
                                                                                                        gTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCa                                                                                        <  1:43849/89‑1 (MQ=255)
                                                                                                        gTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCa                                                                                        >  2:43849/1‑89 (MQ=255)
                                                                                                        gTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGacaac                                              >  3:498901/1‑131 (MQ=255)
                                                                                                        gTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGacaac                                              <  4:498901/131‑1 (MQ=255)
                                                                                                          gACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAAc                                                                                                         <  4:299281/70‑1 (MQ=255)
                                                                                                          gACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAAc                                                                                                         >  3:299281/1‑70 (MQ=255)
                                                                                                           aCTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCg                                                                                                                    >  4:28217/1‑58 (MQ=255)
                                                                                                           aCTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCg                                                                                                                    <  3:28217/58‑1 (MQ=255)
                                                                                                            cTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATgg                         >  1:895501/1‑148 (MQ=255)
                                                                                                              gTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGAt                                                                           <  1:574655/96‑1 (MQ=255)
                                                                                                              gTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGAt                                                                           >  2:574655/1‑96 (MQ=255)
                                                                                                                 ttCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTgg                                                                                                      >  1:363626/1‑66 (MQ=255)
                                                                                                                 ttCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTgg                                                                                                      <  2:363626/66‑1 (MQ=255)
                                                                                                                  tCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGaca                                                >  1:763296/1‑119 (MQ=255)
                                                                                                                  tCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGaca                                                <  2:763296/119‑1 (MQ=255)
                                                                                                                     aCAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATGACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGt                <  4:352213/148‑1 (MQ=255)
                                                                                                                     aCAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCaa                                                                       >  1:776465/1‑93 (MQ=255)
                                                                                                                     aCAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCaa                                                                       <  2:776465/93‑1 (MQ=255)
                                                                                                                          cAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGcaca                     <  3:889487/138‑1 (MQ=255)
                                                                                                                          cAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGcaca                     >  4:889487/1‑138 (MQ=255)
                                                                                                                            gggTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTa                                   <  3:53962/122‑1 (MQ=255)
                                                                                                                            gggTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTa                                   >  4:53962/1‑122 (MQ=255)
                                                                                                                              gTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGaca                                                <  3:564726/107‑1 (MQ=255)
                                                                                                                              gTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGaca                                                >  4:564726/1‑107 (MQ=255)
                                                                                                                                aaTGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGa     >  3:41064/1‑148 (MQ=255)
                                                                                                                                 acgCTCATGTGATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAtt                                                               >  1:592739/3‑89 (MQ=255)
                                                                                                                                 acgCTCATGTCATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCAtt                                                               <  2:592739/87‑1 (MQ=255)
                                                                                                                                   gCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGATgg  <  4:41064/148‑1 (MQ=255)
                                                                                                                                      catgacatgaGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGacaaca                                             <  1:810504/102‑1 (MQ=255)
                                                                                                                                      catgacatgaGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGacaaca                                             >  2:810504/1‑102 (MQ=255)
                                                                                                                                      catgacatgaGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTa                             >  2:2124/1‑118 (MQ=255)
                                                                                                                                      catgacatgaGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTa                             <  1:2124/118‑1 (MQ=255)
                                                                                                                                      catgacatgaGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGa     <  4:853166/142‑1 (MQ=255)
                                                                                                                                      catgacatgaGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGa     >  3:853166/1‑142 (MQ=255)
                                                                                                                                        tgacatgaGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGAt            >  4:576645/1‑133 (MQ=255)
                                                                                                                                        tgacatgaGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGAt            <  3:576645/133‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                            
CCACCACGGCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGATGG  >  NC_000913/1133019‑1133297

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