Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,517,521 T→C 8.7% L104L (TTA→TTG pdeA ← putative c‑di‑GMP phosphodiesterase PdeA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,517,5210TC 8.7% 433.5 / 23.1 151L104L (TTA→TTGpdeAputative c‑di‑GMP phosphodiesterase PdeA
Reads supporting (aligned to +/- strand):  ref base T (76/60);  new base C (9/4);  total (87/64)
Fisher's exact test for biased strand distribution p-value = 3.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

Only 100 of 154 total aligned reads displayed.

GATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCGCCAGCCAAAGATAAACAAGAT  >  NC_000913/2517372‑2517670
                                                                                                                                                     |                                                                                                                                                     
gATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGt                                                                                                                                                       <  1:450943/150‑1 (MQ=255)
gATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGt                                                                                                                                                       <  1:161185/150‑1 (MQ=255)
 atatCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTa                                                                                                                                                      >  2:1288729/1‑150 (MQ=255)
 atatCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTa                                                                                                                                                      >  1:523350/1‑150 (MQ=255)
     cTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAatat                                                                                                                                                  <  2:97311/150‑1 (MQ=255)
       ttAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTg                                                                                                                                                <  1:1835910/150‑1 (MQ=255)
           cGGGAAATCAAAGAAGCTCCCCACAAGATACATGCGGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGCCGTTGCGAGACGTATCGGCTGGGCAATCCATAACGCCCACGCGTCCCCAGCTGCCAGCGGAATATTGCGCa                                                                                                                                            >  2:1740965/1‑150 (MQ=255)
            gggAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAg                                                                                                                                           >  1:1009590/1‑150 (MQ=255)
             ggAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAgg                                                                                                                                          <  1:2439841/150‑1 (MQ=255)
             ggAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAgg                                                                                                                                          >  2:884484/1‑150 (MQ=255)
              gAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGc                                                                                                                                         <  2:800238/150‑1 (MQ=255)
              gAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGc                                                                                                                                         <  1:724566/150‑1 (MQ=255)
               aaaTCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGca                                                                                                                                        <  2:2256750/150‑1 (MQ=255)
               aaaTCAAACAAGCACCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGca                                                                                                                                        >  2:1300707/1‑150 (MQ=255)
                   cAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGcacaca                                                                                                                                    >  1:1675533/1‑150 (MQ=255)
                   cAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGcacaca                                                                                                                                    >  1:2027731/1‑150 (MQ=255)
                    aaaGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACacaa                                                                                                                                   >  1:1198590/1‑150 (MQ=255)
                     aagaagCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCTGCGCAATATTACGCTGGcagctggg                                                                                                                                  <  2:612780/150‑7 (MQ=255)
                      agaagCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACacaaca                                                                                                                                 >  2:214634/1‑150 (MQ=255)
                       gaagCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCCGCGGAATATTGCGCAGGcacaccacac                                                                                                                                >  2:900180/1‑150 (MQ=255)
                        aagCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACg                                                                                                                               <  1:375723/150‑1 (MQ=255)
                         agCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCCAACCAAGCCAGAACAGACGATGCCAGACATATCGGCTGGTCAAGCCATAACGCCAACGCGTCACCAGCTGACAGCGTAATATTGCGCATGCACACAACAAGa                                                                                                                              >  2:1087063/1‑150 (MQ=255)
                               ccACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACg                                                                                                                        <  2:1583887/150‑1 (MQ=255)
                                 acaAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATTGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTc                                                                                                                      <  1:2070420/150‑1 (MQ=255)
                                   aaGATACATGCTGCAGTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTg                                                                                                                    <  2:510129/150‑1 (MQ=255)
                                            gctgcATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGCCGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCCGGCACACAACACGACACACGTCTGGCTACCgcg                                                                                                           >  2:2222114/1‑148 (MQ=255)
                                             ctgcATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGc                                                                                                          >  1:1061895/1‑150 (MQ=255)
                                              tgcATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGGGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTTCCGAgga                                                                                                         <  2:1851634/150‑3 (MQ=255)
                                              tgcATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCa                                                                                                         <  2:976141/150‑1 (MQ=255)
                                                 aTTTGATGCCGATCGGCGTCACCAAACAAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGa                                                                                                      >  2:1998811/1‑150 (MQ=255)
                                                     gATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGc                                                                                                  <  1:884484/150‑1 (MQ=255)
                                                       tGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAg                                                                                                >  2:2203245/1‑150 (MQ=255)
                                                       tGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAg                                                                                                >  1:1748771/1‑150 (MQ=255)
                                                        gCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGt                                                                                               >  2:603388/1‑150 (MQ=255)
                                                         ccGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTc                                                                                              >  1:1678624/1‑150 (MQ=255)
                                                           gATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCTGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCt                                                                                             <  2:2056201/149‑1 (MQ=255)
                                                                 cGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTaaaaa                                                                                      <  2:1982887/150‑1 (MQ=255)
                                                                   tCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAag                                                                                    <  1:1833001/150‑1 (MQ=255)
                                                                   tCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAag                                                                                    <  1:1724810/150‑1 (MQ=255)
                                                                    cACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAaga                                                                                   >  1:2106064/1‑150 (MQ=255)
                                                                    cACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAaga                                                                                   >  2:651475/1‑150 (MQ=255)
                                                                     accaaaccaaGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAagag                                                                                  <  1:35086/150‑1 (MQ=255)
                                                                      ccaaaccaaGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAagaga                                                                                 >  1:2041479/1‑150 (MQ=255)
                                                                           ccaaGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGa                                                                            >  1:709297/1‑150 (MQ=255)
                                                                              aGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAg                                                                         <  2:889195/150‑1 (MQ=255)
                                                                                   gAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCca                                                                    <  1:2203245/150‑1 (MQ=255)
                                                                                    aaCAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCcac                                                                   <  2:1767634/150‑1 (MQ=255)
                                                                                        gACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGAt                                                               >  1:918809/1‑150 (MQ=255)
                                                                                          cGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCCCAGATTg                                                             >  1:2458174/1‑150 (MQ=255)
                                                                                          cGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTg                                                             >  2:386371/1‑150 (MQ=255)
                                                                                          cGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTg                                                             >  2:395160/1‑150 (MQ=255)
                                                                                           gTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTgg                                                            <  1:2219953/150‑1 (MQ=255)
                                                                                              gCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGc                                                         >  1:1029298/1‑150 (MQ=255)
                                                                                               ccAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCa                                                        <  2:713876/150‑1 (MQ=255)
                                                                                                    cATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACa                                                   >  1:1706821/1‑150 (MQ=255)
                                                                                                       atCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATcc                                                >  2:1557915/1‑150 (MQ=255)
                                                                                                       atCGGCTGGTCAATCCATAACGCCAACGCGTCAACAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACCTCTGCCTACCGACCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATcc                                                >  2:460728/1‑150 (MQ=255)
                                                                                                         cGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCg                                              <  2:1009590/150‑1 (MQ=255)
                                                                                                          ggCTGGTCAATCAATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGc                                             >  2:338243/1‑150 (MQ=255)
                                                                                                              ggTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATaa                                         >  1:1271991/1‑150 (MQ=255)
                                                                                                              ggTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATaa                                         >  1:849376/1‑150 (MQ=255)
                                                                                                              ggTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATaa                                         >  1:170384/1‑150 (MQ=255)
                                                                                                                  aaTCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACg                                     >  2:2322041/1‑150 (MQ=255)
                                                                                                                    tCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACggg                                   <  2:1678624/150‑1 (MQ=255)
                                                                                                                    tCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACggg                                   <  1:1581567/150‑1 (MQ=255)
                                                                                                                      cATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACgggcc                                 >  2:214621/1‑148 (MQ=255)
                                                                                                                      cATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGaa                                 >  1:304713/1‑150 (MQ=255)
                                                                                                                          aCGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGac                             <  2:709297/150‑1 (MQ=255)
                                                                                                                           cGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGCGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTCGTGCAAAACATCCCGCATAAAACGGGACCGacc                            >  2:1196601/1‑149 (MQ=255)
                                                                                                                             ccAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACAcc                          >  2:1416969/1‑150 (MQ=255)
                                                                                                                                  gcgTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCgcc                     <  1:1486341/150‑1 (MQ=255)
                                                                                                                                  gcgTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCgcc                     >  1:1918801/1‑150 (MQ=255)
                                                                                                                                    gTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCgccag                   <  1:420668/150‑1 (MQ=255)
                                                                                                                                     tCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCgccagc                  <  1:1288729/150‑1 (MQ=255)
                                                                                                                                     tCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCgccagc                  <  1:2135614/150‑1 (MQ=255)
                                                                                                                                      ccccAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACAGGAACGACACCGCgccagcc                 >  1:2382437/1‑150 (MQ=255)
                                                                                                                                       gttgtgtGCCTGCGCAGTATTACGCTGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCgccagcca                >  2:2328153/7‑150 (MQ=255)
                                                                                                                                       cccAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCgccagcca                >  1:1036002/1‑150 (MQ=255)
                                                                                                                                       cccAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCgccagcca                >  1:1805842/1‑150 (MQ=255)
                                                                                                                                       cccAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCgccagcca                >  1:2252612/1‑150 (MQ=255)
                                                                                                                                         cAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCGCCAGCCaaa              <  2:944226/150‑1 (MQ=255)
                                                                                                                                         cAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCGCCAGCCaaa              <  1:2284551/150‑1 (MQ=255)
                                                                                                                                          gtgtGCCTGCGCAATATTACGCTGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCGCCAGCCAAAg             >  1:103880/4‑150 (MQ=255)
                                                                                                                                          aGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCGCCAGCCAAAg             <  2:1911335/150‑1 (MQ=255)
                                                                                                                                           tgtGCCTGCGCAATATTACGCTGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCGCCAGCCAAAGa            >  1:422671/3‑150 (MQ=255)
                                                                                                                                           tgtGCCTGCGCAATATTACGCTGGCACACAACACGACAAACGTCTAGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACcgcg                       >  1:1978719/3‑139 (MQ=255)
                                                                                                                                           tgtGCCTGCGCAATATTACGCTGGCACACAACACGACAAACGTCTAGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACcgcg                       >  1:1978726/3‑139 (MQ=255)
                                                                                                                                           tgtGCCTGCGCAATATTACGCTGGCACACAACACGACAAACGTCTAGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACcgcg                       <  2:1978726/137‑1 (MQ=255)
                                                                                                                                           tgtGCCTGCGCAATATTACGCTGGCACACAACACGACAAACGTCTAGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACcgcg                       <  2:1978719/137‑1 (MQ=255)
                                                                                                                                           gCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCGCCAGCCAAAGa            >  2:307851/1‑150 (MQ=255)
                                                                                                                                           gCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCGCCAGCCAAAGa            <  1:1416969/150‑1 (MQ=255)
                                                                                                                                            cTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCGCCAGCCAAAGAt           <  1:1998811/150‑1 (MQ=255)
                                                                                                                                             tGCCTGCGCAATATTACGCTGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCGCCAGCCAAAGATa          >  1:1607142/1‑150 (MQ=255)
                                                                                                                                              gCCTGCGCAATATTACGCTGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCGCCAGCCAAAGATaa         >  2:768909/1‑150 (MQ=255)
                                                                                                                                                cAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCGCCAGCCAAAGATAAAc       <  2:170384/150‑1 (MQ=255)
                                                                                                                                                 aGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCGCCAGCCAAAGATAAACa      >  2:1446957/1‑150 (MQ=255)
                                                                                                                                                   cGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCGCCAGCCAAAGATAAACAAg    <  2:2123336/150‑1 (MQ=255)
                                                                                                                                                   cGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCGCCAGCCAAAGATAAACAAg    <  2:1198590/150‑1 (MQ=255)
                                                                                                                                                    gTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCGCCAGCCAAAGATAAACAAGa   <  2:932301/150‑1 (MQ=255)
                                                                                                                                                     tAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCGCCAGCCAAAGATAAACAAGAt  >  2:1489528/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GATATCTTTAGCGGGAAATCAAAGAAGCTCCCCACAAGATACATGCTGCATTTGATGCCGATCGGCGTCACCAAACCAAGCCAGAACAGACGTTGCCAGACATATCGGCTGGTCAATCCATAACGCCAACGCGTCCCCAGCTGCCAGCGTAATATTGCGCAGGCACACAACACGACAAACGTCTGGCTACCGAGCATGACCGCAGTCTGTAAAAAAGAGAGATGAAAGTTCCACAGATTGGTGCAAAACATCCCGCATAAAACGGGAACGACACCGCGCCAGCCAAAGATAAACAAGAT  >  NC_000913/2517372‑2517670

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: