Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 788,584 C→A 7.6% 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,5840CA 7.6% 1118.1 / 50.7 337intergenic (+31/‑182)ettA/gdhAADP/ATP ratio sensor/translation inhibitor/NADP‑dependent glutamate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (166/145);  new base A (15/11);  total (181/156)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 344 total aligned reads displayed.

ACTCGAACGTGGTGTTCTTCGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTACATTCCTG  >  NC_002947/788435‑788733
                                                                                                                                                     |                                                                                                                                                     
  tCGAACGTGGTGTTCTTCGAAGGTAACTACAACGAGTACGAAGACGACCGCAAGAAGCGTCTGGGTGATGCTGATGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGTGGCCCCa                                                                                                                                                     >  1:3719088/1‑150 (MQ=255)
        gtggtgTTCTTCGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGc                                                                                                                                               <  2:3259712/150‑1 (MQ=255)
           gtgtTCTTCGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGcccc                                                                                                                                            >  2:6295258/1‑150 (MQ=255)
            tgtTCTTCGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCa                                                                                                                                           >  1:1110288/1‑150 (MQ=255)
             gtTCTTCGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCAt                                                                                                                                          >  1:2318852/1‑150 (MQ=255)
                 ttcGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCTAGTGGCCCCAttttt                                                                                                                                      <  1:3895016/150‑1 (MQ=255)
                 ttcGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCAttttt                                                                                                                                      <  2:5784967/150‑1 (MQ=255)
                 ttcGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCAttttt                                                                                                                                      >  2:5608119/1‑150 (MQ=255)
                 ttcGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCAttttt                                                                                                                                      >  1:5715777/1‑150 (MQ=255)
                  tcGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTtg                                                                                                                                     >  2:6568/1‑150 (MQ=255)
                   cGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGTGGCCCCAAGTGGCCCCATTTTtgt                                                                                                                                    <  1:5309061/150‑1 (MQ=255)
                   cGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTtgt                                                                                                                                    <  1:5288993/150‑1 (MQ=255)
                     aaGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAACCTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGc                                                                                                                                  <  2:4884465/150‑1 (MQ=255)
                      aGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGcc                                                                                                                                 >  1:5986465/1‑150 (MQ=255)
                        gTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTg                                                                                                                               >  1:3309551/1‑150 (MQ=255)
                        gTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCACTTGGGGCCCCATTCTTGTGcatc                                                                                                                               <  2:2252927/150‑4 (MQ=39)
                          aaCTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGc                                                                                                                             >  2:371468/1‑150 (MQ=255)
                            cTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCgg                                                                                                                           >  1:4159363/1‑150 (MQ=255)
                             tACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGa                                                                                                                          >  2:732644/1‑150 (MQ=255)
                              acacCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGAc                                                                                                                         <  2:5333262/150‑1 (MQ=255)
                               cacCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGAcc                                                                                                                        >  2:3066446/1‑150 (MQ=255)
                                  cGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCt                                                                                                                     >  2:4264689/1‑150 (MQ=255)
                                   gAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCta                                                                                                                    <  2:6835086/150‑1 (MQ=255)
                                    aGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCtat                                                                                                                   <  1:814112/150‑1 (MQ=255)
                                       aCGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATaaa                                                                                                                <  2:4044604/150‑1 (MQ=255)
                                        cGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAAt                                                                                                               <  1:749691/150‑1 (MQ=255)
                                         gAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATg                                                                                                              <  1:6190915/150‑1 (MQ=255)
                                         gAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATg                                                                                                              >  1:4612647/1‑150 (MQ=255)
                                          aaGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATgg                                                                                                             <  2:237984/150‑1 (MQ=255)
                                            gCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATggcg                                                                                                           <  1:1825030/150‑1 (MQ=255)
                                             ccGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATggcga                                                                                                          <  2:2236238/150‑1 (MQ=255)
                                               gACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATggcgagg                                                                                                        <  2:801168/150‑1 (MQ=255)
                                                   gCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGaa                                                                                                    <  1:1190115/150‑1 (MQ=255)
                                                    cAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAg                                                                                                   <  2:1751792/150‑1 (MQ=255)
                                                     aagaagCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAgg                                                                                                  >  2:4455482/1‑150 (MQ=255)
                                                     aagaagCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAgg                                                                                                  >  2:3465879/1‑150 (MQ=255)
                                                         agCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTa                                                                                              <  1:5060997/150‑1 (MQ=255)
                                                         agCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTa                                                                                              <  1:2023584/150‑1 (MQ=255)
                                                           cGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTAtt                                                                                            <  1:1340545/150‑1 (MQ=255)
                                                             tCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCGCAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTAtttt                                                                                          <  2:2052269/150‑1 (MQ=255)
                                                             tCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTAtttt                                                                                          <  1:2628458/150‑1 (MQ=255)
                                                             tCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTAtttt                                                                                          <  2:2719690/150‑1 (MQ=255)
                                                               tGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGt                                                                                        >  2:3590455/1‑150 (MQ=255)
                                                                gggTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGGGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGta                                                                                       <  2:61192/150‑1 (MQ=255)
                                                                             gcCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTGTGTATACACTTACATaa                                                                          >  2:4990166/1‑150 (MQ=255)
                                                                              cccAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATaaa                                                                         >  1:5880628/1‑150 (MQ=255)
                                                                               ccAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATaaaa                                                                        >  2:6163319/1‑150 (MQ=255)
                                                                                cAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGtat                                                                                      <  2:4747969/135‑1 (MQ=255)
                                                                                cAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGtat                                                                                      >  1:4747969/1‑135 (MQ=255)
                                                                                   ccGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCa                                                                    <  1:2689239/150‑1 (MQ=255)
                                                                                   ccGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCa                                                                    <  2:1561491/150‑1 (MQ=255)
                                                                                     gCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCAcc                                                                  >  2:3043518/1‑150 (MQ=255)
                                                                                         cGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAaat                                                              >  1:6378358/1‑150 (MQ=255)
                                                                                           tgtACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTTGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAaataa                                                            <  2:6715935/150‑1 (MQ=255)
                                                                                           tgtACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAaataa                                                            <  1:1010992/150‑1 (MQ=255)
                                                                                               cGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCTAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTc                                                        >  1:2957242/1‑150 (MQ=255)
                                                                                                 gCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCaa                                                      >  1:2833859/1‑150 (MQ=255)
                                                                                                   tcagaaaaCTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCaaaa                                                    >  2:5166792/6‑150 (MQ=255)
                                                                                                   acaAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCaaaa                                                    <  1:3465879/150‑1 (MQ=255)
                                                                                                     aaGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCaaaaaa                                                  >  2:5397262/1‑150 (MQ=255)
                                                                                                         aaCTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCAATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGa                                              <  1:1704261/150‑1 (MQ=255)
                                                                                                          aCTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCATGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGACGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGAc                                             <  1:6096378/150‑1 (MQ=255)
                                                                                                               cccAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTc                                        >  1:1518830/1‑150 (MQ=255)
                                                                                                                ccAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCg                                       >  2:245207/1‑150 (MQ=255)
                                                                                                                ccAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCg                                       >  2:4661883/1‑150 (MQ=255)
                                                                                                                  aGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCACGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGcc                                     <  1:3935732/150‑1 (MQ=255)
                                                                                                                   gTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGccc                                    <  1:4451761/150‑1 (MQ=255)
                                                                                                                        ccGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTc                               >  1:1683462/1‑150 (MQ=255)
                                                                                                                        ccGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTc                               >  1:543482/1‑150 (MQ=255)
                                                                                                                          ggTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTCACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCg                             <  1:450179/150‑1 (MQ=255)
                                                                                                                                 tGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGAt                      >  1:3516754/1‑150 (MQ=255)
                                                                                                                                 tGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGAt                      >  1:3006471/1‑150 (MQ=255)
                                                                                                                                 tGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGACCTATTCAGGTGCGAt                      >  2:4346756/1‑150 (MQ=255)
                                                                                                                                  gATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTTTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGAtt                     <  1:4455519/150‑1 (MQ=255)
                                                                                                                                  gATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGAtt                     <  2:1956251/150‑1 (MQ=255)
                                                                                                                                   aTGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCGGCGGACCGCAAAAAATGGCGCGGCCTAGGTGTAATTTTTATACCCTTACATAAAAAGCACCAAATAATTTGAAAAAAAAGGCATTTAGCCCTATTGCGGTGCGATTGc                   >  1:2392056/1‑150 (MQ=255)
                                                                                                                                    ggCACAAAAATGGGGCCACTTGGGGCCCCATTCTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTGACATAAAACGCACCAGATAATTTCAAAAAAACGACAGTTCGCCCTATTGCGGGGCGATTGc                   >  1:4942839/2‑150 (MQ=37)
                                                                                                                                       acaAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGcttc                >  1:4827483/1‑150 (MQ=255)
                                                                                                                                       acaAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGcttc                >  1:5125857/1‑150 (MQ=255)
                                                                                                                                       acaAAAATGGGGCCACTTGGGGCCCCATTCTTGTGCATGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGcttc                >  2:2887793/1‑150 (MQ=38)
                                                                                                                                       acaAAAATGGGGCCACTTGGGGCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACATACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGcttct               >  1:2239614/1‑150 (MQ=255)
                                                                                                                                        caAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGcttct               <  1:3066446/150‑1 (MQ=255)
                                                                                                                                         aaaaaTGGGGCCACTTGGGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTc                                                        <  2:5091242/105‑1 (MQ=38)
                                                                                                                                         aaaaaTGGGGCCACTTGGGGCCCCATTCTTGTGCATGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGcttctt              >  2:3369807/4‑150 (MQ=39)
                                                                                                                                         aaGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGcttctt              >  2:1671795/1‑150 (MQ=255)
                                                                                                                                         aaGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAAGGGTGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGcttctt              <  1:6848170/150‑1 (MQ=255)
                                                                                                                                          aGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTTTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTg             <  1:4058242/150‑1 (MQ=255)
                                                                                                                                          aGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTg             <  2:4933405/150‑1 (MQ=255)
                                                                                                                                          aGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTg             <  2:1091764/150‑1 (MQ=255)
                                                                                                                                          aGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTg             <  2:4194648/150‑1 (MQ=255)
                                                                                                                                           gAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTgg            >  1:2509892/1‑150 (MQ=255)
                                                                                                                                           aaaTGGGGCCACTTGGGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTgg            >  2:2540409/2‑150 (MQ=255)
                                                                                                                                            aaTGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGt           >  1:6593450/1‑150 (MQ=255)
                                                                                                                                            aaTGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGt           <  2:782913/150‑1 (MQ=255)
                                                                                                                                             aTGGGGCCCCAAGTGGACCCATTTTTGTGACTGGCGGACAGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTAAAAAAAACGACATTTCGACCTATTCCGGTGCGATTGCTTCTTGGTa          >  1:5531453/1‑150 (MQ=255)
                                                                                                                                                gggCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTACAt       >  2:1190621/1‑150 (MQ=255)
                                                                                                                                                  gCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTACATTc     >  2:1371356/1‑150 (MQ=255)
                                                                                                                                                   ccccAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTACATTcc    <  2:4524262/150‑1 (MQ=255)
                                                                                                                                                    cccAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTACATTCCt   <  1:915980/150‑1 (MQ=255)
                                                                                                                                                     ccAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTACATTCCTg  >  1:1990360/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACTCGAACGTGGTGTTCTTCGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTACATTCCTG  >  NC_002947/788435‑788733

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: