Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,323,197 T→C pseudogene (34/426 nt) phnE ← defective phosphonate ABC transporter permease;transport; Transport of small molecules: Anions; membrane channel protein component of Pn transporter; membrane component of ABC superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,323,1970TC99.2% 466.5 / ‑4.7 132pseudogene (34/426 nt)phnEdefective phosphonate ABC transporter permease;transport; Transport of small molecules: Anions; membrane channel protein component of Pn transporter; membrane component of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (52/78);  total (52/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

Only 100 of 132 total aligned reads displayed.

CGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCAC  >  NC_000913/4323101‑4323293
                                                                                                |                                                                                                
 gTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCgcgc                                                                                              <  1:1658344/100‑1 (MQ=255)
  tCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCt                                                                                             <  1:1016733/100‑1 (MQ=255)
  tCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCt                                                                                             <  1:787946/100‑1 (MQ=255)
   cTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCtt                                                                                            >  2:3055376/1‑100 (MQ=255)
   cTTTGATCAGCGTAAACGGAGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACAGCCCAGCTCAAAAGAGAGAACCAGCCgc                                                                                                >  2:1487201/1‑96 (MQ=255)
    tttGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTtg                                                                                           <  1:1435400/100‑1 (MQ=255)
     ttGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTtg                                                                                           <  2:309345/99‑1 (MQ=255)
     tggATCGGCGTGAGCGGGGCCAGTGCGGCGCCCGGCCACGAGGCGACCAACACGGCGAGTACAAGGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTtg                                                                                           <  2:2369195/97‑1 (MQ=255)
      tGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTtggg                                                                                         <  2:1884593/100‑1 (MQ=255)
      tGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTtggg                                                                                         >  1:683310/1‑100 (MQ=255)
        aTCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTtgggtg                                                                                       <  2:2548811/100‑1 (MQ=255)
        aTCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTtgggtg                                                                                       >  1:281495/1‑100 (MQ=255)
           aGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACTAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTgggg                                                                                    >  1:1643626/1‑100 (MQ=255)
           aGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTgggg                                                                                    <  2:2533547/100‑1 (MQ=255)
            gCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGc                                                                                   <  2:1211003/100‑1 (MQ=255)
            gCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGc                                                                                   >  1:2957296/1‑100 (MQ=255)
            gCGTAAGCGGGGCCATTTCCGCGCCCGGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGTTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGc                                                                                   <  1:2136479/100‑1 (MQ=255)
              gTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCga                                                                                 >  2:1816657/1‑100 (MQ=255)
               tAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCgat                                                                                >  1:1844544/1‑100 (MQ=255)
                 aGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCgatgg                                                                              >  2:827098/1‑100 (MQ=255)
                 aGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCgatgg                                                                              >  1:220324/1‑100 (MQ=255)
                 aGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCgatgg                                                                              <  1:2046423/100‑1 (MQ=255)
                  gCGGGGCCGTTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCgatggt                                                                             <  2:1609660/100‑1 (MQ=255)
                  gCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCgatggt                                                                             <  1:181558/100‑1 (MQ=255)
                  gCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCgatggt                                                                             <  2:1050588/100‑1 (MQ=255)
                    ggggCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCgatggtga                                                                           >  2:462945/1‑100 (MQ=255)
                    ggggCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCgatggtga                                                                           <  2:1880378/100‑1 (MQ=255)
                      ggCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCgatggtgatg                                                                         <  2:1990631/100‑1 (MQ=255)
                         cATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGtt                                                                      >  1:247099/1‑100 (MQ=255)
                         cATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGtt                                                                      <  2:2592281/100‑1 (MQ=255)
                          aTTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGttt                                                                     <  1:658362/100‑1 (MQ=255)
                           tttCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTg                                                                    <  1:2717900/100‑1 (MQ=255)
                           tttCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTg                                                                    <  2:167614/100‑1 (MQ=255)
                            ttcgccgcCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTgn                                                                   <  1:3003405/95‑2 (MQ=255)
                            ttCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGc                                                                   <  2:2644956/100‑1 (MQ=255)
                            ttCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGc                                                                   <  1:1382691/100‑1 (MQ=255)
                             tCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCa                                                                  >  1:1990256/1‑100 (MQ=255)
                              ccGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCAt                                                                 <  1:1668984/100‑1 (MQ=255)
                               cgcgCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGAAGGGTGGCATg                                                                >  2:800160/1‑100 (MQ=255)
                                   cccTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCGGCGCTTGGGTGGGGCGATGGTGATTGTTTGCATGTTTg                                                            <  2:1903287/100‑1 (MQ=255)
                                   cccTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGttag                                                            <  1:2861108/100‑3 (MQ=255)
                                   cccTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTg                                                            >  1:2807992/1‑100 (MQ=255)
                                   cccTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTg                                                            >  2:3082022/1‑100 (MQ=255)
                                   cccTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTg                                                            <  2:1367139/100‑1 (MQ=255)
                                     cTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGc                                                          <  1:599818/100‑1 (MQ=255)
                                     cTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGc                                                          >  1:2859218/1‑100 (MQ=255)
                                         cACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCg                                                      >  2:2850918/1‑100 (MQ=255)
                                          aCGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGTTGTTGGTTTGCATGTTTGGCTCCg                                                      <  2:755722/99‑1 (MQ=255)
                                          aCGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCgg                                                     >  2:2262280/1‑100 (MQ=255)
                                          aCGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCgg                                                     <  1:689367/100‑1 (MQ=255)
                                          aCGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCgg                                                     >  1:2897352/1‑100 (MQ=255)
                                            gagaCGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTc                                                   >  1:1353553/1‑100 (MQ=255)
                                            gagaCGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTc                                                   <  1:2977084/100‑1 (MQ=255)
                                             agaCGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCt                                                  >  1:2868053/1‑100 (MQ=255)
                                               acgacCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGt                                                <  2:28850/100‑1 (MQ=255)
                                                 gacCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAg                                              >  2:2782091/1‑100 (MQ=255)
                                                 gacCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAg                                              <  1:1043649/100‑1 (MQ=255)
                                                  acCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGTCGAGGGCGATGGTTTGCATGTTTGGCTCCGGTCTGTAgg                                             >  2:2164278/1‑100 (MQ=255)
                                                  acCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAgg                                             <  2:171420/100‑1 (MQ=255)
                                                  acCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAgg                                             <  2:2824909/100‑1 (MQ=255)
                                                    cAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGcc                                           >  2:1385962/1‑100 (MQ=255)
                                                      acacGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCgg                                         <  1:1472158/100‑1 (MQ=255)
                                                      acacGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCgg                                         <  2:1725445/100‑1 (MQ=255)
                                                      acacGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCgg                                         <  1:1390970/100‑1 (MQ=255)
                                                        acGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGAt                                       <  2:683307/100‑1 (MQ=255)
                                                          ggCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATa                                      <  2:1675481/99‑1 (MQ=255)
                                                             gAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGt                                  <  2:3041844/100‑1 (MQ=255)
                                                             gAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGt                                  <  2:2274348/100‑1 (MQ=255)
                                                              aGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGt                                  <  1:1487215/99‑1 (MQ=255)
                                                              aGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTACATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTg                                 <  2:114116/100‑1 (MQ=255)
                                                                 aCAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCt                               <  2:2722377/99‑1 (MQ=255)
                                                                  cAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTg                             <  2:647247/100‑1 (MQ=255)
                                                                  cAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTg                             <  2:605287/100‑1 (MQ=255)
                                                                  cAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTg                             <  1:1919918/100‑1 (MQ=255)
                                                                  cAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTg                             >  1:1220735/1‑100 (MQ=255)
                                                                  cAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTg                             >  1:802151/1‑100 (MQ=255)
                                                                   aaCGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGc                            <  2:619657/100‑1 (MQ=255)
                                                                      ggtccAGCGCAGAAGCGAGAACCAGCCGCGCTGGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCAcc                         <  2:2488067/97‑1 (MQ=255)
                                                                       gCCCAGGTCAGAAGCGAGAACCAGCCGCGCTGGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCg                        <  2:2856634/100‑1 (MQ=255)
                                                                       gCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGTTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCg                        <  1:1000175/100‑1 (MQ=255)
                                                                         ccAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCa                      <  2:378016/100‑1 (MQ=255)
                                                                         ccAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCa                      >  2:2387686/1‑100 (MQ=255)
                                                                              tCAGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCgg                 <  1:250220/100‑1 (MQ=255)
                                                                                aGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCa               >  1:2195080/1‑100 (MQ=255)
                                                                                aGAAGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCa               <  1:2100374/100‑1 (MQ=255)
                                                                                  aaGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATc             <  1:37139/100‑1 (MQ=255)
                                                                                  aaGCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATc             <  2:220323/100‑1 (MQ=255)
                                                                                    gCGAGAACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCaa           <  1:901353/100‑1 (MQ=255)
                                                                                      gagaACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCCACg         >  2:126221/1‑100 (MQ=255)
                                                                                      gagaACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACg         >  2:778086/1‑100 (MQ=255)
                                                                                      gagaACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACg         >  2:1245948/1‑100 (MQ=255)
                                                                                        gaACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGcc       <  1:2850945/100‑1 (MQ=255)
                                                                                        gaACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGcc       <  1:445942/100‑1 (MQ=255)
                                                                                        gaACCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGcc       >  1:748505/1‑100 (MQ=255)
                                                                                          aCCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTg     >  2:1422202/1‑100 (MQ=255)
                                                                                          aCCAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTg     >  2:2280476/1‑100 (MQ=255)
                                                                                           ccAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGc    <  2:2195057/100‑1 (MQ=255)
                                                                                           ccAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGc    >  2:228616/1‑100 (MQ=255)
                                                                                            cAGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGca   <  2:2859189/100‑1 (MQ=255)
                                                                                             aGCCGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGcac  >  2:2059101/1‑100 (MQ=255)
                                                                                                |                                                                                                
CGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCAC  >  NC_000913/4323101‑4323293

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