Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,967,035 A→G 5.7% F467S (TTC→TCC)  ptsP ← PEP‑protein phosphotransferase enzyme I; GAF domain containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,967,0350AG 5.7% 227.5 / 2.9 71F467S (TTC→TCC) ptsPPEP‑protein phosphotransferase enzyme I; GAF domain containing protein
Reads supporting (aligned to +/- strand):  ref base A (33/33);  new base G (2/2);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

CGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGC  >  NC_000913/2966892‑2967180
                                                                                                                                               |                                                                                                                                                 
cGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGaa                                                                                                                                                  <  5:53603/145‑1 (MQ=255)
 gCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATgg                                                                                                                                               <  3:75986/147‑1 (MQ=255)
     aTCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCTTACGCAAGGTGACGGGTTTATCATTGAACATTTGAAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGAtt                                                                                                                                           <  1:255732/147‑1 (MQ=255)
      tCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTc                                                                                                                                         <  6:30402/148‑1 (MQ=255)
      tCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTc                                                                                                                                         <  7:222938/148‑1 (MQ=255)
            tGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCg                                                                                                                                   <  1:59633/148‑1 (MQ=255)
                aGGCAGCTTCTTATCTGCTCCGACATCAAGCGTACTCAAGGTGAAGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCTACGGGAAACCACTTTGCAGCATTAATGGGATTTCAGTGCGATaa                                                                                                                               <  3:237527/148‑1 (MQ=255)
                aGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATaa                                                                                                                               <  6:263615/148‑1 (MQ=255)
                      ctgctTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAAt                                                                                                                                                 <  5:49114/124‑1 (MQ=255)
                      ctgctTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAAt                                                                                                                                                 >  6:49114/1‑124 (MQ=255)
                      ctgctTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGAcc                                                                                                                          >  7:148790/1‑147 (MQ=255)
                      ctgctTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCg                                                                                                                         <  1:287262/148‑1 (MQ=255)
                           tATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGc                                                                                                                     >  7:295201/1‑147 (MQ=255)
                                 cTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCa                                                                                                                <  1:2955/146‑1 (MQ=255)
                                   ccGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATAc                                                                                                            >  7:200949/1‑148 (MQ=255)
                                    cGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATa                                                                                                             <  8:324085/146‑1 (MQ=255)
                                     gACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACgg                                                                                                          <  8:295201/148‑1 (MQ=255)
                                          ccAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGcc                                                                                                     >  1:277089/1‑148 (MQ=255)
                                           cAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGcc                                                                                                     >  8:302401/1‑147 (MQ=255)
                                              cGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGt                                                                                                 <  7:302401/148‑1 (MQ=255)
                                                  cGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGttttt                                                                                             >  3:111039/1‑148 (MQ=255)
                                                     aatgTGACGTGTTTATCATTGAACATTTGCAGCATCCACTGATACTGCGCCACCGGTTCTTGTTGCGACGGGAAACCACTTTGCAGCTTGAATGGGATTTCAGTGCGATAAAGTCCGATGCCATCAATACGGCTGCCCAGTTTttctt                                                                                          <  6:156294/145‑1 (MQ=255)
                                                      aGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTttctt                                                                                          >  1:318659/1‑147 (MQ=255)
                                                      aGGTGACGGGGTTATCATTGAACAATTGCAGCATCCCCTGATACTGCGCCCCCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTttcttc                                                                                         <  8:116842/148‑1 (MQ=255)
                                                       ggTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCa                                                                                        <  8:200949/148‑1 (MQ=255)
                                                        gTGACGGGTTTATCAATGAACATATGCAGCATCCACTGATACAGCGCAACATGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACAGATGCAATCAAAACGGCTGCCCAGTTTATcatcat                                                                                       >  6:216601/1‑148 (MQ=255)
                                                               gTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCgg                                                                                <  6:315156/148‑1 (MQ=255)
                                                                tttATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCggg                                                                               >  2:131328/1‑148 (MQ=255)
                                                                 ttATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAAt                                                                                                              >  3:28775/1‑116 (MQ=255)
                                                                 ttATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAAt                                                                                                              <  4:28775/116‑1 (MQ=255)
                                                                      tttGAACATTTGCAGCATCCCCTGATTCTGCTCCACCTGTTCTTCTTCCGACGGGAAAACACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTaaa                                                                         <  5:216601/147‑1 (MQ=255)
                                                                      aTTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTaaa                                                                         >  1:85118/1‑148 (MQ=255)
                                                                                        ccccTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGc                                                             <  1:7337/142‑1 (MQ=255)
                                                                                        ccccTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGc                                                             >  2:7337/1‑142 (MQ=255)
                                                                                        ccccTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACt                                                       >  7:31805/1‑148 (MQ=255)
                                                                                           cTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACtttt                                                    <  8:31805/148‑1 (MQ=255)
                                                                                           cTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACttt                                                     >  5:172634/1‑147 (MQ=255)
                                                                                           cTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACttt                                                     >  2:99533/1‑147 (MQ=255)
                                                                                            tGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGAAAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTta                                                   >  5:62427/1‑148 (MQ=255)
                                                                                                    cgcCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTTATACGCTCa                                           >  1:230297/1‑148 (MQ=255)
                                                                                                            gTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGGATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTTATACGCTCACCGCtttt                                   >  8:150940/1‑148 (MQ=255)
                                                                                                            gTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGGATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATAGAGCATGACTTTAATACGCTCACCGCttat                                   >  8:146372/1‑146 (MQ=255)
                                                                                                            gTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTTATACGCTCACCGCtttt                                   >  6:84320/1‑148 (MQ=255)
                                                                                                             ttcttcttcCGACGGGAAACCACTTTGCAGCATGGATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTTATACGCTCACCGCtttt                                   <  7:150940/147‑1 (MQ=255)
                                                                                                             ttcttcttcCGACGGGAAACCACTTTGCAGCATGGATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTTATACGCTCACCGCtttt                                   <  7:146372/147‑1 (MQ=255)
                                                                                                             ttcttcttcCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTTATACGCTCACCGCtttt                                   <  8:96295/147‑1 (MQ=255)
                                                                                                              tcttcttcCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTTATACGCTCACCGCttttt                                  >  3:11022/1‑147 (MQ=255)
                                                                                                                 tcttcCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGc                                                             <  8:12080/117‑1 (MQ=255)
                                                                                                                 tcttcCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGc                                                             >  7:12080/1‑117 (MQ=255)
                                                                                                                     ccGACGGGAAACCCCTTTGCAGCATGAATGGGATTTCAGTGCGTTAAAGACCGATGCCATCCATTCGGCTGCGCAGTTTTTCTTCATGTTGCGGGCTTAACCCAGCATTGAGCATGACTTTTATACGCTCACGGCTTTTTAACGgctc                          >  4:172850/1‑146 (MQ=255)
                                                                                                                     ccGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTTATACGCTCACCGCTTTTTAACTg                             <  3:167224/145‑1 (MQ=255)
                                                                                                                      cGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTTATACGCTCACCGCTTTTTAACTGcgcg                         >  8:240914/1‑148 (MQ=255)
                                                                                                                       gACGGGAAACCACTTTGCAGCATGTATGGGATTTCAGTGAGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCTGCATTGAGCATGACTTTTATACGCTCACCGCTTTTTAACTGCGCgg                        <  6:62427/148‑1 (MQ=255)
                                                                                                                         cGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTTATACGCTCACCGCTTTTTa                                 <  5:291232/137‑1 (MQ=255)
                                                                                                                         cGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTTATACGCTCACCGCTTTTTa                                 >  6:291232/1‑137 (MQ=255)
                                                                                                                         cGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTTATACGCTCACCGCTTTTTAACTGCGCggg                       <  4:111039/147‑1 (MQ=255)
                                                                                                                            gAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTTATACGCTCACCGCTTTTTAACTGCGCGGGTaaa                   <  5:84320/148‑1 (MQ=255)
                                                                                                                             aaaCCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAAt                  <  6:300899/148‑1 (MQ=255)
                                                                                                                               aCCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTg                <  2:85118/148‑1 (MQ=255)
                                                                                                                               aCCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTg                <  4:251643/148‑1 (MQ=255)
                                                                                                                                ccACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACtttt                                                    <  4:35726/111‑1 (MQ=255)
                                                                                                                                ccACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACtttt                                                    >  3:35726/1‑111 (MQ=255)
                                                                                                                                ccACATTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCAAACAAAACGGCTGCCAAGTATATCTTCATGTTACTGGCTTAAACCAGCAATGAGCATGACTTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGa               >  6:166230/1‑148 (MQ=255)
                                                                                                                                  aCTTTGCAGCATGAATGGGATTTCAGT‑‑GATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGt            >  4:175578/1‑147 (MQ=255)
                                                                                                                                     ttGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCCTCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCCTGCATGTTATACGCTCACCGCTTTTAAACTGCGCGGGTAAATTGACGTc           >  4:200703/1‑147 (MQ=255)
                                                                                                                                            aTGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCg   >  1:288680/1‑148 (MQ=255)
                                                                                                                                             tGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGc  >  1:208817/1‑148 (MQ=255)
                                                                                                                                              gAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTtat                                                  >  7:332447/1‑99 (MQ=255)
                                                                                                                                              gAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTtat                                                  <  8:332447/99‑1 (MQ=255)
                                                                                                                                              gAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTTATACGCTCAGCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGc  >  6:276661/1‑147 (MQ=255)
                                                                                                                                              gAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTTATACGCTCACCGCTTTTTAACTGCGCGGGTaa                    <  7:118840/129‑1 (MQ=255)
                                                                                                                                              gAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTTATACGCTCACCGCTTTTTAACTGCGCGGGTaa                    >  8:118840/1‑129 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
CGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGC  >  NC_000913/2966892‑2967180

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