Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 788,586 A→T 6.8% intergenic (+33/‑180) 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,5860AT 6.8% 845.2 / 32.2 264intergenic (+33/‑180)ettA/gdhAADP/ATP ratio sensor/translation inhibitor/NADP‑dependent glutamate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (144/102);  new base T (7/11);  total (151/113)
Fisher's exact test for biased strand distribution p-value = 1.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 271 total aligned reads displayed.

TCGAACGTGGTGTTCTTCGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTACATTCCTGAA  >  NC_002947/788437‑788735
                                                                                                                                                     |                                                                                                                                                     
  gAACGTGGTGTTCTTCGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAg                                                                                                                                                     >  2:4810264/1‑150 (MQ=255)
     cGTGGTGTTCTTCGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTgg                                                                                                                                                  <  1:1419505/150‑1 (MQ=255)
      gtggtgTTCTTCGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGc                                                                                                                                                 <  1:5514469/150‑1 (MQ=255)
          tgtTCTTCGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCa                                                                                                                                             <  2:2512689/150‑1 (MQ=255)
           gtTCTTCGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCAt                                                                                                                                            >  1:4319646/1‑150 (MQ=255)
           gtTCTTCGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCACTTGGGGCCCCAt                                                                                                                                            <  2:2982495/150‑1 (MQ=255)
            ttcttcGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCAtt                                                                                                                                           >  1:41270/1‑150 (MQ=255)
             tcttcGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTTGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCAttt                                                                                                                                          >  1:2642697/1‑150 (MQ=255)
             tcttcGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCACTTGGGGCCCCAttc                                                                                                                                          <  1:782798/150‑2 (MQ=255)
               ttcGAAGGTAACTACACCGAGTATGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCAttttt                                                                                                                                        >  1:1342773/1‑150 (MQ=255)
               ttcGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCAttttt                                                                                                                                        <  2:1825253/150‑1 (MQ=255)
               ttcGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCAttttt                                                                                                                                        <  1:4868605/150‑1 (MQ=255)
               ttcGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCAttttt                                                                                                                                        >  2:2903057/1‑150 (MQ=255)
               ttcGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAAAATGGGGCCACTTGGGGCCCCAttctt                                                                                                                                        <  2:1899218/150‑1 (MQ=255)
                tcGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTtg                                                                                                                                       <  1:2762689/150‑1 (MQ=255)
                 cGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGATACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTtgt                                                                                                                                      <  2:2804846/150‑1 (MQ=255)
                 cGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTtgt                                                                                                                                      <  2:1115117/150‑1 (MQ=255)
                 cGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTtgt                                                                                                                                      <  1:3082198/150‑1 (MQ=255)
                  gAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTtgtg                                                                                                                                     >  1:3275797/1‑150 (MQ=255)
                  gAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGACTGGGAGATGCTGCAGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGAAAGCCGGTTTTCTGATGCACAAGAATGTGGCCCCAAGTGGCCCCATTTTtgtg                                                                                                                                     >  1:459684/1‑150 (MQ=255)
                         aCTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCg                                                                                                                              <  2:1451518/150‑1 (MQ=255)
                          cTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCgg                                                                                                                             >  2:4061850/1‑150 (MQ=255)
                               ccGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGc                                                                                                                        <  1:4972410/150‑1 (MQ=255)
                                      cGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAAt                                                                                                                 >  1:5206771/1‑150 (MQ=255)
                                      cGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAAt                                                                                                                 >  1:3970252/1‑150 (MQ=255)
                                       gAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATg                                                                                                                >  2:5015892/1‑150 (MQ=255)
                                        aaGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATgg                                                                                                               <  1:449488/150‑1 (MQ=255)
                                           ccGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATggcga                                                                                                            >  2:5317234/1‑150 (MQ=255)
                                           ccGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATggcga                                                                                                            <  2:2046485/150‑1 (MQ=255)
                                             gACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATggcgagg                                                                                                          <  1:1528327/150‑1 (MQ=255)
                                             gACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATggcgagg                                                                                                          >  2:621665/1‑150 (MQ=255)
                                                cGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATggcgaggcga                                                                                                       >  2:544673/1‑150 (MQ=255)
                                                 gCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGaa                                                                                                      >  2:3556661/1‑150 (MQ=255)
                                                 gCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGaa                                                                                                      <  1:3788643/150‑1 (MQ=255)
                                                  cAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAg                                                                                                     <  1:1080251/150‑1 (MQ=255)
                                                  cAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAg                                                                                                     >  2:5279293/1‑150 (MQ=255)
                                                   aagaagCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAgg                                                                                                    >  2:4457002/1‑150 (MQ=255)
                                                    agaagCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGgt                                                                                                   <  2:3210634/150‑1 (MQ=255)
                                                          gTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTAttt                                                                                             >  1:1016389/1‑150 (MQ=255)
                                                            cTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTg                                                                                           <  2:172375/150‑1 (MQ=255)
                                                             tGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGt                                                                                          >  1:3029201/1‑150 (MQ=255)
                                                              gggTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGta                                                                                         <  1:486118/150‑1 (MQ=255)
                                                              gggTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGta                                                                                         >  1:424958/1‑150 (MQ=255)
                                                              gggTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGta                                                                                         <  1:1560317/150‑1 (MQ=255)
                                                               ggTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGtat                                                                                        >  1:2943545/1‑150 (MQ=255)
                                                                   aTGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATacac                                                                                    <  2:2177740/150‑1 (MQ=255)
                                                                   aTGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATacac                                                                                    <  2:3124401/150‑1 (MQ=255)
                                                                        gctgcCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACa                                                                               >  2:3828109/1‑150 (MQ=255)
                                                                            cccAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATaaa                                                                           >  1:815235/1‑150 (MQ=255)
                                                                              cAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAAc                                                                         >  2:4226108/1‑150 (MQ=255)
                                                                                 ccGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCa                                                                      <  2:1541362/150‑1 (MQ=255)
                                                                                    cACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCa                                                                   >  2:1038148/1‑150 (MQ=255)
                                                                                     aCCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCaa                                                                  <  2:3666935/150‑1 (MQ=255)
                                                                                         tgtACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCTAGGCGAAGGTGTATTTTGTATACACTTACATAAAAAGCACCAaataa                                                              >  2:2370563/1‑150 (MQ=255)
                                                                                            aCGGCACAAGCTACTGGCCCATGAAGCAGGTTTTCTGATTCACATGAATGTGGCCCCAACTGGCCCAATTTTTGTGCCTGGCTGTCCGCTATAAATGGAGAGGCGTAGGTGTATTTTGTATAAACTTACATAAAACGCACCATATAAttt                                                           <  2:362415/150‑1 (MQ=255)
                                                                                             cGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTc                                                          <  1:5187463/150‑1 (MQ=255)
                                                                                               gCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCaa                                                        <  1:2462459/150‑1 (MQ=255)
                                                                                                cacaAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCaaa                                                       <  1:5062719/150‑1 (MQ=255)
                                                                                                cacaAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCaaa                                                       <  1:3690793/150‑1 (MQ=255)
                                                                                                 acaAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATTGCGAGGCGAAGGTGTATTTTTTATACACTTACATAAAACGCACCAAATAATTTAaaaa                                                      >  1:2426477/1‑150 (MQ=255)
                                                                                                   aaGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCaaaaaa                                                    >  1:4681013/1‑150 (MQ=255)
                                                                                                    aGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGACTGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCaaaataa                                                   <  2:4096958/150‑1 (MQ=255)
                                                                                                     gAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAAc                                                  <  1:58356/150‑1 (MQ=255)
                                                                                                            gCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACAttt                                           >  2:4864924/1‑150 (MQ=255)
                                                                                                              ccAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCg                                         >  1:5269272/1‑150 (MQ=255)
                                                                                                               cAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGc                                        >  2:1990640/1‑150 (MQ=255)
                                                                                                                aGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGcc                                       <  2:4626472/150‑1 (MQ=255)
                                                                                                                aGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGcc                                       >  1:2131775/1‑150 (MQ=255)
                                                                                                                aGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGcc                                       >  1:3162146/1‑150 (MQ=255)
                                                                                                                aGTAAGCCGGTTTTCTGATGCACAAGAAGGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGcc                                       <  2:4359196/150‑1 (MQ=255)
                                                                                                                    aGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTAt                                   >  2:2768373/1‑150 (MQ=255)
                                                                                                                          ttttCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGt                             <  1:1811943/150‑1 (MQ=255)
                                                                                                                               tGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGAt                        >  1:3519685/1‑150 (MQ=255)
                                                                                                                               tGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGAt                        <  1:4456789/150‑1 (MQ=255)
                                                                                                                                 aTGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACTCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCAGGTGCGATTg                      >  2:1270646/1‑150 (MQ=255)
                                                                                                                                 aTGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTg                      <  1:1990640/150‑1 (MQ=255)
                                                                                                                                 aTGCACAAGAATGGGGCCCCAAGGGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGATGGTTTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTg                      <  1:511958/150‑1 (MQ=255)
                                                                                                                                     acaAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGcttc                  <  2:5063388/150‑1 (MQ=255)
                                                                                                                                     acaAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGcttc                  >  1:1882568/1‑150 (MQ=255)
                                                                                                                                     acaAAAATGGGGCCACTTGGGGCCCCATTCTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGcttc                  >  1:887830/1‑150 (MQ=39)
                                                                                                                                      caAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGcttct                 >  2:1943760/1‑150 (MQ=255)
                                                                                                                                       aaGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGcttctt                >  1:2854273/1‑150 (MQ=255)
                                                                                                                                       aaGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGcttctt                >  2:2417112/1‑150 (MQ=255)
                                                                                                                                       aaGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGcttctt                >  1:4912744/1‑150 (MQ=255)
                                                                                                                                       aaGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGcttctt                >  1:5437199/1‑150 (MQ=255)
                                                                                                                                        ataaTGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACAGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACATACATAAAACGCACCAAATAATTTCAAAAAAAAGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTg               >  2:4446209/3‑150 (MQ=255)
                                                                                                                                        aaaaTGGGGCCACTTGGGGCCCCATTCTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTg               >  2:3291577/3‑150 (MQ=255)
                                                                                                                                        aGAATGGGGCCCCACGTGGCCCCATTTTTGTGCCTGGCGGACCGCTCTAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTg               <  1:1192788/150‑1 (MQ=255)
                                                                                                                                         gAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTgg              >  2:3384755/1‑150 (MQ=255)
                                                                                                                                         aaaTGGGGCCACTTGGGGCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGt             >  1:2813877/2‑150 (MQ=255)
                                                                                                                                             ggggCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTACa          >  2:3701421/1‑150 (MQ=255)
                                                                                                                                               ggCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTACAtt        >  1:225289/1‑150 (MQ=255)
                                                                                                                                               ggCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTACAtt        >  2:807507/1‑150 (MQ=255)
                                                                                                                                                gCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTACATTc       <  2:1433182/150‑1 (MQ=255)
                                                                                                                                                gCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTACATTc       >  2:691047/1‑150 (MQ=255)
                                                                                                                                                 ccccAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAAGGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTACATTcc      <  1:1372217/150‑1 (MQ=255)
                                                                                                                                                  cccAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTACATTCCt     >  2:3230019/1‑150 (MQ=255)
                                                                                                                                                     aaGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTa            >  2:1826734/1‑140 (MQ=255)
                                                                                                                                                     aaGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTACATTCCTGaa  <  1:3493847/150‑1 (MQ=255)
                                                                                                                                                     aaGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTACATTCCTGaa  >  2:59875/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCGAACGTGGTGTTCTTCGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTACATTCCTGAA  >  NC_002947/788437‑788735

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: