Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,542,076 C→T 6.5% intergenic (+56/+47) nadE → / ← PP_4870 NH(3)‑dependent NAD(+) synthetase/azurin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,542,0760CT 6.5% 524.2 / 19.8 190intergenic (+56/+47)nadE/PP_4870NH(3)‑dependent NAD(+) synthetase/azurin
Reads supporting (aligned to +/- strand):  ref base C (78/94);  new base T (4/8);  total (83/107)
Fisher's exact test for biased strand distribution p-value = 5.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.06e-02

Only 100 of 193 total aligned reads displayed.

GCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTCGAAC  >  NC_002947/5541927‑5542222
                                                                                                                                                     |                                                                                                                                                  
 cATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCCCGAATACGAGGGTGAACCCACCCACGCa                                                                                                                                                   >  2:4362/1‑150 (MQ=255)
 cATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTCAACCGCCAATACGAGGGTGAAcccacccaccca                                                                                                                                                   >  2:100768/1‑150 (MQ=255)
      gCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGATGGTGAACCCACCGACGCaatagc                                                                                                                                              >  2:4463933/1‑145 (MQ=255)
        aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                            <  1:7097956/150‑1 (MQ=25)
        aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                            <  2:148373/150‑1 (MQ=25)
        aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                            <  2:8338593/150‑1 (MQ=255)
        aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                            <  2:3333077/150‑1 (MQ=255)
        aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                            <  1:4420434/150‑1 (MQ=255)
        aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                            <  2:5313266/150‑1 (MQ=255)
        aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                            <  2:6531424/150‑1 (MQ=255)
        aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                            <  2:7983939/150‑1 (MQ=255)
         gCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCggg                                                                                                                                           <  1:4477474/150‑1 (MQ=255)
          ccGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGt                                                                                                                                          <  1:5117890/150‑1 (MQ=25)
            ggTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGa                                                                                                                                        <  1:2374991/150‑1 (MQ=25)
              ttGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCACAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAAc                                                                                                                                      <  1:928026/150‑1 (MQ=255)
                gACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTCGCCCCCCAAAACGAGGGGGTACCCACCCCCGCAGTAGCGGGTGAAccc                                                                                                                                    >  2:2190563/1‑150 (MQ=255)
                    cggCAGCGTTCGACAACATCTTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATACGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTc                                                                                                                                <  1:7944844/149‑1 (MQ=255)
                    aGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCCGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGAACCCGCTc                                                                                                                                <  2:7777743/150‑1 (MQ=21)
                     ggCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAAAACGAGGGTGAAACCACCGACGCAATAGCGGGTGGAACCGCGcc                                                                                                                               >  2:6958904/1‑150 (MQ=255)
                            tccGAAATCATCGAGGCCACCTACCGCACGCCCCACCACACCGGGGAACTGCGGTTCGCTCCTTGATATGTGGCCTGAAGGAGGGGTTTAACCCTGTATTACTTGGGTCAACCCACAGACGCTTTCGCGGGTGAACCCGCTCCTGCattt                                                                                                                        <  2:6736596/148‑4 (MQ=255)
                                acacCAGCGTCGCCCCCTACCGAACCACCCCCCCCCCGCGCGACCGGCCGTCCGCTCCTTGATCCGTGGCCGGCGGGAGGGGGTACCCCCGGGAATCCGCGGGTGACCCCACCGACGCATCGGCGGGTGACCCCGCTCCTGCAGGGTTCg                                                                                                                    <  2:4422933/150‑1 (MQ=255)
                                aCATCATCGTCGCCACCTCCCCCAAGCCCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATACGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCg                                                                                                                    <  2:558340/150‑1 (MQ=255)
                                  atcatcGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGGGGCCTGCAGGCGCGGGTTCACCCGCGAATACGCGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGcc                                                                                                                  <  2:3434550/150‑1 (MQ=255)
                                    catcGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTCACCCGCGAATACGAGGGTGAACCCACCGACGCAGTCGCGGGTGAACCCGCGCCGGCAGGGTTCGCCGt                                                                                                                >  2:4315773/1‑150 (MQ=255)
                                      tggtcgCCACCCACCGCACGACCCCCCACACGCGCGACCTGCCGTGCGCTCCTTTATACGGGGCCGGCAGGAGCGGGTTCACCCGCGAATCCGCGGGTGACCCCACCGACGCCTTCGGGGGTGACCCCGCTCCTGCAGGGTTCGCCGTCa                                                                                                              <  2:140031/148‑1 (MQ=255)
                                      tcgtcgCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTCACCCCCCCAGACCAGGGTTAACCCACCCACCCAGTAGCGGGTGAACCCCCGACAGCAGGGTTCGCCggcc                                                                                                              >  2:3250229/1‑147 (MQ=255)
                                      tcgtcgCCACCTACCGCAAGACCCACCACACGCGCGAACTGCGGTTCGCTCCTTGATAAGTGGCCGGCAGGCGCGGTTTCACCCGCGAATCCGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCa                                                                                                              <  1:3632585/150‑1 (MQ=38)
                                      tcgtccCCACCCCCCCCACGACCCACCCCACGCGCGAACTGCGGTTCGCCCCTTTATACGTGGCCGGCCGGTGGGGTTTCACCGGCTAATACGCGGGTTAACCCCCCGACGCCTTCGCGGTTGAACCCGCTCCTGCAGGGTTCGCCGTCa                                                                                                              <  2:3915323/150‑1 (MQ=255)
                                        gtcgCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCAGCGAAGACGAGGGTGAACCCACCGACGCAATCGCGGGTGAACCCCCTCCTGCAGGGTTCGCCGGCAGc                                                                                                            >  2:92140/1‑150 (MQ=255)
                                           gCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAAAACGAGGGTGAACCCACCGACGCAATCGCGGGTGAAACCGCTCCAGCAGGGGTCGCCGTCAGccgt                                                                                                         >  2:1285104/1‑148 (MQ=255)
                                            ccACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGCGGGTGAACCCACCGCCGCCTTCGCGGGTGACCCCGCGCCTGCCGGGTTCGCCGTCAGCCttt                                                                                                        >  2:5862935/1‑150 (MQ=255)
                                             cACCTACCGCAAGACCCCGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGCGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTg                                                                                                       <  2:6387194/150‑1 (MQ=255)
                                                  aCCGCAAGACCCAGCACAAGCGCGAACTGCGGTTCGCTCCTTGATACGTGGCCTGCCGGTGGGGGTTCACCCGCTAATCCGCGGGTCAACCCACCGACGCATTCGGGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTAc                                                                                                  <  1:2800553/150‑1 (MQ=255)
                                                    cGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCAGCGCCAGCAGGGTTCGCCGGCAGCCGATGATTaagt                                                                                                >  2:4963554/1‑147 (MQ=38)
                                                    cGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTCACCCCCGAATACGAGGGTGAACCCACCCACGCAGTCGCGGGGGACCCCCCGCCGGCCGGGTTCGCCGGCACCCGTTGATTAcat                                                                                                >  2:140616/1‑148 (MQ=255)
                                                      cAACACCCACCACACGCGCTAACTGCCTTTCGCTCCTTTATAAGGGGCCTGCAGGAGCGGTTTCCCCCGCTAATACGCGGGTGAACCCCCCGACGCTTGCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGa                                                                                              <  1:5683578/150‑1 (MQ=255)
                                                          aCCCCGCACAAGCGCGCCCGGCGGTCCGCACCTTGATACGGGGCCTGCAGGAGCGGGTTCACCCGCGAATACGCGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGAt                                                                                          <  2:7522606/150‑1 (MQ=25)
                                                                cacaAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCtt                                                                                    <  2:7893357/150‑1 (MQ=255)
                                                                        cgAACTGCGGTTCGCTCCTTTATACGTGCCGTGCTGGCGCGGGTTCACCCGCTAATAGGGGGGTGAACCCACCGACGCTGTCGCGGGTGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTtca                                                                            <  2:3671456/150‑1 (MQ=255)
                                                                        cgAACTGCCGTTCGCTCCTTTATAAGTGGCCTGCAGGAGCGGGTTCACCCGCTAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTtca                                                                            <  2:1330933/150‑1 (MQ=255)
                                                                        cgAACTGCCGTTCGCTCCTTGTTAAGTGGCCTGCTGGCGCGGGTTCCCCCGCGAATACGAGGGTGAACCCACCTACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTtca                                                                            <  2:1565416/150‑1 (MQ=39)
                                                                        cgAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTtca                                                                            <  2:8692609/150‑1 (MQ=255)
                                                                                cGTTCGCTCCTTGATATGTGGCCGGCTGGAGCGGGTTCACCCGCGAAAACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCgag                                                                    <  2:5428868/150‑1 (MQ=255)
                                                                                cGTTCGCTCCTTGATACGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCgag                                                                    <  2:8606895/150‑1 (MQ=255)
                                                                                  ttCGCTCCTTGATAAGTGGCCTGCTGGCGGGGGTTCCCCGGCTAATCGGCGGGTGACCCCACCGACGCTTTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGAt                                                                  <  2:1138863/150‑1 (MQ=255)
                                                                                       tCCTTGATAAGTGGCCTGCAGGAGCGGGTTAAACCCCCAAAAACAGGGTGAAACCAACCACGCCATAGCGGGTGAACCCCCGCCTGCAGGGTTCGCCGTCAGCCGGTGATTACGTGTCGAGGACCGTGCCGTTCCTCCACGCGATGGGGc                                                             >  2:3848030/1‑150 (MQ=255)
                                                                                                   ggccggccGGCGGGGGTCCCCCCGCTAATCGGCGGGTCACCCCCCCGCCGCATGCGGGGGTGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGc                                                 <  1:449048/150‑1 (MQ=255)
                                                                                                   ggCCGGCTGGCGGGGGTCCCCCCGCCACTAGGTGGGTCACCCCACCGCCGCTTTCGGGGGTCACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGc                                                 <  1:62871/150‑1 (MQ=255)
                                                                                                     ccTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCag                                               <  2:7084500/150‑1 (MQ=255)
                                                                                                       tGCTGGCGCGGGTTCCCCCGCGATTACGCGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCagaa                                             <  1:8803503/150‑1 (MQ=255)
                                                                                                         tattcgCGGGTGAACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCagaaga                                           >  2:2247363/6‑150 (MQ=255)
                                                                                                             cgggggTCCCCCCGCAAAACCGCGGGTCACCCCCCCGCCGCATGCGCGGGTCACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCt                                       <  1:2941298/147‑1 (MQ=255)
                                                                                                               cGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGt                                     >  2:5931037/1‑150 (MQ=21)
                                                                                                                gggTTCACCCGCGAATACGCGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTa                                    <  2:3643978/150‑1 (MQ=255)
                                                                                                                gggTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTa                                    >  1:6667205/1‑150 (MQ=255)
                                                                                                                gggTTCACCCGATAAAACGTGGGCCAACCCACCGACGCATTGGCGGGTCAACCCGCTCCTGCAGGGTTCGCCGACAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTa                                    <  1:2517961/150‑1 (MQ=255)
                                                                                                                 ggTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAg                                   <  2:8817010/150‑1 (MQ=255)
                                                                                                                  gTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGt                                  >  2:2719582/1‑150 (MQ=255)
                                                                                                                  gTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGt                                  >  1:481584/1‑150 (MQ=255)
                                                                                                                   ttCACCCGCTATTACGTGGTTGAACCCACCGACGCTTTCGGGGGTGAACCCGCTCCTGCGGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTc                                 <  1:1307747/150‑1 (MQ=255)
                                                                                                                   ttCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTca                                >  1:3649508/1‑150 (MQ=255)
                                                                                                                   tacccccGCCAATACGCGGGTCACCCCACCGCCGATTTCGGGGGTTAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTc                                 <  1:464642/146‑1 (MQ=255)
                                                                                                                   ctcccccGCGAATACGCGGGTCACCCCACCGCCGCCTCCGGGGGTCACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTc                                 <  1:100768/146‑1 (MQ=255)
                                                                                                                    tCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTca                                >  1:3863070/1‑150 (MQ=255)
                                                                                                                        ccGCGAATGCGTCGGTGGATTCACCCTCGTATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                            >  1:7362138/1‑150 (MQ=37)
                                                                                                                        ccGCGAATGCGTCGGTGAGTTCACCCTCGTATTCGCGGGTGAACCCGCTCCTGCAGGGTTTGCCGTCAGCCGTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                            >  1:7260398/1‑150 (MQ=25)
                                                                                                                        ccGCGAATACGCGGGTGAACCCACCGACGCTTTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                            <  2:7403689/150‑1 (MQ=255)
                                                                                                                        ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                            >  1:6887333/1‑150 (MQ=255)
                                                                                                                        ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                            >  1:7959346/1‑150 (MQ=255)
                                                                                                                        ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                            >  1:8537599/1‑150 (MQ=255)
                                                                                                                        ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                            >  1:8800361/1‑150 (MQ=255)
                                                                                                                        ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                            >  1:2033453/1‑150 (MQ=255)
                                                                                                                        ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                            >  2:6131921/1‑150 (MQ=255)
                                                                                                                        ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                            >  2:1705786/1‑150 (MQ=255)
                                                                                                                        ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                            >  2:3587516/1‑150 (MQ=255)
                                                                                                                        ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                            >  1:3313097/1‑150 (MQ=255)
                                                                                                                        ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                            >  1:3989995/1‑150 (MQ=255)
                                                                                                                         cgcgAATCCGCGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCaccacc                           <  2:465110/150‑1 (MQ=255)
                                                                                                                         cgcgAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGc                                                 <  2:5930010/128‑1 (MQ=255)
                                                                                                                         cgcgAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCaccacc                           >  2:8519528/1‑150 (MQ=255)
                                                                                                                         cgcgAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCaccacc                           >  1:3615859/1‑150 (MQ=255)
                                                                                                                         cgcgAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCaccacc                           >  1:1480735/1‑150 (MQ=255)
                                                                                                                         cgcgAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCaccacc                           >  2:4087092/1‑150 (MQ=255)
                                                                                                                          gcgAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCg                          >  2:8631683/1‑150 (MQ=255)
                                                                                                                                ccgcgggTCAACCCCCCGACGCAGGCGGGGGTCACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTcc                    <  1:1062807/146‑1 (MQ=255)
                                                                                                                                 cGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCa                   <  2:4986494/150‑1 (MQ=255)
                                                                                                                                   aGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTGTGATTACGTGACGATGACCTTGCCTGTCATCCTCGAGATGTGGCCCGGGCAGGTGCAGACGAAGCTGTAGTCACCACCGGCTTCCAGc                 >  2:3007714/1‑150 (MQ=255)
                                                                                                                                     ggTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCtt               <  2:2280236/150‑1 (MQ=255)
                                                                                                                                      gttaaCCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGAGTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTg              <  2:6281814/147‑1 (MQ=255)
                                                                                                                                      gTGACCCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTg              <  2:1783959/150‑1 (MQ=255)
                                                                                                                                      gTGAACCCACCGTCGCTTTCGCGGGTGCCCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTg              <  1:8398156/150‑1 (MQ=255)
                                                                                                                                      gTGAACCCACCGACGCATTCGCGGGTGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTg              <  1:1314463/150‑1 (MQ=255)
                                                                                                                                      gTGAACCCACCCCCGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTg              <  1:4951029/150‑1 (MQ=255)
                                                                                                                                         aaCCCACCTCCGATTTCGGGGGTCACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGaa           <  1:6176517/150‑1 (MQ=255)
                                                                                                                                         aaCCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGaa           >  1:264623/1‑150 (MQ=255)
                                                                                                                                          aCCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAg          >  2:4572134/1‑150 (MQ=255)
                                                                                                                                             cACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAgtgt       >  1:717951/1‑150 (MQ=255)
                                                                                                                                              aCCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGCTTACTTGACGATGACCTTGCCTTTCATCCTCGAGCTGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTc      >  2:6187574/1‑150 (MQ=255)
                                                                                                                                               ccGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTCg     >  2:714569/1‑150 (MQ=255)
                                                                                                                                                  aCGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTCGAac  <  2:4877157/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTCGAAC  >  NC_002947/5541927‑5542222

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