Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,966,360 A→G 6.5% intergenic (‑169/+33) 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,3600AG 6.5% 229.4 / 5.1 77intergenic (‑169/+33)flhB/cheZflagellar biosynthesis protein FlhB/chemotaxis protein CheZ
Reads supporting (aligned to +/- strand):  ref base A (40/32);  new base G (3/2);  total (43/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-01

ATTAAGCGTGAATGATGCCAGAGCGCAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGAC  >  NC_000913/1966213‑1966508
                                                                                                                                                   |                                                                                                                                                    
aTTAAGCGTGAATGATGCCAGAGCGCAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATg                                                                                                                                                    <  2:828696/150‑1 (MQ=255)
 ttgaGCGTGAATGATGCCAGAGCGCAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGa                                                                                                                                                   >  1:1506379/4‑150 (MQ=255)
      cGTGAATGATGCCAGAGCGCAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTgg                                                                                                                                              <  1:1155175/150‑1 (MQ=255)
      cGTGAATGATGCCAGAGCGCAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTgg                                                                                                                                              <  1:962050/150‑1 (MQ=255)
                gCCAGAGCGCAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATGTGGCGTGACCACGTcata                                                                                                                                    >  1:613360/1‑148 (MQ=255)
                  cAGAGCGCAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATc                                                                                                                                  >  1:1353194/1‑150 (MQ=255)
                       cgcAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCa                                                                                                                             <  1:1105519/150‑1 (MQ=255)
                         cAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATGTGGCGTGACCACGTCATATCAGGcggt                                                                                                                           <  1:1348256/150‑4 (MQ=255)
                         cAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAAt                                                                                                                           <  1:610480/150‑1 (MQ=255)
                         cAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAAt                                                                                                                           >  2:527413/1‑150 (MQ=255)
                           aaGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATAc                                                                                                                         <  2:435760/150‑1 (MQ=255)
                            aGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACa                                                                                                                        <  2:488490/150‑1 (MQ=255)
                             gCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATGTGGCGTGACCACGTCATATCAGGcggtctga                                                                                                                       <  2:613360/150‑8 (MQ=38)
                               gTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAAt                                                                                                                     >  1:715482/1‑150 (MQ=255)
                               gTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCCGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAAt                                                                                                                     >  2:1383450/1‑150 (MQ=255)
                                ttCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATc                                                                                                                    <  1:527413/150‑1 (MQ=255)
                                ttCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATc                                                                                                                    >  1:109437/1‑150 (MQ=255)
                                 tCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCa                                                                                                                   >  1:835226/1‑150 (MQ=255)
                                  cAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCaa                                                                                                                  >  2:807888/1‑150 (MQ=255)
                                    aTGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCGACATCAGGCAATACAAATCaaaa                                                                                                                <  2:1080338/150‑1 (MQ=255)
                                        tttGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCa                                                                                                            <  1:1332719/150‑1 (MQ=255)
                                        tttGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCa                                                                                                            <  2:349589/150‑1 (MQ=255)
                                          tGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAg                                                                                                          >  1:1263135/1‑150 (MQ=255)
                                          tGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAg                                                                                                          >  2:1007397/1‑150 (MQ=255)
                                            aGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGAc                                                                                                        <  1:800130/150‑1 (MQ=255)
                                             gTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACt                                                                                                       >  1:226022/1‑150 (MQ=255)
                                             gTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACt                                                                                                       >  1:369019/1‑150 (MQ=255)
                                              tAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTa                                                                                                      <  1:606427/150‑1 (MQ=255)
                                                 ggggCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATcc                                                                                                   <  2:356079/150‑1 (MQ=255)
                                                 ggggCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATcc                                                                                                   <  1:1598499/150‑1 (MQ=255)
                                                 ggggAAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATcc                                                                                                   >  2:814798/1‑150 (MQ=255)
                                                  gggCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCca                                                                                                  >  1:475625/1‑150 (MQ=255)
                                                  gggCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCca                                                                                                  <  1:841014/150‑1 (MQ=255)
                                                    gCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCcaac                                                                                                <  2:1004313/150‑1 (MQ=255)
                                                            ggCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTc                                                                                        >  2:405770/1‑150 (MQ=255)
                                                               gggATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCAc                                                                                     >  2:739653/1‑150 (MQ=255)
                                                                ggATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCAcc                                                                                    <  1:539079/150‑1 (MQ=255)
                                                                       gggCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCt                                                                             >  2:445810/1‑150 (MQ=255)
                                                                          cTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGAc                                                                          <  2:1604603/150‑1 (MQ=255)
                                                                                ctgcCACACCTCAAGCATAGTGTGCGTTTGTAGGATGCGGCGTCATAGCATTATCATACCGCCTGATATTACTTGGTAACGCCACATAATGCAAGACAAATCAAAATCCAAGAATAGCCAACAAATATTCACCCTGATCCTGAATGGCTa                                                                    >  2:966460/1‑150 (MQ=255)
                                                                                 tgcCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTAc                                                                   >  1:896164/1‑150 (MQ=255)
                                                                                  gcCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTAcc                                                                  >  1:1457633/1‑150 (MQ=255)
                                                                                      cacaTCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCacac                                                              >  2:816962/1‑150 (MQ=255)
                                                                                         aTCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCgg                                                           >  2:1257290/1‑150 (MQ=255)
                                                                                           cAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCt                                                         >  2:1285991/1‑150 (MQ=255)
                                                                                            aaGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCtt                                                        >  1:959684/1‑150 (MQ=255)
                                                                                            aaGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCtt                                                        <  2:835226/150‑1 (MQ=255)
                                                                                                 tAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTg                                                   <  2:1353194/150‑1 (MQ=255)
                                                                                                 tAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTg                                                   <  1:1350776/150‑1 (MQ=255)
                                                                                                  aGTGTGCGTTTGTAGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCCGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCCCCTGATCCTGAATGGCTACCACACCGGCTTTGCTgg                                                  >  2:354405/1‑150 (MQ=255)
                                                                                                   gtgtGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGt                                                 >  2:1262636/1‑150 (MQ=255)
                                                                                                    tgtgCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTa                                                >  2:1327742/1‑150 (MQ=255)
                                                                                                          tttGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGAcc                                          <  2:959684/150‑1 (MQ=255)
                                                                                                             gTCGGATGCGGCGTCATCGCCTTATCCGACCGCCTGATATGACGTGGTCACGCCACAGCAGGCAAGACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGCCCTgg                                       >  2:1366747/1‑149 (MQ=255)
                                                                                                                 gATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTc                                   >  1:258769/1‑150 (MQ=255)
                                                                                                                      ggcgTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTa                              >  1:65621/1‑150 (MQ=255)
                                                                                                                      ggcgTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTa                              <  2:109437/150‑1 (MQ=255)
                                                                                                                         gTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGc                           >  2:760431/1‑150 (MQ=255)
                                                                                                                          tCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCa                          <  2:715482/150‑1 (MQ=255)
                                                                                                                                ccTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAActct                    >  2:1070625/1‑150 (MQ=255)
                                                                                                                                   tATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGt                 <  1:739653/150‑1 (MQ=255)
                                                                                                                                   tATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTAACCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGt                 >  2:236983/1‑150 (MQ=255)
                                                                                                                                    aTCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGtt                >  1:606237/1‑150 (MQ=255)
                                                                                                                                     tCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGttt               <  1:1285991/150‑1 (MQ=255)
                                                                                                                                      cAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGtttt              >  2:938617/1‑150 (MQ=255)
                                                                                                                                         aCCGCCTGATATGACGTGGTCACGCCACATCAGGCAATAGAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTTGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCAc           <  2:413522/150‑1 (MQ=255)
                                                                                                                                         aCCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCAc           <  1:445810/150‑1 (MQ=255)
                                                                                                                                         aCCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCAc           >  1:533769/1‑150 (MQ=255)
                                                                                                                                          ccGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACg          <  2:258769/150‑1 (MQ=255)
                                                                                                                                           cGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGt         >  1:750582/1‑150 (MQ=255)
                                                                                                                                           cGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGt         >  2:1212144/1‑150 (MQ=255)
                                                                                                                                             ccTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGttt       <  2:475625/150‑1 (MQ=255)
                                                                                                                                              cTGATGTGGCGTGACCACGTCATATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGCCCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGtttt      >  2:487167/1‑150 (MQ=255)
                                                                                                                                              cTGATGTGGCGTGACCACGTCATATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGtttt      >  1:1091448/1‑150 (MQ=255)
                                                                                                                                                gATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTg                                                                        <  1:1120934/82‑1 (MQ=255)
                                                                                                                                                gATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTg                                                                        >  2:1120934/1‑82 (MQ=255)
                                                                                                                                                gATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGCCCTGAGGTCCATTAAGCAAACTCTGGGGTTCACGTTTTgg    >  2:915237/1‑150 (MQ=255)
                                                                                                                                                  tatGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGAc  >  2:731235/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ATTAAGCGTGAATGATGCCAGAGCGCAAAGCGTTCAATGGTTTGAGTAAGGGGCAAAACAGGCGGGATTTAGGGCTTTTGCTGCCACACATCAAGCATAGTGTGCGTTTGTCGGATGCGGCGTCATCGCCTTATCAGACCGCCTGATATGACGTGGTCACGCCACATCAGGCAATACAAATCAAAATCCAAGACTATCCAACAAATCGTCCACCTGATCCTGACTGGCTACCACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGAC  >  NC_000913/1966213‑1966508

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