Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,966,377 C→T 5.1% intergenic (‑186/+16) 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,3770CT 5.1% 372.7 / 5.8 116intergenic (‑186/+16)flhB/cheZflagellin export apparatus, substrate specificity protein/chemotaxis regulator, protein phosphatase for CheY
Reads supporting (aligned to +/- strand):  ref base C (59/51);  new base T (3/3);  total (62/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.30e-01

Only 100 of 124 total aligned reads displayed.

GCCAGAGCGCAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGACGCGACTCCTGTTCCGGG  >  NC_000913/1966229‑1966525
                                                                                                                                                    |                                                                                                                                                    
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)
       cacAAAGCGTTCAATGGTGTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATTGTGTGCGTTTGGCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATGTGGCGTGACCACGTCATATCAGGcg                                                                                                                                              <  2:1259556/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)
           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                                                                                                                           <  2:1003464/150‑1 (MQ=255)
                           gAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAgc                                                                                                                          >  1:497257/1‑149 (MQ=255)
                           gAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGa                                                                                                                          >  1:19733/1‑150 (MQ=255)
                               aaGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTAt                                                                                                                      >  1:606990/1‑150 (MQ=255)
                                  gggCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCca                                                                                                                   >  1:1529980/1‑150 (MQ=255)
                                   ggCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCcaa                                                                                                                  >  1:742304/1‑150 (MQ=255)
                                     cAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCTTCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCGACATCAGTCAATACAAAACAAAATCCACGACTATCaaaca                                                                                                                >  2:1817417/1‑150 (MQ=255)
                                     cAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCcaaca                                                                                                                <  1:767017/150‑1 (MQ=255)
                                         aCAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATc                                                                                                            <  1:1389315/150‑1 (MQ=255)
                                          cAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCg                                                                                                           <  2:1870753/150‑1 (MQ=255)
                                          cAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCg                                                                                                           >  2:1635761/1‑150 (MQ=255)
                                           aGGCGGGAGTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGt                                                                                                          <  2:258942/150‑1 (MQ=255)
                                            ggCGGGATATAGGGCTTTTGCTGCCAGACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCGTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTc                                                                                                         <  2:100402/150‑1 (MQ=255)
                                             gCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTcc                                                                                                        >  2:26084/1‑150 (MQ=255)
                                             gCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTcc                                                                                                        <  1:1931495/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)
                                                        ggCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACACATCGTCCACCTGATCCTg                                                                                             >  1:552213/1‑150 (MQ=255)
                                                                 tgcCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTAc                                                                                    >  2:713247/1‑150 (MQ=255)
                                                                 tgcCACACATCAAGCATAGTGTGCGTTCGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTAc                                                                                    >  1:342954/1‑150 (MQ=255)
                                                                  gcCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTAcc                                                                                   >  1:2017183/1‑150 (MQ=255)
                                                                    cacacaTCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACcac                                                                                 <  2:157099/150‑1 (MQ=255)
                                                                        caTCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCg                                                                             >  2:415145/1‑150 (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)
                                                                                   gggtgCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGAGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGt                                                                  <  1:1244200/148‑1 (MQ=255)
                                                                                    tgtgCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTa                                                                 >  1:648637/1‑150 (MQ=255)
                                                                                      tgCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACAACACCGGCTTTGCTGGTAta                                                               >  2:1301953/1‑149 (MQ=255)
                                                                                       gCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCg                                                              >  1:1607122/1‑150 (MQ=255)
                                                                                       gCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCg                                                              >  1:10728/1‑150 (MQ=255)
                                                                                       gCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCg                                                              <  1:1073826/150‑1 (MQ=255)
                                                                                         gTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATACTGACTGGCTACCACACCGGCTTTGCGGGTATCGAc                                                            >  1:1055971/1‑150 (MQ=255)
                                                                                          tttGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCCCACCGGCTTTGCTGGTATCGAcc                                                           >  2:141825/1‑150 (MQ=255)
                                                                                          tttGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGAcc                                                           <  1:1581017/150‑1 (MQ=255)
                                                                                             gTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACTCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGa                                                        <  2:310247/150‑1 (MQ=255)
                                                                                                ggATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGGCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGt                                                     <  2:1455941/150‑1 (MQ=255)
                                                                                                    gcggcgTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCAt                                                 <  1:1913114/150‑1 (MQ=255)
                                                                                                     cggcgTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGCCCTGAGGTCCAtt                                                >  2:2039379/1‑150 (MQ=255)
                                                                                                     cggcgTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCAtt                                                >  2:217934/1‑150 (MQ=255)
                                                                                                     cggcgTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCAACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCAtt                                                >  2:301705/1‑150 (MQ=255)
                                                                                                     cggcgTCATCGCCGTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAACTCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCAtt                                                <  2:588786/150‑1 (MQ=255)
                                                                                                      ggcgTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCAAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTa                                               >  1:593890/1‑150 (MQ=255)
                                                                                                        cgTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAg                                             <  1:998725/150‑1 (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)
                                                                                                                               tGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTg                      <  2:192518/150‑1 (MQ=255)
                                                                                                                               tGATATGACGTGGTCACGCCACATAAGGCAATACAAATCAAAATCCAAGACTATGCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCGGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTg                      <  2:141634/150‑1 (MQ=255)
                                                                                                                                gATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAGTCCAAGACTATCCAACAAAGCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGCCCTGAGGTCCATTCAGCAAACTCTGGTTTTGACGTTTTgg                     >  2:152562/1‑150 (MQ=255)
                                                                                                                                gATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTgg                     >  2:388005/1‑150 (MQ=255)
                                                                                                                                gATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTgg                     <  1:220211/150‑1 (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)
                                                                                                                                             cACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGACGCGACTCCTGt        >  2:548165/1‑150 (MQ=255)
                                                                                                                                               cGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGACGCGCCTCCTGTTc      >  1:882470/1‑150 (MQ=255)
                                                                                                                                                 ccACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGACGCGACTCCTGTTCCg    <  2:598075/150‑1 (MQ=255)
                                                                                                                                                  cacaTCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTCCCACCCCGGCTTTGCTGGTATCGCCCTGAGGTCCAGTAAGCACACTCTGGTTTTCCCGTTTTGGACGCGACTCCTGTTCCgg   >  2:517280/1‑150 (MQ=255)
                                                                                                                                                  cacaTCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGACGCGACTCCTGTTCCgg   >  2:1352921/1‑150 (MQ=255)
                                                                                                                                                  aacaTCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGACGCGACTCCTGTTCCgg   <  2:1785526/149‑1 (MQ=255)
                                                                                                                                                   acaTCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGACGCGACTCCTGTTCCggg  <  2:679988/150‑1 (MQ=255)
                                                                                                                                                   acaTCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGACGCGACTCCTGTTCCggg  <  2:1981371/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCCAGAGCGCAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGACGCGACTCCTGTTCCGGG  >  NC_000913/1966229‑1966525

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