Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,233,497 G→A 5.1% intergenic (+130/+100) PP_5626 → / ← PP_3711 repressor‑like DNA‑binding domain‑containing protein/diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,233,4970GA 5.1% 333.5 / 4.1 111intergenic (+130/+100)PP_5626/PP_3711repressor‑like DNA‑binding domain‑containing protein/diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (54/50);  new base A (2/4);  total (56/55)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 114 total aligned reads displayed.

ATCGATCGGTGGAAGGCTGAGGGAGGAAGGCACAAGCTGGCTTACTGAAAGCTCAGGGCTGCTCGGTTTGGCTGTAGACCTGCTCAGGTTCACTCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGCCACGGGTCGGGTGTTCC  >  NC_002947/4233348‑4233645
                                                                                                                                                     |                                                                                                                                                    
atcgatcgGTGGAAGGCTGAGGGAGGAAGGCACAAGCTGGCTTACTGAAAGCTCAGGGCTGCTCGGTTTGGCTGTAGACCTGCTCAGGTTCACTCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTgg                                                                                                                                                      <  2:1047465/150‑1 (MQ=255)
  cgatcgGTGGAAGGCTGAGGGAGGAAGGCACAAGCTGGCTTACTGAAAGCTCAGGGCTGCTCGGTTTGGCTGTAGACCTGCTCAGGTTCACTCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTgggg                                                                                                                                                    <  2:1784721/150‑1 (MQ=255)
     tcgGTGGAAGGCTGAGGGAGGAAGGCACAAGCTGGCTTACTGAAAGCTCAGGGCTGCTCGGTTTGGCTGTAGACCTGCTCAGGTTCACTCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCtt                                                                                                                                                 >  1:1350067/1‑150 (MQ=255)
      cgGTGGAAGGCTGAGGGAGGAAGGCACAAGCTGGCTTACTGAAAGCTCAGGGCTGCTCGGTTTGGCTGTAGACCTGCTCAGGTTCACTCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCttt                                                                                                                                                <  2:1151578/150‑1 (MQ=255)
        gTGGAAGGCTGAGGGAGGAAGGCACAAGCTGGCTTACTGAAAGCTCAGGGCTGCTCGGTTTGGCTGTAGACCTGCTCAGGTTCACTCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGa                                                                                                                                              <  2:161897/150‑1 (MQ=255)
         tGGAAGGCTGAGGGAGGAAGGCACAAGCTGGCTTACTGAAAGCTCAGCGCTGCTCGGTTTGGCTGTAGACCTGCTCAGGTTCACTCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGAc                                                                                                                                             <  2:523313/150‑1 (MQ=255)
          ggAAGGCTGAGGGAGGAAGGCACAAGCTGGCTTACTGAAAGCTCAGGGCTGCTCGGTTTGGCTGTAGACCTGCTCAGGTTCACTCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGAcc                                                                                                                                            <  2:768171/150‑1 (MQ=255)
              gcatggtcatgccgatggacttcattttctagacacataccccttcagtgtcaaaacccacaaggttttgactgctcctgctcctgctcctgctcctgctcctgctcctgctcctgctttaggTCAAAGCCTTTGGGGCTTTGACCTaaa                                                                                                                                        <  2:1987119‑M1/27‑1 (MQ=255)
                 tGAGGGAGGAAGGCACAAGCTGGCTTACTGAAAGCTCAGGGCTGCTCGGTTTGGCTGTAGACCTGCTCAGGTTCACTCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAagca                                                                                                                                     >  2:2163761/1‑150 (MQ=255)
                 tGAGGGAGGAAGGCACAAGCTGGCTTACTGAAAGCTCAGGGCTGCTCGGTTTGGCTGTAGACCTGCTCAGGTTCACTCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTCGACCTAAagca                                                                                                                                     >  2:1221715/1‑150 (MQ=255)
                    gggaggAAGGCACAAGCTGGCTTACTGAAAGCTCAGGGCTGCTCGGTTTGGCTGTAGACCTGCTCAGGTTCACTCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAagcagga                                                                                                                                  <  1:1923296/150‑1 (MQ=255)
                            ggCAAAAGCTGGGGTACTTAAAGCTCCGGGCTGCTCGGTTTGGCAGTAGACCTCCTCAGCTTCACGCACAATAAATCAGCACCGAGCTTTGCTCTTTCAATCAAGTGAATCAAAGCTTTTTGGTCGTTGACTTAAagccggagtaggagc                                                                                                                          >  2:193555/1‑150 (MQ=255)
                                   gCTGGCTTACTGAAAGCTCAGGGCTGCTCGGTTTGGCTGTAGACCTGCTCAGGTTCACTCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTCGACCTAAagcaggagcaggagcaggagca                                                                                                                   <  1:1221715/150‑1 (MQ=255)
                                       gCTTACTGAAAGCTCAGGGCTGCTCGGTTTGGCTGTAGACCTGCTCAGGTTCACTCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAagcag                                                                                                                                    >  1:361931/1‑129 (MQ=255)
                                        cTTACTGAAAGCTCAGGGCTGCTCGGTTTGGCTGTAGACCTGCTCAGGTTCACTCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTaaa                                                                                                                                        <  2:1866048/124‑1 (MQ=255)
                                        cTTACTGAAAGCTCAGGGCTGCTCGGTTTGGCTGTAGACCTGCTCAGGTTCACTCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTaaa                                                                                                                                        >  1:1866048/1‑124 (MQ=255)
                                         ttACTGAAAGCTCAGGGCTGCTCGGTTTGGCTGTAGACCTGCTCAGGTTCACTCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCttt                                                                                                                                                >  2:2140368/1‑115 (MQ=255)
                                           aCTGAAAGCTCAGGGCTGCTCGGTTTGGCTGTAGACCTGCTCAGGTTCACTCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAagcaggagcaggagcaggagcaggagcagg                                                                                                           >  2:1488466/1‑150 (MQ=255)
                                            cTGAAAGCTCAGGGCTGCTCGGTTTGGCTGTAGACCTGCTCAGGTTCACTCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAagcaggagcaggagcaggagcaggagcagga                                                                                                          >  1:502737/1‑150 (MQ=255)
                                               aaaGCTCAGGGCTGCTCGGTTTGGCTGTAGACCTGCTCAGGTTCACTCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGAcc                                                                                                                                            >  2:500360/1‑113 (MQ=255)
                                               aaaGCTCAGGGCTGCTCGGTTTGGCTGTAGACCTGCTCAGGTTCACTCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGAcc                                                                                                                                            <  1:500360/113‑1 (MQ=255)
                                                aaGCTCAGGGCTGCTCGGTTTGGCTGTAGACCTGCTCAGGTTCACTCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAagcaggagcaggagcag                                                                                                                        <  1:1509173/132‑1 (MQ=255)
                                                aaGCTCAGGGCTGCTCGGTTTGGCTGTAGACCTGCTCAGGTTCACTCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAagcaggagcaggagcag                                                                                                                        >  2:1509173/1‑132 (MQ=255)
                                                  gCTCAGGGCTGCTCGGTTTGGCTGTAGACCTGCTCAGGTTCACTCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAagcaggagcaggagcaggagcaggagcaggagcagga                                                                                                    >  2:1770855/1‑150 (MQ=255)
                                                   cTCAGGGCTGCTCGGTTTGGCTGTAGACCTGCTCAGGTTCACTCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAagcaggagcaggagcaggagcaggagcaggagcaggag                                                                                                   <  2:1444998/150‑1 (MQ=255)
                                                      aGGGCTGCTCGGTTTGGCTGTAGACCTGCTCAGGTTCACTCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAagcaggagcaggagcaggagcaggagcaggagcaggagcag                                                                                                <  1:549367/150‑1 (MQ=255)
                                                      aGGGCTGCTCGGTTTGGCTGTAGACCTGCTCAGGTTCACTCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAagcaggagcaggagcaggagcaggagcaggagcaggagcag                                                                                                >  2:549367/1‑150 (MQ=255)
                                                           tgctCGGTTTGGCTGTAGACCTGCTCAGGTTCACTCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAACGCCTTTGGGGCTTTGACCTAAAGCAggagcagaagcaggagcaggagcaggagcaggagcaggagca                                                                                           <  2:807225/150‑1 (MQ=255)
                                                           tgctCGGTTTGGCTGTAGACCTGCTCAGGTTCACTCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAagcag                                                                                                                                    >  1:2048836/1‑109 (MQ=255)
                                                           tgctCGGTTTGGCTGTAGACCTGCTCAGGTTCACTCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAagcag                                                                                                                                    <  2:2048836/109‑1 (MQ=255)
                                                             ctCGGTTTGACTGTAGACCTGCTCACGTTCACTCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGt                                                                                         <  1:2226583/150‑1 (MQ=255)
                                                                        tGTAGACCTGCTCAGGTTCACTCCTCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTCAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTtgt                                                                              >  1:1639471/1‑150 (MQ=255)
                                                                        tGTAGACCTGCTCAGGTTCACTCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAagcaggagcaggagcag                                                                                                                        <  2:622045/108‑1 (MQ=255)
                                                                              ccTGCTCAGGTTCACTCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGttt                                                                        <  2:2283781/150‑1 (MQ=255)
                                                                               cTGCTCAGGTTCACTCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAagcaggagcaggagcaggagcaggagcaggag                                                                                                         >  2:118907/1‑116 (MQ=255)
                                                                               cTGCTCAGGTTCACTCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAagcaggagcaggagcaggagcaggagcaggag                                                                                                         <  1:118907/116‑1 (MQ=255)
                                                                               cTGCTCAGGTTCACTCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGATTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGtttt                                                                       >  1:2131911/1‑150 (MQ=255)
                                                                                      ggTTCACTCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTg                                                                >  1:5682/1‑150 (MQ=255)
                                                                                       gTTCACTCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAagcaggagcaggagcaggagcaggagcag                                                                                                            <  2:386616/105‑1 (MQ=255)
                                                                                       gTTCACTCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAagcaggagcaggagcaggagcaggagcag                                                                                                            >  1:386616/1‑105 (MQ=255)
                                                                                       gTTCACTCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGa                                                               >  2:1838544/1‑150 (MQ=255)
                                                                                        ttCACTCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGaa                                                              >  1:1138733/1‑150 (MQ=255)
                                                                                             tCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAagcaggagcaggagcaggagcaggagcaggagcaggagcag                                                                                                <  1:2122949/111‑1 (MQ=255)
                                                                                             tCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTg                                                                      <  1:2079050/137‑1 (MQ=255)
                                                                                             tCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTg                                                                      >  2:2079050/1‑137 (MQ=255)
                                                                                             tCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGt                                                         >  2:879567/1‑150 (MQ=255)
                                                                                             tCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGt                                                         >  1:1329969/1‑150 (MQ=255)
                                                                                               cccATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGAc                                                                                                                                             <  2:928258/64‑1 (MQ=255)
                                                                                               cccATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGAc                                                                                                                                             >  1:928258/1‑64 (MQ=255)
                                                                                               cccATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAagcaggagcag                                                                                                                              >  2:2098299/1‑79 (MQ=255)
                                                                                               cccATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAagcaggagcag                                                                                                                              <  1:2098299/79‑1 (MQ=255)
                                                                                               cccATCAAACAGCACCGAGCTTTGCTATTTCAATCAAGAGCATCAAAGCCTATGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTAt                                                       >  2:370439/1‑150 (MQ=255)
                                                                                                ccATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGAAGGAGCAGGAAAAATCAAATGCTTGTTGAATTTGACAGGGAAGGGGTAtg                                                      >  2:2059419/1‑150 (MQ=255)
                                                                                                  aTCGAACAGCAGCGGGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTGGGAGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTAtgtg                                                    >  2:1259026/1‑150 (MQ=255)
                                                                                                  aTCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTAtgtg                                                    >  1:1216644/1‑150 (MQ=255)
                                                                                                   tCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATgtgt                                                   <  1:1838544/150‑1 (MQ=255)
                                                                                                     aaaCAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCt                                                 >  1:1209476/1‑150 (MQ=255)
                                                                                                         aGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGaa                                             >  2:2114236/1‑150 (MQ=255)
                                                                                                             ccGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAagcaggagcaggagcaggagcagga                                                                                                                >  2:1187425/1‑79 (MQ=255)
                                                                                                             ccGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAagcaggagcaggagcaggagcagga                                                                                                                <  1:1187425/79‑1 (MQ=255)
                                                                                                              cGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGa                                        >  1:1847011/1‑150 (MQ=255)
                                                                                                               gAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATTTGTCTAGAAAATGaa                                       <  1:879567/150‑1 (MQ=255)
                                                                                                               gAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGaa                                       >  1:1316572/1‑150 (MQ=255)
                                                                                                                aGCTTTACTCTTTCAATCAAGAGCATCAAAACCATTGAGTCTTGGACCTTAAGCAAGATCATGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAg                                      <  1:2059419/150‑1 (MQ=255)
                                                                                                                 gCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGt                                     <  1:2163761/150‑1 (MQ=255)
                                                                                                                   tttGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTtg                                                                               <  1:2329313/106‑1 (MQ=255)
                                                                                                                   tttGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTtg                                                                               >  2:2329313/1‑106 (MQ=255)
                                                                                                                      gCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGAAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGACTAGAAAATGAAGTCCATc                                >  2:517619/1‑150 (MQ=255)
                                                                                                                      gCTCTTTCAATCAAGAGAAACAAAGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGTAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATc                                <  2:1639471/150‑1 (MQ=255)
                                                                                                                       ctctTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCg                               >  2:145044/1‑150 (MQ=255)
                                                                                                                       ctctTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCg                               <  2:1216644/150‑1 (MQ=255)
                                                                                                                        tctTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCgg                              >  2:401980/1‑150 (MQ=255)
                                                                                                                          tttCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAggg                                                           >  2:1185481/1‑119 (MQ=255)
                                                                                                                          tttCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAggg                                                           <  1:1185481/119‑1 (MQ=255)
                                                                                                                             tcctgctttaggTCAAAGCCCCAAAGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTgg                                                                            <  2:1753467‑M1/87‑1 (MQ=255)
                                                                                                                             tcctgctttaggTCAAAGCCCCAAAGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTgg                                                                            >  1:1753467‑M1/13‑99 (MQ=255)
                                                                                                                             tcctgctttaggTCAAAGCCCCAAAGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGa                         <  1:1987119‑M1/138‑1 (MQ=255)
                                                                                                                                tgctttaggTCAAAGCCCCAAAGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTtgt                                                                              <  1:915636‑M1/85‑1 (MQ=255)
                                                                                                                                tgctttaggTCAAAGCCCCAAAGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTtgt                                                                              >  2:915636‑M1/10‑94 (MQ=255)
                                                                                                                                tcaaGAGCATCAAAGCCTTTGGGGCTTTGACCTAAagcaggagcaggagcaggagcaggagcaggag                                                                                                         >  1:962875/1‑67 (MQ=255)
                                                                                                                                tcaaGAGCATCAAAGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTAtgt                                                     >  1:2327594/1‑119 (MQ=255)
                                                                                                                                tcaaGAGCATCAAAGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTAtgt                                                     <  2:2327594/119‑1 (MQ=255)
                                                                                                                                tcaaGAGCATCAAAGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCa                      >  1:368449/1‑150 (MQ=255)
                                                                                                                                 caaGAGCATCAAAGCCTTTGGGGCTTTGACCTAAagcaggagcaggagcaggagcaggagcag                                                                                                            >  2:1667898/1‑63 (MQ=255)
                                                                                                                                 caaGAGCATCAAAGCCTTTGGGGCTTTGACCTAAagcaggagcaggagcaggagcaggagcag                                                                                                            <  1:1667898/63‑1 (MQ=255)
                                                                                                                                      aggTCAAAGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTa                                                        <  2:837753‑M1/107‑1 (MQ=255)
                                                                                                                                      aggTCAAAGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTa                                                        >  1:837753‑M1/4‑110 (MQ=255)
                                                                                                                                       cATCAAAGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTCGGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGCCAcg               >  1:1776763/1‑150 (MQ=255)
                                                                                                                                         TCAAAGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGCCAcggg             <  1:2277963‑M1/150‑1 (MQ=255)
                                                                                                                                             AGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGCCAcgggtcgg         <  1:1030600‑M1/150‑1 (MQ=255)
                                                                                                                                               CCTTTGGGGCTTTGACCTAAAGCAGGAGCAGTAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGCCAcgggtcgggt       >  2:1461874‑M1/1‑150 (MQ=255)
                                                                                                                                               CCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGaa                                       <  1:1420398‑M1/118‑1 (MQ=255)
                                                                                                                                               CCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGaa                                       >  2:1420398‑M1/1‑118 (MQ=255)
                                                                                                                                               CCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGCCAcgggtcgggt       <  1:958779‑M1/150‑1 (MQ=255)
                                                                                                                                               CCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGCCAcgggtcgggt       >  2:1209476‑M1/1‑150 (MQ=255)
                                                                                                                                                CTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGt                                                                                         <  1:2328760‑M1/67‑1 (MQ=255)
                                                                                                                                                CTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGt                                                                                         >  2:2328760‑M1/1‑67 (MQ=255)
                                                                                                                                                 tttGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGCCAcg               <  2:2236906/140‑1 (MQ=255)
                                                                                                                                                 tttGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGCCAcg               >  1:2236906/1‑140 (MQ=255)
                                                                                                                                                    ggggCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGCCACGGGTCGGGTGTTcc  <  1:1488466/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATCGATCGGTGGAAGGCTGAGGGAGGAAGGCACAAGCTGGCTTACTGAAAGCTCAGGGCTGCTCGGTTTGGCTGTAGACCTGCTCAGGTTCACTCCCCATCAAACAGCACCGAGCTTTGCTCTTTCAATCAAGAGCATCAAAGCCTTTGGGGCTTTGACCTAAAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGGAGCAGTCAAAACCTTGTGGGTTTTGACACTGAAGGGGTATGTGTCTAGAAAATGAAGTCCATCGGCATGACCATGCCACGGGTCGGGTGTTCC  >  NC_002947/4233348‑4233645

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