Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 381,856 T→A 7.2% D230V (GAC→GTC)  PP_0317 ← methyl‑accepting chemotaxis transducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947381,8560TA 7.2% 205.2 / 5.1 69D230V (GAC→GTC) PP_0317methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base T (31/33);  new base A (3/2);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.42e-01

TCTTCTGCCGCCGAGGCCAGCTGGGTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGCCGATCACCAGGGCGAAAGCGCCGGCATA  >  NC_002947/381709‑381992
                                                                                                                                                   |                                                                                                                                         
tcttctGCCGCCGAGGCCAGCTGGGTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTc                                                                                                                                          >  7:91519/1‑149 (MQ=255)
tcttctGCCGCCGAGGCCAGCTGGGTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTc                                                                                                                                          >  6:358511/1‑149 (MQ=255)
  ttctGCCGCCGAGGCCAGCTGGGTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCgg                                                                                                                                        >  3:276678/1‑149 (MQ=255)
       ccgccgAGGCCAGCTGGGTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACt                                                                                                                                   <  8:91519/149‑1 (MQ=255)
       ccgccgAGGCCAGCTGGGTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACt                                                                                                                                   <  8:329912/149‑1 (MQ=255)
       ccgccgAGGCCAGCTGGGTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACt                                                                                                                                   <  8:138201/149‑1 (MQ=255)
       ccgccgAGGCCAGCTGGGTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACt                                                                                                                                   >  1:71363/1‑149 (MQ=255)
       ccgccgAGGCCAGCTGGGTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCAGGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGa                                                                                                                                     <  3:278420/147‑1 (MQ=255)
        cgccgAGGCCAGCTGGGTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACtt                                                                                                                                  <  2:71363/149‑1 (MQ=255)
         gccgAGGCCAGCTGGGTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTg                                                                                                                                 >  2:14656/1‑149 (MQ=255)
          ccgAGGCCAGCTGGGTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGc                                                                                                                                <  2:309623/149‑1 (MQ=255)
               gCCAGCTGGGTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGgc                                                                                                                           >  1:196110/1‑149 (MQ=255)
               gCCAGCTGGGTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGgc                                                                                                                           >  1:236600/1‑149 (MQ=255)
                 cAGCTGGGTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGgcgc                                                                                                                         <  2:236600/149‑1 (MQ=255)
                       ggTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCAttt                                                                                                                     <  4:333420/147‑1 (MQ=255)
                       ggTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCAttt                                                                                                                     >  3:333420/1‑147 (MQ=255)
                       ggTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTACTCGGCGCAttt                                                                                                                     <  1:177663/147‑1 (MQ=255)
                         gggCGGAGCGGGCGATCTGCTGCAAGGTATCGCGCTGCTTGTCCTGCATGGGCGGCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGt                                                                                                                 <  3:88044/148‑1 (MQ=255)
                             ggAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGgtt                                                                                                             >  4:207788/1‑149 (MQ=255)
                                   ggCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGa                                                                                                       >  3:131671/1‑149 (MQ=255)
                                            ctgcAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGa                                                                                              <  5:358511/149‑1 (MQ=255)
                                              gcAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGAc                                                                                             >  2:68597/1‑148 (MQ=255)
                                                   gTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCCTGActcaggc                                                                                       >  5:39489/1‑143 (MQ=255)
                                                        gcgcAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCggg                                                                                  >  8:283397/1‑149 (MQ=255)
                                                        gcgcAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCggg                                                                                  >  5:291364/1‑149 (MQ=255)
                                                           cAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTca                                                                               <  6:242602/149‑1 (MQ=255)
                                                            aGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTcag                                                                              <  2:196110/149‑1 (MQ=255)
                                                                  tCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGaa                                                                        <  6:27550/149‑1 (MQ=255)
                                                                  gccTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGaa                                                                        <  5:44751/148‑1 (MQ=255)
                                                                   ccTGCATGCGCGCCATGGCGCGCAACAGCCTGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATc                                                                                                     <  7:590/119‑1 (MQ=255)
                                                                   ccTGCATGCGCGCCATGGCGCGCAACAGCCAGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATc                                                                                                     >  8:590/1‑119 (MQ=255)
                                                                    cTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGa                                                                                           >  8:84229/1‑128 (MQ=255)
                                                                    cTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAg                                                                       <  5:124174/148‑1 (MQ=255)
                                                                       gttGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGca                                                                    >  8:28992/3‑148 (MQ=255)
                                                                         tGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCAt                                                                                                                       >  5:40388/1‑95 (MQ=255)
                                                                         tGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAg                                                                 <  2:115874/149‑1 (MQ=255)
                                                                         tGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAg                                                                 <  1:14656/149‑1 (MQ=255)
                                                                         tGCGCGCCATGGCGCGCAACAGCAGGGCGGCTTCGTCGGTGCCCTCTGCGCCGAAGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTAGAGTGCTTCTTCGATCGGCTTGACGAAACTGCGCGTCAGAAGGAATGCACAg                                                                 >  8:305147/1‑149 (MQ=255)
                                                                            gcgcCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGc                                                               >  2:181321/1‑148 (MQ=255)
                                                                              gcCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCg                                                              >  7:221254/1‑147 (MQ=255)
                                                                               ccATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAg                                                           <  2:32818/149‑1 (MQ=255)
                                                                                  tGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTc                                                        <  6:291364/149‑1 (MQ=255)
                                                                                   ggCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCAttt                                                                                                                     >  5:362121/1‑87 (MQ=255)
                                                                                   ggCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCAttt                                                                                                                     <  6:362121/87‑1 (MQ=255)
                                                                                    gcgcgcAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCa                                                       <  4:276678/148‑1 (MQ=255)
                                                                                        gcAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACg                                                 >  6:124540/1‑149 (MQ=255)
                                                                                          aaCAGCCGGGCGGCTTCGTCGGTGCCATCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTTAGAGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCATAACGGt                                               >  8:153961/1‑149 (MQ=255)
                                                                                           tgaGCCGGGCGGCTTGGTCGGTGGCCTGTGCGCGGATGTCGTGTGTTAGGTCGCCGTCGGCGACTTGCTCGTCGCTTTTGAGTGCTTCTTGGATCGGGTTGACGATACTTCGGGTCAGCAGGAAGGCATTGGCGAAGGTCAGTACGTTg                                              <  4:131671/147‑1 (MQ=255)
                                                                                           aCAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTg                                              >  1:264440/1‑149 (MQ=255)
                                                                                             aGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTggc                                            <  4:214851/149‑1 (MQ=255)
                                                                                              gCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTggcg                                           <  6:83241/149‑1 (MQ=255)
                                                                                                   gcggcTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAAc                                      <  8:221254/149‑1 (MQ=255)
                                                                                                          cgtcgGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAg                                                           <  4:60155/122‑1 (MQ=255)
                                                                                                          cgtcgGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAgg                               >  2:146588/1‑149 (MQ=255)
                                                                                                                     tctGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGCCGATCACCAg                    >  2:228680/1‑149 (MQ=255)
                                                                                                                              aTGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGCCGATCACCAGGGCGAAAgc           >  7:269272/1‑149 (MQ=255)
                                                                                                                                gACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGCCGATCACCAGGGCGAAAgcgc         >  8:89954/1‑149 (MQ=255)
                                                                                                                                 aCGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGCCGATCACCAGGGCGAAAGCGcc        <  1:181321/149‑1 (MQ=255)
                                                                                                                                  cGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGCCGATCACCAGGGCGAAAGCGcc        <  3:207788/148‑1 (MQ=255)
                                                                                                                                   gTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACCGGTGGCGGCAACCAGCAGGCCGATCACCAGGGCGAAAGCGCCg       <  1:68597/149‑1 (MQ=255)
                                                                                                                                       gTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTg                                              <  1:59009/105‑1 (MQ=255)
                                                                                                                                       gTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTg                                              >  6:32890/1‑105 (MQ=255)
                                                                                                                                       gTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTg                                              >  2:59009/1‑105 (MQ=255)
                                                                                                                                       gTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTg                                              <  5:32890/105‑1 (MQ=255)
                                                                                                                                       agcaagTCGCCGACGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCTACCAGCAGGCCGATCACCAGGGCGAAAGCGCCGGCATa  >  2:84890/6‑149 (MQ=255)
                                                                                                                                        gcaagTCGCCGACGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGCCGATCACCAGGGCGAAAGCGCCGGCAt   <  6:55071/143‑1 (MQ=255)
                                                                                                                                        gcaagTCGCCGACGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGCCGATCACCAGGGCGAAAGCGCCGGCAt   >  5:55071/5‑147 (MQ=255)
                                                                                                                                          aagTCGCCGACGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCa                                                       <  1:321487/92‑1 (MQ=255)
                                                                                                                                          aagTCGCCGACGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCa                                                       >  2:321487/3‑94 (MQ=255)
                                                                                                                                                   |                                                                                                                                         
TCTTCTGCCGCCGAGGCCAGCTGGGTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGCCGATCACCAGGGCGAAAGCGCCGGCATA  >  NC_002947/381709‑381992

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