Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 895,829 G→A 7.6% R197C (CGC→TGC)  PP_0778 ← group 1 family glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947895,8290GA 7.6% 187.5 / 4.4 65R197C (CGC→TGC) PP_0778group 1 family glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base G (34/26);  new base A (3/2);  total (37/28)
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.90e-01

ACCATGATCAGGCGCAGGCGCTTTTGTGGATACGTCTTGCGTAATGCCTCGAAGCAAGGCCGTAGCAAGCCAAGGTTCTTTTCCGCCGCCAGCCGCCCGACATGCAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATCGACCCCTCGGGCCAGCAACTCCAGGCGTTCGAAACCGCGACGCTCCAGCTCCAGGCGCTGGCTGGCACTCGGCACCAAGGTGATTGCCGTGCGCCGATG  >  NC_002947/895684‑895970
                                                                                                                                                 |                                                                                                                                             
aCCATGATCAGGCGCAGGCGCTTTTGTGGATACGTCTTGCGTAATGCCTCGAAGCAAGGCCGTAGCAAGCCAAGGTTCTTTTCCGCCGCCAGCCGCCCGACATGCAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAgg                                                                                                                                           >  2:792364/1‑150 (MQ=255)
     gATCAGGCGCAGGCGCTTTTGTGGATACGTCTTGCGTAATGCCTCGAAGCAAGGCCGTAGCAAGCCAAGGTTCTTTTCCGCCGCCAGCCGCCCGACATGCAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAGGGCCTg                                                                                                                                      >  2:718534/1‑150 (MQ=255)
      aTCAGGCGCAGGCGCTTTTGTGGATACGTCTTGCGTAATGCCTCGAAGCAAGGCCGTAGCAAGCCAAGGTTCTTTTCCGCCGCCAGCCGCCCGACATGCAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAGGGCCTgg                                                                                                                                     <  1:201583/150‑1 (MQ=255)
           gcgcaggcgcTTTTGTGGATACGTCTTGCGTAATGCCTCGAAGCAAGGCCGTAGCAAGCCAAGGTTCTTTTCCGCCGCCAGCCGCCCGACATGCAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAGGGc                                                                                                                                         <  2:926183/141‑1 (MQ=255)
           gcgcaggcgcTTTTGTGGATACGTCTTGCGTAATGCCTCGAAGCAAGGCCGTAGCAAGCCAAGGTTCTTTTCCGCCGCCAGCCGCCCGACATGCAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAGGGc                                                                                                                                         >  1:926183/1‑141 (MQ=255)
                     ttttGTGGATACGTCTTGCGTAATGCCTCGAAGCAAGGCCGTAGCAAGCCAAGGTTCTTTTCCGCCGCCAGCCGCCCGACATGCAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGa                                                                                                                      <  2:700117/150‑1 (MQ=255)
                       ttGTGGATACGTCTTGCGTAATGCCTCGAAGCAAGGCCGCATCAAGCCAAGGTTCTTTTCCGCCGCCAGCCGCCCGACATGCAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGAtt                                                                                                                    <  1:976861/150‑1 (MQ=255)
                       ttGTGGATACGTCTTGCGTAATGCCTCGAAGAACGGCCGTAGCAAGCCAAGGTTCTTTTCCGCCGCCAGCCGCCCGACATGCAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGAtt                                                                                                                    <  2:850753/150‑1 (MQ=255)
                          tgGATACGTCTTGCGTAATGCCTCGAAGCAAGGCCGTAGCAAGCCAAGGTTCTTTTCCGCCGCCAGCCGCCCGACATGCAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGaa                                                                                                                 <  1:500598/150‑1 (MQ=255)
                             aTACGTCTTGCGTAATGCCTCGAAGCAAGGCCGTAGCAAGCCAAGGTTCTTTTCCGCCGCCAGCCGCCCGACATGCAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGAAcag                                                                                                              >  2:432488/1‑150 (MQ=255)
                                 gTCTTGCGTAATGCCTCGAAGCAAGGCCGTAGCAAGCCAAGGTTCTTTTCCGCCGCCAGCCGCCCGACATGCAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGAAcaggcag                                                                                                          >  2:645452/1‑150 (MQ=255)
                                           aTGCCTCGAAGCAAGGCCGTAGCAAGCCAAGGTTCTTTTCCGCCGCCAGCCGCCCGACATGCAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATCGAccc                                                                                                >  2:59397/1‑150 (MQ=255)
                                              ccTCGAAGCAAGGCCGTAGCAAGCCAAGGTTCTTTTCCGCCGCCAGCCGCCCGACATGCAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCACAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATCGACCCCTc                                                                                             >  2:799253/1‑150 (MQ=255)
                                               cTCGAAGCAAGGCCGTAGCAAGCCAAGGTTCTTTTCCGCCGCCAGCCGCCCGACATGCAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAGAGCCTGGTGGCGCCGGGCCGGGTTGAACAGGCAAGCATCGACCCCTCg                                                                                            >  1:198770/1‑150 (MQ=255)
                                               cTCGAAGCAAGGCCGTAGCAAGCCAAGGTTCTTTTCCGCCGCCAGCCGCCCGACATGCAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTACCAATAGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATCGACCCCTCg                                                                                            >  1:79106/1‑150 (MQ=255)
                                                tCGAAGCAAGGCCGTAGCAAGCCAAGGTTCTTTTCCGCCGCCAGCCGCCCGACATGCAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATc                                                                                                     >  1:1213227/1‑140 (MQ=255)
                                                tCGAAGCAAGGCCGTAGCAAGCCAAGGTTCTTTTCCGCCGCCAGCCGCCCGACATGCAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATc                                                                                                     <  2:1213227/140‑1 (MQ=255)
                                                tCGAAGCAAGGCCGTAGCAAGCCAAGGTTCTTTTCCGCCGCCAGCCGCCCGACATGCAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATCGACCCCTCgg                                                                                           <  1:718534/150‑1 (MQ=255)
                                                tCGAAGCAAGGCCGTAGCAAGCCAAGGTTCTTTTCCGCCGCCAGCCGCCCGACATGCAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATCGACCCCTCgg                                                                                           >  1:188422/1‑150 (MQ=255)
                                                tCGAAGCAAGGCCGTAGCAAGCCAAGGTTCTTTTCCGCCGCCAGCCGCCCGACATGCAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATCGACCCCTCgg                                                                                           >  1:15248/1‑150 (MQ=255)
                                                tCGAAGCAAGGCCGTAGCAAGCCAAGGTTCTTTTCCGCCGCCAGCCGCCCGACATGCAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATCGACCCCTCgg                                                                                           >  1:1135423/1‑150 (MQ=255)
                                                tCGAAGCAAGGCCGTAGCAAGCCAAGGTTCTTTTCCGCCGCCAGCCGCCCGACATGCAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCACAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATCGACCCCTCgg                                                                                           >  2:452987/1‑150 (MQ=255)
                                                tCGAAGCAAGGACGTAGCAAGCCAAGGTTCTTTTCCGCCGCCAGCCGACCGCGAAGTGTTACCGCAAGCTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATCGACCCCTCgg                                                                                           <  1:792364/150‑1 (MQ=255)
                                                  gAAGCAAGGCCGTAGCAAGCCAAGGTTCTTTTCCGCCGCCAGCCGCCCGACATGCAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATCGACCCCTCGGGc                                                                                         <  2:1135423/150‑1 (MQ=255)
                                                  gAAGCAAGGCCGTAGCAAGCCAAGGTTCTTTTCCGCCGCCAGCCGCCCGACATGCAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCACAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATCGACCCCTCGGGc                                                                                         <  1:452987/150‑1 (MQ=255)
                                                      cAAGGCCGTAGCAAGCCAAGGTTCTTTTCCGCCGCCAGCCGCCCGACATGCAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATCGACCCCTCGGGCCAGc                                                                                     >  1:30215/1‑150 (MQ=255)
                                                             gTAGCAAGCCAAGGTTCTTTTCCACCGCCAGCCGCCCGACATGCAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAGGACCTAGCTGCGCCGGGCCGGATTGAACAGGCAGACCTCGACCCCTCGGGACAGCAACTCCa                                                                              >  2:412968/1‑150 (MQ=255)
                                                              tAGCAAGCCCAAGCTCTTTTCCCGCTTCAGCCTCCCGACATGCAGTACGGCGATGTCCTTCGGGCCAAACAGCCAGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATCGACCCCTCGGGCCAGCAACTCCAg                                                                             <  1:347893/150‑1 (MQ=255)
                                                                    gCCAAGGTTCTTTTCCGCCGCCAGCCGCCCGACATGCAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATCGACCCCTCGGGCCAGCAACTCCAGGCGTTc                                                                       >  1:522514/1‑150 (MQ=255)
                                                                    gCCAAGGTTCTTTTCCGCCGCCAGCCGCCCGACATGCAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCACAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATCGACCCCTCGGGCCAGCAACTCCAGGCGTTc                                                                       <  1:799253/150‑1 (MQ=255)
                                                                     ccAAGGTTCTTTTCCGCCGCCAGCCGCCCGACATGCAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATCGACCCCTCGGGCCAGCAACTCCAGGCGTTCg                                                                      >  2:597315/1‑150 (MQ=255)
                                                                       aaGGTTCTTTTCCGCCGCCAGCCGCCCGACATGCAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATc                                                                                                     >  2:1028623/1‑117 (MQ=255)
                                                                       aaGGTTCTTTTCCGCCGCCAGCCGCCCGACATGCAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATc                                                                                                     <  1:1028623/117‑1 (MQ=255)
                                                                          gTTCTTTTCCGCCGCCAGCCGCCCGACATGCAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAgg                                                                                                                                           <  1:738863/76‑1 (MQ=255)
                                                                          gTTCTTTTCCGCCGCCAGCCGCCCGACATGCAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAgg                                                                                                                                           >  2:738863/1‑76 (MQ=255)
                                                                          gTTCTTTTCCGCCGCCAGCCGCCCGACATGCAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATCGACCCCTCGGGCCAGCAACTCCAGGCGTTCGAAAcc                                                                 >  2:1062326/1‑150 (MQ=255)
                                                                              ttttCCGCCGCCAGCCGCCCGACATGCAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATCGACCCCTCGGGCCAGCAACTCCAGGCGTTCGAAACCGCGa                                                             >  1:533468/1‑150 (MQ=255)
                                                                              ttatcctccGCCAGCCGCCGGACCTGCAGTAGGGCGATGTCGTCCGGGCCTAGCCCCCAGCGTTCGCGCAGGGCCTGGATGCGGCGGGCCGAACTGAGCATGCAGGCATCGACCCCTCGGGTCAGCAACTCCAGGCGTTCGAAACCGCGa                                                             <  2:198770/147‑1 (MQ=255)
                                                                               tttCCGCCGCCAGCCGCCCGACATGCCGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATCGACCCCTCGCGCCAGCAACTCCAGGCGTTCGAAACCGCGAc                                                            >  1:280213/1‑150 (MQ=255)
                                                                               tttCCGCCGCCAGCCGCCCGACATGCAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATCGACCCCTCGGGCCAGCAACTCCAGGCGTTCGAAACCGCGAc                                                            <  1:1062326/150‑1 (MQ=255)
                                                                               tttCCGCCGCCAGCCGCCCGACATGCAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAGGGCCTGGATCACTCGGCCCGGATTGAACAGGCAGGCATCGACCCCTCGGGCCAGCAACTCCAGGCGTTCGAAACCGCGAc                                                            >  1:22603/1‑150 (MQ=255)
                                                                                                 cGACATGCAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAGGGCCTGGCTGCGCCGGGAGGGCATGAACATGCAGGCATCGACCCCTCGGGCCAGCAACTCCAGGCGTTCGAAACCGCGACGCTCCAGCTCCAGGCgct                                          >  2:43992/1‑150 (MQ=255)
                                                                                                  gACATGCAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATCGACCCCTCGGGCCAGCAACTCCAg                                                                             <  2:511288/114‑1 (MQ=255)
                                                                                                  gACATGCAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATCGACCCCTCGGGCCAGCAACTCCAg                                                                             >  1:511288/1‑114 (MQ=255)
                                                                                                   aCATGCAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATCGACCCCTCGGGCCAGCAACTCCAGGCGTTCGAAACCGCGACGCTCCAGCTCCAGGCgctgg                                        <  1:59397/150‑1 (MQ=255)
                                                                                                        cAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATCGAcc                                                                                                 >  2:599205/1‑88 (MQ=255)
                                                                                                        cAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATCGAcc                                                                                                 <  1:599205/88‑1 (MQ=255)
                                                                                                        cAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATCGACCCCTCGGGCCAGCAACTCCAGGCGTTCGAAACCGCGACGCTCCAGCTACAGGCGCtggctggc                                   >  2:891090/1‑150 (MQ=255)
                                                                                                                cGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAGGGCGTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATCGACCCCTCGGGCCAGCa                                                                                    >  1:1225800/1‑93 (MQ=255)
                                                                                                                cGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATCGACCCCTCGGGCCAGCa                                                                                    <  2:1225800/93‑1 (MQ=255)
                                                                                                                  aTGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATCGACCCCTCGGGCCAGCAACTCCAGGCGTTCGAAACCGCGACGCTCCAGCTCCAGGCgctggctg                                     <  2:46985/138‑1 (MQ=255)
                                                                                                                  aTGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATCGACCCCTCGGGCCAGCAACTCCAGGCGTTCGAAACCGCGACGCTCCAGCTCCAGGCgctggctg                                     >  1:46985/1‑138 (MQ=255)
                                                                                                                      cctccGGGCCTAGCCCCCAGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATCGACCCCTCGGGCCAGCAACTCCAGGCGTTCGAAACCGCGACGCTCCAGCTCCAGGCGCTGGCTGGCACTCGGCACCAAgg                     >  1:1091438/1‑150 (MQ=255)
                                                                                                                       ctccGGGCCTAGCCCCCAGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATCGACCCCTCGGGCCAGCAACTCCAGGCGTTCGAAACCGCGAACAGCTCTCTCCGGGCGCTGGCTGGCACTCGGCACCAAGGt                    >  2:426559/1‑150 (MQ=255)
                                                                                                                        tccGGGCCTAGCCCCCAGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATCGACCCCTCGGGCCAGCAACTCCAGGCGTTCGAAACCGCGAc                                                            >  2:677073/1‑109 (MQ=255)
                                                                                                                        tccGGGCCTAGCCCCCAGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATCGACCCCTCGGGCCAGCAACTCCAGGCGTTCGAAACCGCGAc                                                            <  1:677073/109‑1 (MQ=255)
                                                                                                                                tAGCCCCCAGCTTTCGAGCAGGGCCGAGAAGCGCAGGGCGAATTTGAACAGGCAGGCATCGACCCCTCGGGCCAGCAACTCCAGGCGTTCGAAACCGCGACGCTCCAGCTCCAGGCGCTGGCTGGCACTCGGCACCAAGGTGATTGCCGt           <  2:79106/150‑1 (MQ=255)
                                                                                                                                   cccccAGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATCGACCCCTCGGGCCAGCAACTCCAGGCGTTCGAAACCGCGACGCTCCAGCTCCAGGCGCTGGCTGGCACTCGGCACCAAGGTGATTGCCGTgcg        >  2:731594/1‑150 (MQ=255)
                                                                                                                                    ccccAGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATCGACCCCTCGGGCCAGCAACTCCAGGCGTTCGAAACCGCGACGCTCCAGCTCCAGGCGCTGGCTGGCACTCGGCACCAAGGTGATTGCCGTgcgc       >  1:948031/1‑150 (MQ=255)
                                                                                                                                        aGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATCGACCCCTCGGGCCAGCAACTCCAGGCGTTCGAAACCGCGACGCTCCAGCTCCAGGCgctggctg                                     <  2:712511/116‑1 (MQ=255)
                                                                                                                                        aGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATCGACCCCTCGGGCCAGCAACTCCAGGCGTTCGAAACCGCGACGCTCCAGCTCCAGGCgctggctg                                     >  1:712511/1‑116 (MQ=255)
                                                                                                                                        aGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATCGACCCCTCGGGCCAGCAACTCCAGGCGTTCGAAACCGCGACGCTCCAGCTCCAGGCGCTGGCTGGCACTCGGCACCAAGGTGATTGCCGTGCGCCGAt   <  2:948031/150‑1 (MQ=255)
                                                                                                                                         gCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATCGACCCCTCGGGCCAGCAACTCCAGGCGTTCGAAACCGCGACgctcc                                                       >  1:522983/1‑97 (MQ=255)
                                                                                                                                         gCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATCGACCCCTCGGGCCAGCAACTCCAGGCGTTCGAAACCGCGACgctcc                                                       <  2:522983/97‑1 (MQ=255)
                                                                                                                                         gCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATCGACCCCTCGGGCCAGCAACTCCAGGCGTTCGAAACCGCGACGCTCCAGCTCCAGGCGCTGGCTGGCACTCGGCACCAAGGTGATTGCCGTGCGCCGATg  <  2:1091438/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
ACCATGATCAGGCGCAGGCGCTTTTGTGGATACGTCTTGCGTAATGCCTCGAAGCAAGGCCGTAGCAAGCCAAGGTTCTTTTCCGCCGCCAGCCGCCCGACATGCAGTACGGCGATGTCCTCCGGGCCTAGCCCCCAGCTTTCGCGCAGGGCCTGGCTGCGCCGGGCCGGATTGAACAGGCAGGCATCGACCCCTCGGGCCAGCAACTCCAGGCGTTCGAAACCGCGACGCTCCAGCTCCAGGCGCTGGCTGGCACTCGGCACCAAGGTGATTGCCGTGCGCCGATG  >  NC_002947/895684‑895970

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