Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,477,875 T→A 5.0% Y226F (TAT→TTT) ‡ nadD ← nicotinate‑nucleotide adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,477,8750TA 5.0% 796.5 / 20.0 234Y226F (TAT→TTT) ‡nadDnicotinate‑nucleotide adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base T (132/90);  new base A (4/8);  total (136/98)
Fisher's exact test for biased strand distribution p-value = 1.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 242 total aligned reads displayed.

GATGACCTGGATGTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGAC  >  NC_002947/5477726‑5478024
                                                                                                                                                     |                                                                                                                                                     
gATGACCTGGATGTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGat                                                                                                                                                       >  2:815055/1‑150 (MQ=255)
   gACCTGGATGTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGAtata                                                                                                                                                    <  1:3795837/150‑1 (MQ=255)
   gACCTGGATGTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGAtata                                                                                                                                                    <  2:2828970/150‑1 (MQ=255)
      cTGGATGTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGtt                                                                                                                                                 >  2:1817141/1‑150 (MQ=255)
        ggATGTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCg                                                                                                                                               <  2:2077232/150‑1 (MQ=255)
          aTGTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTg                                                                                                                                             <  1:2707328/150‑1 (MQ=255)
           tGTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTgc                                                                                                                                            >  2:5349621/1‑150 (MQ=255)
               cTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCt                                                                                                                                        >  2:1883268/1‑150 (MQ=255)
                 gggCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGAACTATATCGTGCGcccaa                                                                                                                                      <  1:3935902/150‑4 (MQ=255)
                 gggCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGGTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGAACTATATCGTGCGcccaa                                                                                                                                      <  2:4657651/150‑4 (MQ=255)
                  ggCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGa                                                                                                                                     >  1:1527535/1‑150 (MQ=255)
                    ccTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGata                                                                                                                                   >  2:4870749/1‑150 (MQ=255)
                       tGACGTCTTCCAGCGCTGCCTTCGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAg                                                                                                                                <  1:1845524/150‑1 (MQ=255)
                           gTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCa                                                                                                                            <  2:3319655/150‑1 (MQ=255)
                           gTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCa                                                                                                                            <  1:1991165/150‑1 (MQ=255)
                              ttCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTa                                                                                                                         >  2:5076197/1‑150 (MQ=255)
                               tCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTATTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTAc                                                                                                                        >  1:4560038/1‑150 (MQ=255)
                                 cAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACgg                                                                                                                      >  2:1643527/1‑150 (MQ=255)
                                  aGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGc                                                                                                                     >  1:3027643/1‑150 (MQ=255)
                                   gcgcTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCg                                                                                                                    >  2:5406185/1‑150 (MQ=255)
                                   gcgcTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATACGCCAGTACGGCg                                                                                                                    <  2:648385/150‑1 (MQ=255)
                                       tGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCg                                                                                                                <  1:4870749/150‑1 (MQ=255)
                                        gCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCgg                                                                                                               >  1:5418845/1‑150 (MQ=255)
                                         ccTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGc                                                                                                              >  1:1748360/1‑150 (MQ=255)
                                           ttGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCAc                                                                                                            >  1:1365118/1‑150 (MQ=255)
                                           ttGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCAc                                                                                                            >  1:2492596/1‑150 (MQ=255)
                                           ttGGTCACTGCGACCAGTTCGTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCAc                                                                                                            <  2:4265622/150‑1 (MQ=255)
                                             ggTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCa                                                                                                          <  2:1533417/150‑1 (MQ=255)
                                             ggTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCa                                                                                                          <  2:474969/150‑1 (MQ=255)
                                              gTCACTGCGACCAGTTCGTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAg                                                                                                         <  2:2900302/150‑1 (MQ=255)
                                                   tGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAAcc                                                                                                    >  1:2720827/1‑150 (MQ=255)
                                                    gCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCt                                                                                                   >  2:1176186/1‑150 (MQ=255)
                                                    gCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCt                                                                                                   <  1:1720145/150‑1 (MQ=255)
                                                       aCCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCAc                                                                                                >  1:627318/1‑150 (MQ=255)
                                                       aCCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCAc                                                                                                <  2:444795/150‑1 (MQ=255)
                                                        ccAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCAcc                                                                                               >  1:3241177/1‑150 (MQ=255)
                                                        ccAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCAcc                                                                                               >  2:1485412/1‑150 (MQ=255)
                                                               ttcGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGc                                                                                        >  2:5620389/1‑150 (MQ=255)
                                                                 cGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCg                                                                                      <  1:1358324/150‑1 (MQ=255)
                                                                   ccGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCt                                                                                    <  2:4466922/150‑1 (MQ=255)
                                                                      tAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGc                                                                                 >  2:1960401/1‑150 (MQ=255)
                                                                      tAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGc                                                                                 <  1:2072084/150‑1 (MQ=255)
                                                                      tAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGc                                                                                 >  1:3292019/1‑150 (MQ=255)
                                                                         aTTTTCTGCTTGGTCATATAACACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCcag                                                                              <  2:2936663/150‑1 (MQ=255)
                                                                         aTTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCcag                                                                              <  1:4293074/150‑1 (MQ=255)
                                                                            ttCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCcagcag                                                                           <  1:671874/150‑1 (MQ=255)
                                                                                gCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGt                                                                       <  1:5566793/150‑1 (MQ=255)
                                                                                   tGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCgg                                                                    <  2:4100280/150‑1 (MQ=255)
                                                                                     gTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGACAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGAt                                                                  <  1:4034985/150‑1 (MQ=255)
                                                                                       cATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCCGAGGCACCGTCAGTTGGGCGCACGATCTAGTGCGTGCGCCTCGATATCGGCCAGTACGGCGTCCGGCACCAGGCACCGCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCt                                                                >  2:3562432/1‑150 (MQ=255)
                                                                                         tataAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATCTAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCtgt                                                              >  1:5504327/1‑150 (MQ=255)
                                                                                         tataAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCtgt                                                              <  1:1118181/150‑1 (MQ=255)
                                                                                         tataAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCtgt                                                              <  1:2540251/150‑1 (MQ=255)
                                                                                          ataAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCtgtg                                                             >  2:1037236/1‑150 (MQ=255)
                                                                                          ataAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCtgtg                                                             >  2:1403157/1‑150 (MQ=255)
                                                                                          ataAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCtgtg                                                             >  1:525410/1‑150 (MQ=255)
                                                                                          ataAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCACGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCtgtg                                                             >  2:4241764/1‑150 (MQ=255)
                                                                                           taAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTgtgt                                                            >  1:4573444/1‑150 (MQ=255)
                                                                                           taAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACAGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTgtgt                                                            >  2:2872929/1‑150 (MQ=255)
                                                                                              aaCTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGc                                                         >  2:4945086/1‑150 (MQ=255)
                                                                                              aaCTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGc                                                         >  1:3009274/1‑150 (MQ=255)
                                                                                               aCTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGAtt                                                                                           <  2:3446494/115‑1 (MQ=255)
                                                                                                 tCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGa                                                      <  2:1641382/150‑1 (MQ=255)
                                                                                                   gTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGaca                                                    >  1:4307134/1‑150 (MQ=255)
                                                                                                       aaCTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGc                                                >  1:4660734/1‑150 (MQ=255)
                                                                                                       aaCTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGATCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGc                                                >  1:5023691/1‑150 (MQ=255)
                                                                                                        aCTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCa                                               >  1:1080182/1‑150 (MQ=255)
                                                                                                         cTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCaa                                              <  1:4809531/150‑1 (MQ=255)
                                                                                                          tCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAg                                             >  2:2759973/1‑150 (MQ=255)
                                                                                                           cGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAgc                                            <  1:5044886/150‑1 (MQ=255)
                                                                                                            gTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAgcg                                           >  2:66100/1‑150 (MQ=255)
                                                                                                            gTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAgcg                                           <  2:5571303/150‑1 (MQ=255)
                                                                                                            gTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAgcg                                           <  1:4016120/150‑1 (MQ=255)
                                                                                                            gTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAgcg                                           >  1:3465496/1‑150 (MQ=255)
                                                                                                              tGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAgcggc                                         <  1:1187247/150‑1 (MQ=255)
                                                                                                                gATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGt                                       >  2:2209906/1‑150 (MQ=255)
                                                                                                                   gAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTg                                    >  1:3698189/1‑150 (MQ=255)
                                                                                                                     ggCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGcc                                  >  2:4856540/1‑150 (MQ=255)
                                                                                                                     ggCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGcc                                  >  1:3830416/1‑150 (MQ=255)
                                                                                                                     ggCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGcc                                  <  2:2511728/150‑1 (MQ=255)
                                                                                                                     ggCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGcc                                  >  2:2451004/1‑150 (MQ=255)
                                                                                                                     ggCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGcc                                  >  1:2199736/1‑150 (MQ=255)
                                                                                                                     ggCGCCAGAGGCACCGGCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGcc                                  >  2:5019046/1‑150 (MQ=255)
                                                                                                                      gcgcCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCa                                 >  1:3627337/1‑150 (MQ=255)
                                                                                                                      gcgcCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCa                                 >  1:4489800/1‑150 (MQ=255)
                                                                                                                        gcCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGa                               >  1:334710/1‑150 (MQ=255)
                                                                                                                        gcCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGa                               >  1:4904165/1‑150 (MQ=255)
                                                                                                                            gagGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGacgaa                           <  1:4633884/150‑1 (MQ=255)
                                                                                                                             agGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGacgaac                          <  2:4104200/150‑1 (MQ=255)
                                                                                                                                cACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGAcgaacgaa                       <  1:3222486/150‑1 (MQ=255)
                                                                                                                                 aCCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGaaa                      >  2:766138/1‑150 (MQ=255)
                                                                                                                                  ccGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAat                     >  1:3600629/1‑150 (MQ=255)
                                                                                                                                   cGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCACGCGGCGTCTGCCAGACGAACGAAata                    >  2:530127/1‑150 (MQ=255)
                                                                                                                                         cgaggCGCACGAACTATATCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGcc                                                                                       >  1:240717/4‑77 (MQ=37)
                                                                                                                                          tGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCccg             <  1:3176839/150‑1 (MQ=255)
                                                                                                                                            ggCGCACGAACTATATCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCccgcc           >  2:3210461/1‑150 (MQ=255)
                                                                                                                                              cgcACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCgccgg         >  2:4424384/1‑150 (MQ=255)
                                                                                                                                               gcACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCgccggg        >  1:2891539/1‑150 (MQ=255)
                                                                                                                                                 aCGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCgccgggcc      >  1:5159176/1‑150 (MQ=255)
                                                                                                                                                    atatAGGTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGa   <  2:3398305/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GATGACCTGGATGTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGAC  >  NC_002947/5477726‑5478024

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