Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 788,581 C→A 11.0% intergenic (+31/‑182) ettA → / → gdhA ADP/ATP ratio sensor/translation inhibitor/NADP‑dependent glutamate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947788,5810CA11.0% 176.7 / 13.3 64intergenic (+31/‑182)ettA/gdhAADP/ATP ratio sensor/translation inhibitor/NADP‑dependent glutamate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (27/30);  new base A (2/5);  total (29/35)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.29e-01

CTCGAACGTGGTGTTCTTCGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTACATTCCT  >  NC_002947/788433‑788729
                                                                                                                                                    |                                                                                                                                                    
cTCGAACGTGGTGTTCTTCGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGccc                                                                                                                                                      >  2:60199/1‑149 (MQ=255)
cTCGAACGAGGTGTTCTTCGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGccc                                                                                                                                                      >  2:315916/1‑149 (MQ=255)
       gtggtgTTCTTCGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACAGCACAAGAAACTGGCCCTGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTgg                                                                                                                                               <  1:220146/149‑1 (MQ=255)
         ggtgTTCTTCGAAGGTAACTACACCGAGTACGAAGCCGACCGCGAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGcc                                                                                                                                             >  4:231873/1‑149 (MQ=255)
                ttcGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCAtttt                                                                                                                                      >  1:148836/1‑149 (MQ=255)
                ttcGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCAtttt                                                                                                                                      >  8:212378/1‑149 (MQ=255)
                  cGAAGGTAAGTACACCGAGTACGAAGCCGACCGCAAGAAGGGTCTGGGTGATGCTGGTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGTCCCAGTAAGCCGGTTTTCTGATGCACAAAAATGGGGCCACTTGGGGCCCCAttctt                                                                                                                                     <  2:224439/148‑1 (MQ=255)
                  cGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCAttttt                                                                                                                                     <  5:108526/148‑1 (MQ=255)
                  cGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTtg                                                                                                                                    >  6:160113/1‑149 (MQ=255)
                   gAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTCTACGGCACAAGAAACTGCCCCAGTAAGCCGGTTTTCTGATGCACAACAATGGGGCCACTTGGGGCCCCATTCTtgt                                                                                                                                   <  5:172066/149‑1 (MQ=255)
                    aaGGTAACTACACCGAGTACGACGCCGCCCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAAAATGGGGCCACTTGGGGCCCCATTCTtgtg                                                                                                                                  <  6:263495/149‑1 (MQ=38)
                         aaCTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTgg                                                                                                                             >  1:68924/1‑149 (MQ=255)
                         aaCTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTgg                                                                                                                             >  3:308565/1‑149 (MQ=255)
                           cTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCTGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCg                                                                                                                           <  5:79940/149‑1 (MQ=255)
                           cTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCGCATTTTTGTGCCTGGCg                                                                                                                           >  6:77166/1‑149 (MQ=255)
                            tACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCg                                                                                                                           >  5:75103/1‑148 (MQ=255)
                                              gACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATggcgag                                                                                                        <  2:148836/149‑1 (MQ=255)
                                                  gCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATg                                                                                                             >  6:315162/1‑140 (MQ=255)
                                                  gCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATg                                                                                                             <  5:315162/140‑1 (MQ=255)
                                                   cAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGaa                                                                                                   <  8:79347/149‑1 (MQ=255)
                                                    aagaagCGTCTGGGTGATGCTGGTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAg                                                                                                  <  6:188466/149‑1 (MQ=255)
                                                      caagCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGACTACGGCACAACAAACACCCCCAGTAACCCGCTATTATGATGCACAAGAATGGGGCCCCAAGTGGCCCCAATTTTGTCCCTCGCGGACCGAAATAAATGGCGAGGCCAAgct                                                                                                >  4:260922/2‑147 (MQ=255)
                                                        agCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGgtg                                                                                               >  6:101334/1‑148 (MQ=255)
                                                            tCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAAAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTAttt                                                                                          >  6:32107/1‑149 (MQ=255)
                                                               gggTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGTGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTg                                                                                        <  3:231873/148‑1 (MQ=255)
                                                                ggTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGt                                                                                       <  4:308565/148‑1 (MQ=255)
                                                                ggTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGt                                                                                       <  2:68924/148‑1 (MQ=255)
                                                                      gctgctgcCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGaa                                                                                                   <  4:247878/130‑1 (MQ=255)
                                                                      gctgctgcCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGaa                                                                                                   >  3:247878/1‑130 (MQ=255)
                                                                      gctgctgcCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGtataagat                                                                                >  5:145140/1‑145 (MQ=255)
                                                                      gctgctgcCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGtat                                                                                     <  6:145140/144‑1 (MQ=255)
                                                                       ctgctgcCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGTCCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACtt                                                                               <  1:60199/149‑1 (MQ=255)
                                                                           tgcCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAAAATGGGGCCACTTGGGGCCCCATTCTtgt                                                                                                                                   <  2:60655/93‑1 (MQ=25)
                                                                            gcCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGgt                                                                                                <  2:143737/127‑1 (MQ=255)
                                                                              ccAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATaa                                                                         >  5:40249/1‑148 (MQ=255)
                                                                               cAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATaaaa                                                                       <  5:101334/149‑1 (MQ=255)
                                                                                aGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAAc                                                                      >  6:172599/1‑149 (MQ=255)
                                                                                 gCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGta                                                                                      >  3:254695/1‑132 (MQ=255)
                                                                                 gCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGta                                                                                      <  4:254695/132‑1 (MQ=255)
                                                                                   cggaCCGTGTACCGCACAAGAAACTGGCCCCGTACTCCGGTTTTCTGATCCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCa                                                                   <  4:185844/146‑1 (MQ=255)
                                                                                   cGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAAAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCa                                                                   <  5:32107/149‑1 (MQ=255)
                                                                                        cGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTAtttt                                                                                         >  2:217162/1‑122 (MQ=255)
                                                                                         gtgtACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCAtttt                                                                                                                                      >  5:274490/1‑76 (MQ=255)
                                                                                         gtgtACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCAtttt                                                                                                                                      <  6:274490/76‑1 (MQ=255)
                                                                                              cGGCACAAGAAACTGGCCCATTAAGCCGGTTTTCTTATGCACAATAAAGGGGCCCCAAGTGTCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGTTGTATTTTGTATACACTTACATAAAACGCACCAAATAAtt                                                         <  5:77166/148‑1 (MQ=255)
                                                                                                 cacaAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCa                                                      >  5:178316/1‑148 (MQ=255)
                                                                                                  acaAGAAACTTGCCCAGTAATCCGGTTTTCTGATGCACAAGAATTGGTCCCCAAGTGGCCACATTTTTGTGCCTGTCTGACCTATATAAATGGCGAGGCGAAGGTGTATTTTGTATACACtt                                                                               <  3:156864/122‑1 (MQ=255)
                                                                                                    aaGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCaaaa                                                   >  6:293719/1‑148 (MQ=255)
                                                                                                     aGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCaaaaa                                                  >  5:274698/1‑148 (MQ=255)
                                                                                                      gAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATggcga                                                                                                         <  5:102954/92‑1 (MQ=255)
                                                                                                      gAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATggcga                                                                                                         >  6:102954/1‑92 (MQ=255)
                                                                                                      gAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCaaaaaaa                                                <  6:75103/149‑1 (MQ=255)
                                                                                                             gCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACAtt                                         <  7:212378/149‑1 (MQ=255)
                                                                                                                 aGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGc                                     <  6:178316/149‑1 (MQ=255)
                                                                                                                  gTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGcc                                    <  6:274698/149‑1 (MQ=255)
                                                                                                                      gCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTAt                                >  1:100090/1‑149 (MQ=255)
                                                                                                                          gTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCg                            <  1:315916/149‑1 (MQ=255)
                                                                                                                               cTGATGCACAAGAATGGGGCCCCAAGGGGCCCCATTTTTGTGCCTGGCGGACCTCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGc                        <  3:260922/148‑1 (MQ=255)
                                                                                                                                  aTGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCg                       <  2:266292/146‑1 (MQ=255)
                                                                                                                                    gCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGcc                                    <  6:287855/131‑1 (MQ=255)
                                                                                                                                    gCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGcc                                    >  5:287855/1‑131 (MQ=255)
                                                                                                                                      acaAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGTGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGct                 >  5:266138/1‑148 (MQ=255)
                                                                                                                                        aaaaaTGGGGCCACTTGGGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCa                                                                   >  6:127598/4‑96 (MQ=38)
                                                                                                                                        aaaaaTGGGGCCACTTGGGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCa                                                                   <  5:127598/93‑1 (MQ=38)
                                                                                                                                          gAATGGGCCCCCAAGTGGCCCCATATTTGTGCCTGGCGGACCGCTATAAATGGCGAGGAGAAGGTGTATTTTGTATACGATTACATAAAACTCACCAAAAAATTTCAAAAAAACGACATTTCGCCCTATAACGCTGCGATTGCTTCTTg            >  6:225419/1‑149 (MQ=255)
                                                                                                                                             tGGGGCCACTTGGGGCCCCATTCTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTCGGTa         >  2:34721/1‑149 (MQ=255)
                                                                                                                                                  ccccAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTAGTCCGGTGCGATTGCTTCTTGGTACATTc    >  2:134313/1‑149 (MQ=255)
                                                                                                                                                    ccAAGTGGCCCCATCTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTACATTCCt  >  3:240843/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTCGAACGTGGTGTTCTTCGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTACATTCCT  >  NC_002947/788433‑788729

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: