Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 2,976,822 T→A 22.3% Y245F (TAT→TTT)  lysA ← diaminopimelate decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077792,976,8220TA22.3% 151.6 / 12.2 73Y245F (TAT→TTT) lysAdiaminopimelate decarboxylase
Reads supporting (aligned to +/- strand):  ref base T (34/20);  new base A (3/13);  total (40/33)
Fisher's exact test for biased strand distribution p-value = 3.46e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGATCGGTGTACCAG  >  NC_007779/2976588‑2977058
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                            
tCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCataata                                                                                                                                                                                                                                          <  1:524687/239‑1 (MQ=255)
             cACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGtt                                                                                                                                                                                                                                     <  1:2028493/231‑1 (MQ=255)
                     tggtggTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGTAATTAATACGCCCGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGGATTTCTCCCCCCCGGCAATCCCCAGAACAAAAAAATATTTGGGATTAACCGccc                                                                                                                                                                                                                      >  2:2113655/1‑237 (MQ=255)
                       gtggtAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCCCACCTGAGTAATTAATACGCCAGACCGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGGGGCCCAAAAGGCGGGCGGTTTGCTCACGCGCGGCCTTTCACAGACCCCAAAAATATTCGGGATCAACCCCttct                                                                                                                                                                                                                    >  2:2028493/1‑238 (MQ=255)
                          gtAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCCTTTGTTTGACGCTCCGCACCTGAGGAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCCAAGGGGGGCCCAAATGGCGGGGCAATTGCTCAACCCCGGCATTCCCCCCACCAAAATAAAGGTCGGGATCCACCGCCCCTTCa                                                                                                                                                                                                                 >  2:1771147/1‑238 (MQ=255)
                                           ccGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGAATTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCCGTTGATAAGGAACa                                                                                                                                                                                               >  1:448830/1‑239 (MQ=255)
                                           ccGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCCCCGCACCTGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCCATTTTCTCACCCCCGGCATTCCACAGACCCCAATAAAGTTCGGTATCCAACCCCCCTTCACCCTGTTTATAAGGAACa                                                                                                                                                                                               >  2:1697003/1‑239 (MQ=255)
                                              ggCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCCCCGCACCTGAGTAATTAATACGCCCGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGGGGCCCAAATGGCGGGCGAGTTGCTCACGCGCCGGCCTCCCCAGACCACAAAAATACTGGGGATCAACCGCCCCTTCCACCTGTTGTTAAGGAACaaaa                                                                                                                                                                                            >  2:1997525/1‑239 (MQ=255)
                                                     cAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGGGGCCCAAATGGCGGGCGAATTGCTCACCCCCGGCATTCCACAGACCAAAAAAATGTTCGGGATCAACCGCCTCCTCACCCTGTTGATAAGGAACCAAAAACCCa                                                                                                                                                                                      >  1:640987/1‑238 (MQ=255)
                                                          cGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCCCCCCATTTGTTTGACGCTCCGCACCTGATTAATTAATACGCCAGACTGCGCTCCCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAAc                                                                                                                                                                                                                           <  1:855908/196‑1 (MQ=255)
                                                               aaCCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGTAATTAATACGCCCGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCAccgccattgaa                                                                                                                                                                           >  1:57968/1‑233 (MQ=255)
                                                                   cgccATCAACCAGCACAAAGTGGCGGCCCCCCATTTGTTGGCCGCTCCGCACCTGGGTAATTAATCCGCCAGACTGCGCTACCGGGAGGGGCCCCGGTTCAATTTCCAGTTTCACAGGGTGGCCAAAAGGGCGGGCGATTTGCTCACGCGCGGCTTTCCACAGACAATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATa                                                                                                                                                                       <  2:986683/236‑1 (MQ=255)
                                                                      cATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTTCTCACGCCCGGCCTTCCACCGACCATAATAATGTTCGGGATCAAACGCCTCCTCACCCTCTTGTTAAGGAACAGAAAACCCACCCCCCGCAGAAATAGc                                                                                                                                                                     >  2:935718/1‑238 (MQ=255)
                                                                                 acaAAGTGGCGGCCCCCAATTTGTTTGACGCTCCGCCCCTGAGTAATTAATACGCCAGCCGGCGCTACCAGGAAGCGCCCCGGTTAAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATc                                                                                                                                                            <  1:1770570/236‑1 (MQ=255)
                                                                                     gggggcggcTCCCCATTTTTTTGACGCCCCGCCCCTGAGAAATTAATACCCCAGCCTGCGCTCCCAGGAAGCGACCCGGTTCATTTCCCAGTTCCACAGGGTGCCCCAAAGGGCGGGCGTTTTGCCCACGCGCGGCATCCCACAGACAAAAATAATGTTCGGTACCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGc                                                                                                                                                                     <  2:1560539/220‑1 (MQ=255)
                                                                                        ggcggcTCCCCATTTGTTTGCCGCTCCGCCCCTGAGTAATTAATCCGCCAGCCTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGCCCCAAATGGCGGGCGATTTGCTCCCGCGCGGCATTCCACAAACCAAAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGa                                                                                                                                                    <  1:1771147/237‑1 (MQ=255)
                                                                                               ccccATTTGTTTGACGCTCCGCAACTGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATcc                                                                                                                                                           >  1:1560539/1‑223 (MQ=255)
                                                                                                 ccATTTTTTGGCCGCTCCCCGCCGGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCAAAATGGCGGGCGATTTCCCCCCGCGCGGCATTCCACAGACCATAATAATGTCCGGTACCACCCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGAcc                                                                                                                                         <  2:299738/239‑1 (MQ=255)
                                                                                                        tttgACGCTCCGCACCTGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATACGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGc                                                                                                                                   >  1:299738/1‑238 (MQ=255)
                                                                                                            aCGCTCCGCACCTGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCAt                                                                                                                              >  1:931566/1‑239 (MQ=255)
                                                                                                              gCTCCGCACCTGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCACTTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCCCGCGCGGCATTCCACAGACCATAAAAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAg                                                                                                                            <  1:1221305/239‑1 (MQ=255)
                                                                                                                       cTGAGTAATTAATACCCCAGACTGCGCTCCCAGGAGGCGCCCCGGTTCAATTTCCAGTTTCACGGGGTGGCCAAAAGGGCGGGCGTTTTGCTCCCGCGCGGCATTCCACAGACCAAAATAATTTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCAccgcccg                                                                                                                                                                               <  2:2197811/179‑1 (MQ=255)
                                                                                                                          aGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGcaccaccacc                                                                                                                  >  2:390712/1‑233 (MQ=255)
                                                                                                                             atttaatTCGCCAGACTGCCCTACCAGGAAGCGACCCGTTTCAATTTCCAGTTACACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAAACCAAAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGtt                                                                                                             <  1:390712/237‑1 (MQ=255)
                                                                                                                              attaatACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATa                                                                                                                                                                       >  1:271955/1‑180 (MQ=255)
                                                                                                                                 aatACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGc                                                                                                                                                                     >  1:2197811/1‑179 (MQ=255)
                                                                                                                                 aatACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCacacac                                                                                                                  >  1:208694/1‑230 (MQ=255)
                                                                                                                                           aCTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGACCGCAGAAATAGCCTGTAAATCCTGACAGAATTCGATGACCTGACGCACCATAGCCCCACACACACGTCCCAGACGGGCATaa                                                                                               >  1:2041898/1‑239 (MQ=255)
                                                                                                                                           aCTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACCCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAACTAGCCTGTAAATCATGACCGAATTCGATGACCTGACGCACCATAGCACCACCCACCTGTTCCAGATGGGCATag                                                                                               >  1:647430/1‑238 (MQ=255)
                                                                                                                                                gcTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAAc                                                                                          >  1:1588202/1‑239 (MQ=255)
                                                                                                                                                  tACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCCCGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGc                                                                                        <  1:1697003/239‑1 (MQ=255)
                                                                                                                                                  acccgggagcggACCCGGTTAAATTCCCGGTTTCCCGGGGGGGCCAAAGGGGGGGGGGTTTTGCCCCCGCGGGGCACTCCACAGACCATAAAAATGTTGGGTACCACCCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGc                                                                                        <  2:2041898/228‑1 (MQ=255)
                                                                                                                                                    ccAGGAAGCGCCCCGTTTCAATTTCCAGTTCCCCGGGGGGGCCAAAATGGCGGGCGATTTGCTCCCGCGCGGCATTCCACAGACCAAAAAAATGTTCGGTACCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAg                                                                                                                                                                      <  2:1440226/159‑1 (MQ=255)
                                                                                                                                                    caaggaggCGCCCCGGTGCAATTCCCAGTTTCACAGGGTGGCCCAAAGGGCGGGCGTTTGCCCCCCGCGCGGCATCCCACAGACCAAAATAATGTCCGGTACCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGaaa                                                                                                                                                                         <  2:1044112/154‑1 (MQ=255)
                                                                                                                                                         aaGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTTCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCAACATAGCACCACCCACCTCATCCAAATGGGCAAAATACACGCCAgcaa                                                                                  >  2:1602609/1‑235 (MQ=255)
                                                                                                                                                              aCCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCAccc                                                                                                                                                                                                              >  1:991536/1‑109 (MQ=255)
                                                                                                                                                              aCCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCAccc                                                                                                                                                                                                              <  2:991536/109‑1 (MQ=255)
                                                                                                                                                              aCCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATc                                                                                                                                                            >  1:1440226/1‑159 (MQ=255)
                                                                                                                                                              aCCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAAtgt                                                                            >  1:1299670/1‑239 (MQ=255)
                                                                                                                                                              aCCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGAAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGa                                                                                                        >  1:1487231/1‑211 (MQ=255)
                                                                                                                                                              aCCCGGTTCAATTTCCAGTTTCACAGGGCGGCCCAAATGGAGGGCGAGTTGCTCACCCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCAGCAGAAATAACCTGTaa                                                                                                                                                               >  1:1044112/1‑156 (MQ=255)
                                                                                                                                                                  ggTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACcacac                                                                                                                    <  1:105929/195‑1 (MQ=255)
                                                                                                                                                                  ggTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACcacac                                                                                                                    >  2:105929/1‑195 (MQ=255)
                                                                                                                                                                              aGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACAcc                                                                                                                 <  1:537224/186‑1 (MQ=255)
                                                                                                                                                                              aGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGa                                                             >  2:77713/1‑238 (MQ=255)
                                                                                                                                                                              aGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGa                                                             >  1:1562754/1‑238 (MQ=255)
                                                                                                                                                                              aGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGGATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACAcc                                                                                                                 >  2:537224/1‑186 (MQ=255)
                                                                                                                                                                                 ttCACAGGGTGGCCCAAAGGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGAc                                                            <  1:77713/236‑1 (MQ=255)
                                                                                                                                                                                 tacccaGGGTGGCCAAAATGCCGGGCGATTTGCTCACGCCCGCCATCCCCCAGACAAAAATAATGTCCGGTACCACCCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGAc                                                            <  1:2113655/232‑1 (MQ=255)
                                                                                                                                                                                    acaGGGTGCCCAAAATGGCGGGCGATTTCCTCACGCGCGGCATCCCACAGACCATAAAAATGTTCGCTACCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTg                                                      <  1:1878463/239‑1 (MQ=255)
                                                                                                                                                                                     caGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCCGAACCAATGTGCATGTGAATGCCGACCAGCtac                                                     >  2:1215119/1‑237 (MQ=255)
                                                                                                                                                                                          tGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAAACAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAgatg                                                >  1:1775123/1‑239 (MQ=255)
                                                                                                                                                                                             cccAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGatgatg                                             >  1:2208869/1‑239 (MQ=255)
                                                                                                                                                                                             ccaaaggggggggcgATTGCCCCACGGGCGGCATCCAACAGACAATAATAATGTTCGGTACCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGatgatg                                             <  2:265967/229‑1 (MQ=255)
                                                                                                                                                                                             ccaaactGGCGGTCGATTTCCCCCCGCGCGCCTTTCAAAAAACAATAAAAAGGTCCGTTACCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAgatgat                                              <  2:1299670/232‑1 (MQ=255)
                                                                                                                                                                                                  aTGGCGGGCGATTTGCTCCCGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGtt                                        <  1:1881053/239‑1 (MQ=255)
                                                                                                                                                                                                  aTGGCGGGCATTTTCCACCCCCCCGGCTTCCCCCAGACCATAATAATGTCCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGtt                                        <  1:326284/239‑1 (MQ=255)
                                                                                                                                                                                                      cgggcgATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTa                                     >  1:265967/1‑238 (MQ=255)
                                                                                                                                                                                                      cgggcgATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGGTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATTGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCCACCAGCTGCAGATGATGACGTTGTa                                     >  2:1269687/1‑238 (MQ=255)
                                                                                                                                                                                                       gggcgATTTGCTCACGCGCGCCATTCCACAGCCCATAATAAGGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATc                                   <  1:229344/239‑1 (MQ=255)
                                                                                                                                                                                                       gggcattttGCTCACCCGGGCCTTCCCCCAGACCAAAAAAATGTGGGGTACCAGCCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTAt                                    <  2:583337/232‑1 (MQ=255)
                                                                                                                                                                                                        ggcgATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGa                                                                                                                                      >  2:533449/1‑139 (MQ=255)
                                                                                                                                                                                                        ggcgATTTGCTCACGCGCGGCATTCCACAGACCAAAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGa                                                                                                                                      <  1:533449/139‑1 (MQ=255)
                                                                                                                                                                                                             tttGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTcc                             >  2:1499720/1‑239 (MQ=255)
                                                                                                                                                                                                                             tttCCGCAAACCAAAATAATGTTCGGTATCACCCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGAt             <  1:307223/238‑1 (MQ=255)
                                                                                                                                                                                                                             ctccccaggACAATAATAAGTTCCGGTACAAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGAt             <  2:1775123/231‑1 (MQ=255)
                                                                                                                                                                                                                             atcccacgaaCCAAAATAATGTCCGTTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGAt             <  2:1588202/230‑1 (MQ=255)
                                                                                                                                                                                                                             aTTCCACAGCCCATAAAAATTTCGGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGAt             <  1:1484868/239‑1 (MQ=255)
                                                                                                                                                                                                                               accacaaccaataagaAGGTTCGTTACCACCCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACAGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGATCg           <  2:1422254/229‑1 (MQ=255)
                                                                                                                                                                                                                                 cacaGACCAAAAAAAGTTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGATCGgt         <  2:647430/239‑1 (MQ=255)
                                                                                                                                                                                                                                    aGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGATCGgtgt       >  2:429680/1‑238 (MQ=255)
                                                                                                                                                                                                                                      aCCATAATAATGTTCGGTATCAACCGCCTCTTCACTCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCa                                                                                                  >  2:596163/1‑145 (MQ=255)
                                                                                                                                                                                                                                      aCCAAAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCa                                                                                                  <  1:596163/145‑1 (MQ=255)
                                                                                                                                                                                                                                        cATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGGACGGGGTaaacg  >  2:46655/1‑235 (MQ=255)
                                                                                                                                                                                                                                        cATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGAAGATGATGACGTTGTATCACGTCCAGTGCGCCGGGCAGATCGGGGTACCaa  >  1:1745289/1‑238 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                            
TCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGATCGGTGTACCAG  >  NC_007779/2976588‑2977058

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