Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,966,374 C→T 5.1% intergenic (‑183/+19) flhB ← / ← cheZ flagellin export apparatus, substrate specificity protein/chemotaxis regulator, protein phosphatase for CheY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,966,3740CT 5.1% 387.1 / 5.4 118intergenic (‑183/+19)flhB/cheZflagellin export apparatus, substrate specificity protein/chemotaxis regulator, protein phosphatase for CheY
Reads supporting (aligned to +/- strand):  ref base C (60/52);  new base T (3/3);  total (63/55)
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.51e-01

Only 100 of 122 total aligned reads displayed.

TGATGCCAGAGCGCAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGACGCGACTCCTGTTCCG  >  NC_000913/1966225‑1966523
                                                                                                                                                     |                                                                                                                                                     
tgatgCCAGAGCGCAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGc                                                                                                                                                       >  1:549379/1‑150 (MQ=255)
 gatgCCAGAGCGCAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGcc                                                                                                                                                      >  1:1161169/1‑150 (MQ=255)
    gCCAGAGCGCAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCcaca                                                                                                                                                   <  2:1934556/150‑1 (MQ=255)
    gCCAGAGCGCAAAGCGTGCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCcaca                                                                                                                                                   <  2:1440034/150‑1 (MQ=255)
     ccAGAGCGCAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCcacaa                                                                                                                                                  >  2:1612959/1‑149 (MQ=255)
       agagCGCAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCa                                                                                                                                                >  1:1415608/1‑150 (MQ=255)
        gagCGCAAAGCGTTCCATCGTTTGAGTAAGGGGCAAACCAGGCGGGATTTAGGGCTTTTGCGGCGACACCTCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTCTCAGACCGCCTGATATGACGTGGTCACGCCACATCAg                                                                                                                                               <  1:198316/150‑1 (MQ=255)
          gcgcAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGc                                                                                                                                             <  1:935969/150‑1 (MQ=255)
           cacAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATGTGGCGTGACCACGTCATATCAGGcg                                                                                                                                            <  2:808912/148‑2 (MQ=255)
           agcAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTATGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGAGGCGACATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCa                                                                                                                                            >  1:390839/2‑150 (MQ=255)
             cAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAAt                                                                                                                                          >  1:192518/1‑150 (MQ=255)
              aaaGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATa                                                                                                                                         <  1:710994/150‑1 (MQ=255)
              aaaGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGGCATCGCCTTATCAGACCGCCGGATGTGGCGTGACCACGTCATATCAGGcggtc                                                                                                                                         <  2:1543613/150‑5 (MQ=255)
               aaGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATaa                                                                                                                                        >  1:1954093/1‑149 (MQ=255)
               aaGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATAc                                                                                                                                        >  1:27878/1‑150 (MQ=255)
                aGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACa                                                                                                                                       >  1:215649/1‑150 (MQ=255)
                   gTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAAt                                                                                                                                    <  1:1757855/150‑1 (MQ=255)
                    ttCAATGGTTTGAGTAAGTTTATTTACATGAGGGATTTAGGGCTTTTGGTGCCACACATCAAGCTTAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCATAACGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATc                                                                                                                                   <  2:75532/150‑1 (MQ=255)
                          ggTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATc                                                                                                                             >  2:349653/1‑150 (MQ=255)
                              tGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAg                                                                                                                         <  1:1198001/150‑1 (MQ=255)
                              tGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAg                                                                                                                         <  1:149539/150‑1 (MQ=255)
                              tGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAg                                                                                                                         <  2:1003464/150‑1 (MQ=255)
                                   aaGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTAt                                                                                                                    >  1:606990/1‑150 (MQ=255)
                                      gggCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCca                                                                                                                 >  1:1529980/1‑150 (MQ=255)
                                         cAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCTTCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCGACATCAGTCAATACAAAACAAAATCCACGACTATCaaaca                                                                                                              >  2:1817417/1‑150 (MQ=255)
                                         cAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCcaaca                                                                                                              <  1:616796/150‑1 (MQ=255)
                                         cAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCcaaca                                                                                                              <  1:767017/150‑1 (MQ=255)
                                             aCAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATc                                                                                                          <  1:1389315/150‑1 (MQ=255)
                                              cAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCg                                                                                                         <  2:1870753/150‑1 (MQ=255)
                                               aGGCGGGAGTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGt                                                                                                        <  2:258942/150‑1 (MQ=255)
                                                ggCGGGATATAGGGCTTTTGCTGCCAGACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCGTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTc                                                                                                       <  2:100402/150‑1 (MQ=255)
                                                 gCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTcc                                                                                                      <  1:1931495/150‑1 (MQ=255)
                                                 gCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTcc                                                                                                      >  2:26084/1‑150 (MQ=255)
                                                     gATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCt                                                                                                  <  2:2007772/150‑1 (MQ=255)
                                                       tttAGGGCTTTTGCTGCCACACCTCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCCTCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCCAGCCTATCCATCAAATCGTCCACGTGa                                                                                                >  2:1396177/1‑150 (MQ=255)
                                                          aGGGCTTTTGCTGCCACACATCAAGCTTAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATcc                                                                                             <  1:727457/150‑1 (MQ=255)
                                                          aGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATcc                                                                                             <  2:1981058/150‑1 (MQ=255)
                                                          aGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATcc                                                                                             <  2:1378269/150‑1 (MQ=255)
                                                            ggCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACACATCGTCCACCTGATCCTg                                                                                           >  1:552213/1‑150 (MQ=255)
                                                                  tgctgcCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACATGATCCTGACTGGc                                                                                     >  2:1914151/1‑150 (MQ=255)
                                                                     tgcCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTAc                                                                                  <  2:1683140/150‑1 (MQ=255)
                                                                     tgcCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTAc                                                                                  >  2:713247/1‑150 (MQ=255)
                                                                      gcCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTAcc                                                                                 >  1:2017183/1‑150 (MQ=255)
                                                                       ccACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACca                                                                                >  2:1160089/1‑150 (MQ=255)
                                                                        cacacaTCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACcac                                                                               <  2:157099/150‑1 (MQ=255)
                                                                             aTCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCCCACCgg                                                                          >  2:60139/1‑150 (MQ=255)
                                                                              tCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGc                                                                         <  2:806882/150‑1 (MQ=255)
                                                                              tCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGc                                                                         >  1:1622455/1‑150 (MQ=255)
                                                                                 aGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCttt                                                                      >  1:1135229/1‑150 (MQ=255)
                                                                                  gCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTg                                                                     <  2:549379/150‑1 (MQ=255)
                                                                                       gggtgCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGAGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGt                                                                <  1:1244200/148‑1 (MQ=255)
                                                                                          tgCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACAACACCGGCTTTGCTGGTAta                                                             >  2:1301953/1‑149 (MQ=255)
                                                                                           gCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCg                                                            >  1:1607122/1‑150 (MQ=255)
                                                                                           gCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCg                                                            <  1:1073826/150‑1 (MQ=255)
                                                                                           gCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCg                                                            >  1:10728/1‑150 (MQ=255)
                                                                                            cGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCTCCTGATATGAAGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGa                                                           >  2:1432444/1‑150 (MQ=255)
                                                                                             gTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATACTGACTGGCTACCACACCGGCTTTGCGGGTATCGAc                                                          >  1:1055971/1‑150 (MQ=255)
                                                                                              tttGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCCCACCGGCTTTGCTGGTATCGAcc                                                         >  2:141825/1‑150 (MQ=255)
                                                                                              tttGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGAcc                                                         <  1:970228/150‑1 (MQ=255)
                                                                                              tttGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGAcc                                                         <  1:1581017/150‑1 (MQ=255)
                                                                                                  tCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAg                                                     <  2:27878/150‑1 (MQ=255)
                                                                                                    ggATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGGCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGt                                                   <  2:1455941/150‑1 (MQ=255)
                                                                                                        gcggcgTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCAt                                               <  1:1913114/150‑1 (MQ=255)
                                                                                                         cggcgTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCAtt                                              >  2:217934/1‑150 (MQ=255)
                                                                                                         cggcgTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCAtt                                              >  2:38486/1‑150 (MQ=255)
                                                                                                         cggcgTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCAACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCAtt                                              >  2:301705/1‑150 (MQ=255)
                                                                                                         cggcgTCATCGCCGTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAACTCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCAtt                                              <  2:588786/150‑1 (MQ=255)
                                                                                                          ggcgTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTa                                             >  1:485781/1‑150 (MQ=255)
                                                                                                          ggcgTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCAAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTa                                             >  1:593890/1‑150 (MQ=255)
                                                                                                            cgTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAg                                           >  1:731799/1‑150 (MQ=255)
                                                                                                               cATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGACCATTAAGCaa                                        >  1:258323/1‑150 (MQ=255)
                                                                                                                aTCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCAAACAAATCGTCCACCTGATCATGACTGGCTACAACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCaaa                                       >  1:1022931/1‑150 (MQ=255)
                                                                                                                    ccTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAActct                                   >  1:688241/1‑150 (MQ=255)
                                                                                                                      ttATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTgg                                 <  1:242530/150‑1 (MQ=255)
                                                                                                                       tATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGt                                <  2:342954/150‑1 (MQ=255)
                                                                                                                          cAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGtttt                             <  1:372464/150‑1 (MQ=255)
                                                                                                                            gACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGAATATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCa                           >  2:1462697/1‑150 (MQ=255)
                                                                                                                             attgCCTGATGTGGCGTGACCACGTCATATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATGCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTCAGTTCCATTAAGAAAACTCTGGTATTCAc                          >  1:306276/4‑150 (MQ=255)
                                                                                                                                 ccTGATATGACGTGGTCACGCCACATCAGGCAATACTATTCAAAATCGTAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCAGTCCGGCTTTGCTGTTATCGACCTGAGGTCCATTAAGCAATCTCTGGTTTTCACGttt                      <  1:97559/150‑1 (MQ=255)
                                                                                                                                 ccTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGttt                      >  1:172104/1‑150 (MQ=255)
                                                                                                                                    gATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTgg                   <  1:220211/150‑1 (MQ=255)
                                                                                                                                    gATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTgg                   >  2:388005/1‑150 (MQ=255)
                                                                                                                                       atGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGAcg                <  1:349653/150‑1 (MQ=255)
                                                                                                                                        tGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGAcgc               >  1:1870695/1‑150 (MQ=255)
                                                                                                                                         gACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCTTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGAcgcg              <  1:1009628/150‑1 (MQ=255)
                                                                                                                                         gACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGAcgcg              >  2:1378943/1‑150 (MQ=255)
                                                                                                                                         gACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGAcgcg              <  1:149548/150‑1 (MQ=255)
                                                                                                                                         gACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACAACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCATACTCTGGTTTTCACGTTTTGGAcgcg              >  1:1726695/1‑150 (MQ=255)
                                                                                                                                         gACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATACATCAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGAcgcg              <  1:1435105/150‑1 (MQ=255)
                                                                                                                                         gACGAGGTCACGCCACATCAGGCAATACAAATCAAAATCCCAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCCACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGAcgcg              <  1:143441/150‑1 (MQ=255)
                                                                                                                                          aCGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGGTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGACGCGa             >  2:1239496/1‑150 (MQ=255)
                                                                                                                                           cGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGACGCGAc            >  2:1790/1‑150 (MQ=255)
                                                                                                                                           cGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGACGCGAc            <  2:10728/150‑1 (MQ=255)
                                                                                                                                            gtgaccACGTCATATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGACGCGACt           >  1:11930/6‑150 (MQ=255)
                                                                                                                                             tgaccACGTCATATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGACGCGACTc          >  2:1940998/5‑150 (MQ=255)
                                                                                                                                               gTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGCTCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGACGCGcctcct        >  1:31335/1‑150 (MQ=255)
                                                                                                                                                 cACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGACGCGACTCCTGt      >  2:548165/1‑150 (MQ=255)
                                                                                                                                                   cGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGACGCGCCTCCTGTTc    >  1:882470/1‑150 (MQ=255)
                                                                                                                                                   cGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGACGCGACTCCTGTTc    <  2:1529979/150‑1 (MQ=255)
                                                                                                                                                     ccACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGACGCGACTCCTGTTCCg  <  2:598075/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGATGCCAGAGCGCAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGACGCGACTCCTGTTCCG  >  NC_000913/1966225‑1966523

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: