Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,966,368 G→A 8.8% intergenic (‑177/+25) flhB ← / ← cheZ flagellar biosynthesis protein FlhB/chemotaxis protein CheZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,966,3680GA 8.8% 271.1 / 11.6 91intergenic (‑177/+25)flhB/cheZflagellar biosynthesis protein FlhB/chemotaxis protein CheZ
Reads supporting (aligned to +/- strand):  ref base G (45/38);  new base A (4/4);  total (49/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

GTGAATGATGCCAGAGCGCAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGACGCGAC  >  NC_000913/1966220‑1966513
                                                                                                                                                    |                                                                                                                                                 
gTGAATGATGCCAGAGCGCAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGt                                                                                                                                                  <  1:1639925/150‑1 (MQ=255)
     tgatgTCAGAGCGCAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGc                                                                                                                                             >  1:248610/1‑150 (MQ=255)
      gatgCCAGAGCGCAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGcc                                                                                                                                            <  1:1517465/150‑1 (MQ=255)
          ccAGAGCGCAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACAt                                                                                                                                        >  1:192103/1‑150 (MQ=255)
          ccAGAGCGCAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACAt                                                                                                                                        >  1:1581402/1‑150 (MQ=255)
           cAGAGCGCAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATc                                                                                                                                       <  2:1550655/150‑1 (MQ=255)
                cgcAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCTGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCa                                                                                                                                  <  2:1312229/150‑1 (MQ=255)
                 gcAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATGTGGCGTGACCACGTCATATCAGGcgg                                                                                                                                 <  1:53672/150‑3 (MQ=255)
                  cAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAAt                                                                                                                                >  2:466991/1‑150 (MQ=255)
                   aaaGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATa                                                                                                                               <  2:1638621/150‑1 (MQ=255)
                    aaGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATGTGGCGTGACCACGTCATATCAGGcggtct                                                                                                                              <  2:762558/150‑6 (MQ=38)
                     aGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATGTGGCGTGACCACGTCATATCAGGcggtctg                                                                                                                             <  2:63482/150‑7 (MQ=38)
                        gTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAAt                                                                                                                          <  2:248610/150‑1 (MQ=255)
                         ttCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATc                                                                                                                         <  2:953108/150‑1 (MQ=255)
                           cAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCaa                                                                                                                       >  1:912528/1‑150 (MQ=255)
                           cAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCaa                                                                                                                       >  2:1555750/1‑150 (MQ=255)
                            aaTGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCaaa                                                                                                                      >  1:245340/1‑150 (MQ=255)
                             aTGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCaaaa                                                                                                                     <  2:1432222/150‑1 (MQ=255)
                             aTGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCaaaa                                                                                                                     >  2:1447695/1‑150 (MQ=255)
                                gTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATcc                                                                                                                  >  2:954870/1‑150 (MQ=255)
                                gTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATcc                                                                                                                  >  2:509906/1‑150 (MQ=255)
                                 tttGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCa                                                                                                                 <  1:706954/150‑1 (MQ=255)
                                   tGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCTATACAAATCAAAATCCAAg                                                                                                               <  1:1732343/150‑1 (MQ=255)
                                    gAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGa                                                                                                              >  1:601665/1‑150 (MQ=255)
                                     aGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGAc                                                                                                             <  1:490273/150‑1 (MQ=255)
                                         aGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATc                                                                                                         <  2:1335812/150‑1 (MQ=255)
                                          ggggCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATcc                                                                                                        <  2:1565409/150‑1 (MQ=255)
                                          ggggCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATcc                                                                                                        >  2:634315/1‑150 (MQ=255)
                                             gCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCcaac                                                                                                     >  1:750721/1‑150 (MQ=255)
                                              cAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCcaaca                                                                                                    <  2:192103/150‑1 (MQ=255)
                                                  aCAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATc                                                                                                <  2:1899742/150‑1 (MQ=255)
                                                       cGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCa                                                                                           >  1:1721662/1‑150 (MQ=255)
                                                          gATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCt                                                                                        <  1:473323/150‑1 (MQ=255)
                                                              tAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATc                                                                                    >  2:1596663/1‑150 (MQ=255)
                                                               aGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATcc                                                                                   <  1:1549174/150‑1 (MQ=255)
                                                                  gCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGa                                                                                >  2:869941/1‑150 (MQ=255)
                                                                     tttGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTg                                                                             <  2:1205117/150‑1 (MQ=255)
                                                                       tgctgcCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGc                                                                           <  1:17311/150‑1 (MQ=255)
                                                                        gctgcCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGCCGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCTTCCACCTGATCCTGCCTGGCt                                                                          >  2:1049197/1‑150 (MQ=255)
                                                                         ctgcCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTa                                                                         <  2:1705841/150‑1 (MQ=255)
                                                                           gcCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTAcc                                                                       <  1:869941/150‑1 (MQ=255)
                                                                             cacacaTCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACcac                                                                     >  1:1196334/1‑150 (MQ=255)
                                                                                 caTCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCg                                                                 >  2:1945140/1‑150 (MQ=255)
                                                                                    cAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCt                                                              >  1:1498025/1‑150 (MQ=255)
                                                                                      aGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCttt                                                            >  1:798913/1‑150 (MQ=255)
                                                                                      aGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCttt                                                            >  1:713074/1‑150 (MQ=255)
                                                                                        cATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGc                                                          <  2:601665/150‑1 (MQ=255)
                                                                                           aGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTgg                                                       <  2:1581402/150‑1 (MQ=255)
                                                                                             tgtgCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTa                                                     >  1:1586731/1‑150 (MQ=255)
                                                                                              gtgCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTAt                                                    >  2:1346296/1‑150 (MQ=255)
                                                                                               tgCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTTGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATc                                                   <  1:700390/150‑1 (MQ=255)
                                                                                                gCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCg                                                  >  1:1664211/1‑150 (MQ=255)
                                                                                                      gTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGa                                            >  1:974812/1‑150 (MQ=255)
                                                                                                           aTGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTcc                                       <  2:750721/150‑1 (MQ=255)
                                                                                                               ggcgTCATCGCCTTATCAGAGCGCCTGATATGAGGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACATATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACGAGAGGTCCATTa                                   <  2:1065514/150‑1 (MQ=255)
                                                                                                               ggcgTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTa                                   <  2:1346036/150‑1 (MQ=255)
                                                                                                               ggcgTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTa                                   <  1:1447695/150‑1 (MQ=255)
                                                                                                               ggcgTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTa                                   <  2:732000/150‑1 (MQ=255)
                                                                                                                gcgTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTaa                                  <  1:954870/150‑1 (MQ=255)
                                                                                                                gcgTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTaa                                  >  2:43784/1‑150 (MQ=255)
                                                                                                                           ttATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTgg                       <  2:1877886/150‑1 (MQ=255)
                                                                                                                           ttATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTgg                       <  2:245340/150‑1 (MQ=255)
                                                                                                                             aTCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCAGCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAACCTCTGGtt                     >  1:776131/1‑150 (MQ=255)
                                                                                                                             aTCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGtt                     >  1:1933513/1‑150 (MQ=255)
                                                                                                                               cAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTgg                       >  2:737042/1‑146 (MQ=255)
                                                                                                                               cAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTgg                       <  1:737042/146‑1 (MQ=255)
                                                                                                                                aGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTc                  >  2:446195/1‑150 (MQ=255)
                                                                                                                                  aCCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCAc                >  1:805136/1‑150 (MQ=255)
                                                                                                                                   ccGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGCCTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCCCACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACg               >  2:1516160/1‑150 (MQ=255)
                                                                                                                                   ccGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCCCCTGCTCCTGACTGGCTACCCCACCGGCTATGCTGGTATCGCCCTGAGGTCCATTAAGCAAACTCTGGTTTTCACg               >  2:1523329/1‑150 (MQ=255)
                                                                                                                                    tgCCTGATGTGGCGTGACCACGTCATATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAg                                 <  2:395208/130‑1 (MQ=38)
                                                                                                                                    tgCCTGATGTGGCGTGACCACGTCATATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAg                                 >  1:395208/2‑131 (MQ=38)
                                                                                                                                    tgCCTGATGTGGCGTGACCACGTCATATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGt              >  1:140760/2‑150 (MQ=255)
                                                                                                                                    cGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGt              >  1:144920/1‑150 (MQ=255)
                                                                                                                                     gCCTGATGTGGCGTGACCACGTCATATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGtt             >  2:1258545/1‑150 (MQ=255)
                                                                                                                                     gCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGtt             >  2:88992/1‑150 (MQ=255)
                                                                                                                                      ccTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGttt            >  1:1733372/1‑150 (MQ=255)
                                                                                                                                       cTGATGTGGCGTGACCACGTCATATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGtttt           >  2:72098/1‑150 (MQ=255)
                                                                                                                                        tGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTg          <  1:1346296/150‑1 (MQ=255)
                                                                                                                                        tGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTg          >  1:947658/1‑150 (MQ=255)
                                                                                                                                         gATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTgg         >  2:191543/1‑150 (MQ=255)
                                                                                                                                           tatGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGAc       <  2:1721662/150‑1 (MQ=255)
                                                                                                                                           tatGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGAc       <  1:466991/150‑1 (MQ=255)
                                                                                                                                             tGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGAcgc     >  1:383368/1‑150 (MQ=255)
                                                                                                                                             tGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGAcgc     >  1:1484594/1‑150 (MQ=255)
                                                                                                                                             tGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGAcgc     >  1:577819/1‑150 (MQ=255)
                                                                                                                                              gACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGc                                <  1:578679/122‑1 (MQ=255)
                                                                                                                                              gACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGc                                >  2:578679/1‑122 (MQ=255)
                                                                                                                                              gACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGAcgcg    >  1:653207/1‑150 (MQ=255)
                                                                                                                                                cGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGACGCGAc  <  2:457129/150‑1 (MQ=255)
                                                                                                                                                cGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTCATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCAATAAGCAAACTCTGGTTTTCACGTTTTGGACGCGAc  >  2:1179343/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GTGAATGATGCCAGAGCGCAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGACGCGAC  >  NC_000913/1966220‑1966513

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