Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,362,263 C→G 8.2% V49L (GTG→CTG)  tpiA ← triose phosphate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,362,2630CG 8.2% 266.1 / 12.6 85V49L (GTG→CTG) tpiAtriose phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base C (39/39);  new base G (3/4);  total (42/43)
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.43e-01

GCTCCGAATGGCCAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCATAGGGCGACGCAT  >  NC_002947/5362116‑5362407
                                                                                                                                                    |                                                                                                                                                 
gCTCCGAATGGCCAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCt                                                                                                                                                 <  1:947927/150‑1 (MQ=255)
  tCCGAATGGCCAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCtga                                                                                                                                               >  2:1028997/1‑150 (MQ=255)
    cGAATGGCCAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCtgatt                                                                                                                                             <  2:1258567/150‑1 (MQ=255)
     gAATGGCCAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCtgattg                                                                                                                                            >  2:1180975/1‑150 (MQ=255)
       aTGGCCAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCtgattgat                                                                                                                                          <  1:1028997/150‑1 (MQ=255)
        tGGCCAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCAGGTGATtgatg                                                                                                                                         <  2:422902/150‑1 (MQ=255)
        tGGCCAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATtgatg                                                                                                                                         <  2:2165278/150‑1 (MQ=255)
          gCCAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCtgattgat                                                                                                                                          >  2:998434/1‑147 (MQ=255)
          gCCAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCtgattgat                                                                                                                                          <  1:998434/147‑1 (MQ=255)
           ccAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCAGGTGATtgatggcc                                                                                                                                      <  1:1130347/150‑4 (MQ=255)
           ccAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAAc                                                                                                                                      >  2:82132/1‑150 (MQ=255)
             aaTCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGTCATCAATCAGGTGATtgatggcct                                                                                                                                     <  2:707816/149‑5 (MQ=255)
             aaTCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGTCATCAATCAGGTGATtgatggcct                                                                                                                                     >  1:707816/1‑145 (MQ=255)
              aTCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGCTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAg                                                                                                                                   <  1:1000825/150‑1 (MQ=255)
                  acacCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCg                                                                                                                               <  2:1178148/150‑1 (MQ=255)
                      cAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTgg                                                                                                                           >  1:1773604/1‑150 (MQ=255)
                      cAACTTGCAACCTGATTCAACCAGCTGCTCGGGAGCAACTTCCCCGTTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTgg                                                                                                                           >  2:1421483/1‑150 (MQ=255)
                           tGCAACCTCCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGCCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGCGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCACTCACCTGATTGATGAACAAGGCCGGTGGAAATa                                                                                                                      <  2:132865/150‑1 (MQ=255)
                             cAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTGCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACg                                                                                                                    <  2:2337555/150‑1 (MQ=255)
                              aaCCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACgg                                                                                                                   >  2:47309/1‑150 (MQ=255)
                               aCCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGc                                                                                                                  <  2:1773604/150‑1 (MQ=255)
                                  tGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGAc                                                                                                               >  1:2397201/1‑150 (MQ=255)
                                   gCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACt                                                                                                              >  1:746345/1‑150 (MQ=255)
                                    cTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAAGCCATCAATCACCTGACTGATAAACAAGGCCGTTGGAAATACGGCGACtt                                                                                                             >  2:162177/1‑150 (MQ=255)
                                     ttCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTc                                                                                                            >  2:1635221/1‑150 (MQ=255)
                                     ttCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTc                                                                                                            >  2:768310/1‑150 (MQ=255)
                                       cAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCa                                                                                                          <  1:155177/150‑1 (MQ=255)
                                        aaCCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTTTTTGAACGGCCTGGGATGAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCAc                                                                                                         >  1:2090057/1‑150 (MQ=255)
                                         aCCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACt                                                                                                        >  2:584461/1‑150 (MQ=255)
                                          ccAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACgg                                                                                                                   >  1:108965/1‑138 (MQ=255)
                                          ccAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACgg                                                                                                                   <  2:108965/138‑1 (MQ=255)
                                           cAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTc                                                                                                       <  2:479848/149‑1 (MQ=255)
                                           cAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTc                                                                                                       >  1:479848/1‑149 (MQ=255)
                                            aGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCg                                                                                                     >  1:316772/1‑150 (MQ=255)
                                                gctCGGGAGCAACTTCCCCGTTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATGCGGCGACTTCCACTCCGCTCg                                                                                                 <  1:1421483/150‑1 (MQ=255)
                                                    gggAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTg                                                                                                                            <  1:650697/119‑1 (MQ=255)
                                                    gggAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTg                                                                                                                            >  2:650697/1‑119 (MQ=255)
                                                            cTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGAt                                                                                    <  2:746345/150‑1 (MQ=255)
                                                             ttCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCAGGTGATtgatg                                                                                                                                         <  2:2303125/97‑1 (MQ=255)
                                                             ttCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCAGGTGATtgatg                                                                                                                                         >  1:2303125/1‑97 (MQ=255)
                                                             ttCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGAtt                                                                                   >  2:114603/1‑150 (MQ=255)
                                                               ccccGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTAc                                                                                 >  1:1020077/1‑150 (MQ=255)
                                                                 ccGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCaa                                                                                            >  1:306452/1‑137 (MQ=255)
                                                                 ccGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCaa                                                                                            <  2:306452/137‑1 (MQ=255)
                                                                 ccGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTc                                                                               >  1:115365/1‑150 (MQ=255)
                                                                   ggTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACTAGGCGCGTCTAAAGACGGCTACTTCCACTCCGCTCGGCAAGGCGAGATTACTCaa                                                                             <  1:47309/150‑1 (MQ=255)
                                                                   ggTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCaa                                                                             <  2:1906843/150‑1 (MQ=255)
                                                                          gcACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTg                                                                      <  1:82132/150‑1 (MQ=255)
                                                                            aCCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGt                                                                    >  2:2151877/1‑150 (MQ=255)
                                                                             cccTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGtc                                                                   <  1:1673650/150‑1 (MQ=255)
                                                                              ccTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGtca                                                                  >  1:1031682/1‑150 (MQ=255)
                                                                                      ggTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTc                                                                                                            <  1:714159/101‑1 (MQ=255)
                                                                                      ggTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTc                                                                                                            >  2:714159/1‑101 (MQ=255)
                                                                                            aCAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTa                                                    <  2:1402029/150‑1 (MQ=255)
                                                                                             cagcagAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTgg                                                                     <  1:1261163/132‑1 (MQ=255)
                                                                                             cagcagAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTgg                                                                     >  2:1261163/1‑132 (MQ=255)
                                                                                              agcagAATTCTGTGCGCCTACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACt                                                                                                        >  2:1844120/1‑97 (MQ=255)
                                                                                              agcagAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACt                                                                                                        <  1:1844120/97‑1 (MQ=255)
                                                                                                  gAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCgg                                              <  2:2090057/150‑1 (MQ=255)
                                                                                                  gAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCgg                                              <  1:114603/150‑1 (MQ=255)
                                                                                                   aaTTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTg                                                                      <  1:2150379/125‑1 (MQ=255)
                                                                                                   aaTTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTg                                                                      >  2:2150379/1‑125 (MQ=255)
                                                                                                        tgtgCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCg                                        >  2:1620396/1‑150 (MQ=255)
                                                                                                             gcCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGCCAAACGGCAAACCGTTGAAGATCGGGTGAGCTCAGCGACGCTAGCGCGGGTACCGTGCAt                                   >  1:2341405/1‑150 (MQ=255)
                                                                                                                 aCAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGa                                                         >  1:281037/1‑124 (MQ=255)
                                                                                                                 aCAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGa                                                         <  2:281037/124‑1 (MQ=255)
                                                                                                                 aCAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTc                               <  2:115365/150‑1 (MQ=255)
                                                                                                                   aGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCa                             <  2:1020077/150‑1 (MQ=255)
                                                                                                                   aGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCa                             <  1:1914465/150‑1 (MQ=255)
                                                                                                                    gTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCAttt                                 >  2:1326205/1‑145 (MQ=255)
                                                                                                                    gTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCAttt                                 <  1:1326205/145‑1 (MQ=255)
                                                                                                                     tAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTg                                                                      <  2:2042336/107‑1 (MQ=255)
                                                                                                                     tAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTg                                                                      >  1:2042336/1‑107 (MQ=255)
                                                                                                                            tttACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGCCTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGc                    >  1:1591306/1‑150 (MQ=255)
                                                                                                                            tttACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGc                    >  2:1952789/1‑150 (MQ=255)
                                                                                                                                  cTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGAt                                                                                    >  1:1280555/1‑81 (MQ=255)
                                                                                                                                  cTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGAt                                                                                    <  2:1280555/81‑1 (MQ=255)
                                                                                                                                      cAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCATAgg           <  2:2147802/150‑1 (MQ=255)
                                                                                                                                      cAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAACGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCATAgg           <  1:1620396/150‑1 (MQ=255)
                                                                                                                                        ggCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCATAGGGc         <  1:768310/150‑1 (MQ=255)
                                                                                                                                        ggCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTCACCATCGAGCGTAGGGc         >  2:441750/1‑150 (MQ=255)
                                                                                                                                         ttcATCAATCAGGTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCATAGGGCg        >  1:2385677/3‑150 (MQ=255)
                                                                                                                                               aatcaCCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCaa                                                                             <  2:1444672/75‑1 (MQ=255)
                                                                                                                                               aatcaCCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCaa                                                                             >  1:1444672/1‑75 (MQ=255)
                                                                                                                                               aatcaCCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCATAGGGCGACGCAt  >  1:724088/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GCTCCGAATGGCCAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCATAGGGCGACGCAT  >  NC_002947/5362116‑5362407

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: