Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,065,193 C→G 12.4% intergenic (+85/‑122) PP_1838 → / → PP_1839 membrane protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,065,1930CG12.4% 254.6 / 23.3 95intergenic (+85/‑122)PP_1838/PP_1839membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (46/37);  new base G (7/5);  total (53/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.08e-01

CCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGC‑TTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGTGGGAAGTCATTTTTGC  >  NC_002947/2065045‑2065340
                                                                                                                                                     |                                                                                                                                                    
ccAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGC‑TTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGc                                                                                                                                                      <  1:3446003/149‑1 (MQ=255)
ccAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGC‑TTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGc                                                                                                                                                      >  2:3446003/1‑149 (MQ=255)
  ggACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGG‑GTGCCCGCGCCCCCAGGCAGTGGGCACGGTTCGGGCCATTTGCAGGGCCGGGGCggg                                                                                                                                                   >  1:1247884/2‑148 (MQ=255)
  aGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGG‑GTGGCCGCGCCCCCAGGCCATGGACAGGGTTGGGGCCGTGTGCCGGGCCGGGGCgcg                                                                                                                                                   >  1:3458291/1‑148 (MQ=255)
   gACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGC‑TTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGc                                                                                                                                                  <  1:3469921/150‑1 (MQ=18)
   gACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGCGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGC‑TTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGc                                                                                                                                                  <  2:4729617/150‑1 (MQ=18)
     ccGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGC‑TTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCgg                                                                                                                                                <  2:3794759/150‑1 (MQ=18)
     ccGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGC‑TTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCgg                                                                                                                                                <  2:4478153/150‑1 (MQ=18)
     ccGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGC‑TTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCgg                                                                                                                                                <  2:4174862/150‑1 (MQ=255)
     ccGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGC‑TTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCgg                                                                                                                                                <  1:1034005/150‑1 (MQ=255)
         gACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGC‑TTGCCCGCTCCCACAGGCATTGCACAGGGGTCGGACCGTGTGCAGGGACGGGGCGAGCGGGccc                                                                                                                                            >  1:3104254/1‑148 (MQ=39)
             aaGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGC‑GTTCCCGCTCCCACAGGCCTTGGACAGGGTTCGGACCTTGTGCAGGTGCAGTGCGAGCGGGCCAGccc                                                                                                                                        >  2:1427903/1‑150 (MQ=255)
               gagCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGA‑TTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCcgc                                                                                                                                      <  1:2677264/150‑1 (MQ=255)
                agCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCGTTGCCC‑CGCCCACAGGCAGTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGGGGGCGGGCCAGCCcgcg                                                                                                                                     >  2:941383/1‑150 (MQ=255)
                     gctgGATAACACCGCCGGCGTGGTGCCCGGCACACCTAATTAAACCGGGGCCGCGTTGCGGCCCATTCGGGGGC‑TTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCTGGTCGTGTGCGAGCGGGCAAGCCCGCGAATgg                                                                                                                                <  2:1448560/150‑1 (MQ=38)
                     gctgGATAACACCGCCGGCGTGGTGCAGGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGC‑TTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACGTTGTGCAGGTCCTGTGCGAGGGGGCAAGCCCGCGAATgg                                                                                                                                <  2:1448591/150‑1 (MQ=255)
                       tgGATAACACCGCCGGCGTGGTGCATGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTAGCGGGG‑TTTGCCGCTCCCAAAGGGATTTGACATGGTTTGGTACTTTTTGAGGTGCTGTGGGAGCGGGGAATCCCGCGTATGGGc                                                                                                                              >  2:1510131/1‑150 (MQ=255)
                               acCGCCGGCGTGGTGCACGGCACACCCAAGTAAACCGGGGCCGCGTTGCGGCCCATACGCGGGC‑TTTACCCCACCCCAAGCCCAAGCACAGCGGTCGGACCATGTGCAGGTACTGAGCGGGCGGGCAAACCCCAGCCTTGGCCCCGcaca                                                                                                                      >  2:2784536/1‑150 (MQ=255)
                                           gtgCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGC‑TTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGGCCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGAc                                                                                                         >  2:3986007/1‑150 (MQ=255)
                                                     cacCGAAGTAAACGGGGGCCGCGTTGCGGCCCATTCGCGGGC‑TTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCTGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTAt                                                                                               <  2:937462/150‑1 (MQ=255)
                                                             tAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGC‑TTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGccc                                                                                       >  2:5433581/1‑150 (MQ=255)
                                                                 ccGGGGCCGCGTTGCGGCCCATTCGCGGGC‑TTTCCCGCTCCCACACGCAGTGCACAGGGTTCGGACCTTGTGCAGGTCCGGTGCGAGCGTGCAAGCCCGCGACTGGGCCCCCCACAATTTT‑CCCGGGCTGAACAATATACCCGccccccg                                                                                   >  2:3743792/1‑147 (MQ=255)
                                                                 ccGGGGCCGCGTTGCGGCCCATTCGCGGGC‑TTGCCCCCGCCCCCAGGCCGTGCGCCGGGGTCGGGCCGTGTGCAGGGCCGGGGCGCGCGGGCCCGCCCGCGCCTGGGCCCCGCACCAGTTT‑ACCCGCCGGAACAGTATAccccccccccg                                                                                   >  2:3411198/1‑147 (MQ=255)
                                                                     ggCCGCGTTGCGGCCCATTCGCGGGC‑TTGCCCGCTCCCACCGGCATTGCACAGGGTTCGGACCTTGTGCTGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCa                                                                               <  2:3674742/150‑1 (MQ=255)
                                                                     gcccggtttGCCGCCCTTGCGGGGGC‑TGCCCCCCACCCCCCGGCATGGCACGGGGTTCGGACCTGGCGCCGGTCGGGGGCCCGGGGCCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCa                                                                               <  1:5154777/143‑1 (MQ=255)
                                                                      gCCGCGTTGCGGCCCATTCGCGGGC‑TTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAg                                                                              <  2:3949793/150‑1 (MQ=255)
                                                                             tggggCCCATGCGCGGGG‑TGCCCCGCCCCCACCGTATTTGCACTGGGTTCGGACGTTGTGCTGGTCGTGGCCGCGCGGCCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCa                                                                       <  1:1427903/147‑1 (MQ=255)
                                                                              ggggCCAATTCGGGGGG‑TGGCCCGCTCCCCCTGGCTTTGCACTGGGTTCGGACCTTGGGCTGGTCGTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAg                                                                      <  1:383873/148‑1 (MQ=255)
                                                                                 gCCCTTCCGCGGGC‑TGGCCCGCCCCCACCGCCTTTGCCCCGGTTTCGGACCTTGGCCTGGGCCGGTGCGAGCGGCCAACCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCa                                                                   <  1:3893690/150‑1 (MQ=255)
                                                                                 gCCCATTCGCGGGC‑TTGCCCGCTCCCAACGGCCTTGCCCAGGGGTCCGACCTTGTGCAGGTCCTGTGCGAGCGGGCAATCCCGCGAAAGGGCCCCCCACAAATTT‑CCCCGACTGAACAATAAACCCGCCCCACCGTCACTTATTCAGGCa                                                                   >  2:4571565/1‑150 (MQ=25)
                                                                                  cccATTCGCGGGC‑TTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCAt                                                                  <  2:3049676/150‑1 (MQ=255)
                                                                                  cccATGGGCGGGC‑TGGCCCGCCCCCCCCGGCTGTGCCGGGGGTTGGGCCCTTGTGCCGGCCGTGGGCGCGGGGGCACGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCAt                                                                  <  1:3411198/150‑1 (MQ=255)
                                                                                   ccATTCGCGGGC‑TTGCCCGCTCCCCCAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATc                                                                 <  2:4820741/150‑1 (MQ=255)
                                                                                                    ccGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCg                                                 >  2:4437681/1‑150 (MQ=18)
                                                                                                    ccGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCg                                                 >  1:4068295/1‑150 (MQ=18)
                                                                                                    ccGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCg                                                 >  2:1694587/1‑150 (MQ=18)
                                                                                                    ccGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCg                                                 >  2:1680653/1‑150 (MQ=255)
                                                                                                      gCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGtt                                               >  2:3689783/1‑150 (MQ=255)
                                                                                                         cccACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCTAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGAc                                            <  2:5103324/150‑1 (MQ=255)
                                                                                                          ccACAGGCATTGCACAGGGTTCGTACCTTGTGCAGGACCTGTGCGAGCGGGCAAGCCCTCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACt                                           <  2:427534/150‑1 (MQ=255)
                                                                                                          ccACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACt                                           >  1:885737/1‑150 (MQ=255)
                                                                                                          ccACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACt                                           >  2:5016117/1‑150 (MQ=255)
                                                                                                           cacaGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTg                                          >  2:4868502/1‑150 (MQ=255)
                                                                                                            acaGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTga                                         >  2:5221035/1‑150 (MQ=255)
                                                                                                            acaGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTga                                         >  1:45022/1‑150 (MQ=255)
                                                                                                            acaGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTga                                         >  1:3335894/1‑150 (MQ=255)
                                                                                                              aggacctGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTgaga                                       >  2:180858/7‑150 (MQ=21)
                                                                                                              agcaaTTCCAAAGTGTTCTCACATTGTCCACGACCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAACGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTgaga                                       <  2:1005684/146‑1 (MQ=255)
                                                                                                              aGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTgaga                                       >  2:740484/1‑150 (MQ=255)
                                                                                                               ggcttggCACTGGGTTGGGACCTTGTCCCGCTCGTGTGCGCGGGGGCACGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGagag                                      <  1:1344939/144‑1 (MQ=255)
                                                                                                                gctttGCCCCGGGTCCGGCCCTGGGGCCGGTCCGGGGCGCGCGGGCACGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAgg                                     <  1:2549648/147‑1 (MQ=255)
                                                                                                                gacctGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAgg                                     >  1:4507160/5‑150 (MQ=25)
                                                                                                                gCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAgg                                     <  2:3390125/150‑1 (MQ=255)
                                                                                                                gCATGGCACTGGGGCCGGCCCTTGGGCCGGGCGGGGGCGAGCGGGCACCCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAgg                                     <  1:2549657/150‑1 (MQ=255)
                                                                                                                       acaGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGtt                              >  2:2763846/1‑150 (MQ=255)
                                                                                                                       acaGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGtt                              >  1:317478/1‑150 (MQ=255)
                                                                                                                       acaGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGtt                              >  2:3656386/1‑150 (MQ=255)
                                                                                                                       acaGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGtt                              >  2:2187062/1‑150 (MQ=255)
                                                                                                                       acaGGGTACGGACCTGGTGCAGGTCCTGTGCGACCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGtt                              <  1:3099403/150‑1 (MQ=255)
                                                                                                                         cgggTTCGGCCCTTGTGCAGGACCTGTGCGCGGGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCa                            <  1:138342/149‑1 (MQ=255)
                                                                                                                         aGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCGAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCa                            >  1:3062762/1‑150 (MQ=255)
                                                                                                                           ggTTCGGCCCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAAt                          <  2:1247884/150‑1 (MQ=255)
                                                                                                                           ggTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAAt                          <  2:1247716/150‑1 (MQ=255)
                                                                                                                           ggTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAAt                          >  1:351084/1‑150 (MQ=255)
                                                                                                                           ggTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAAt                          <  2:4068295/150‑1 (MQ=37)
                                                                                                                             ttCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGa                        >  2:4528805/1‑150 (MQ=255)
                                                                                                                               cGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGaaa                      <  1:3656386/150‑1 (MQ=255)
                                                                                                                               cGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGaaa                      <  1:4490103/150‑1 (MQ=255)
                                                                                                                                ggACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAag                     >  2:4130965/1‑150 (MQ=255)
                                                                                                                                ggACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAag                     >  2:209182/1‑150 (MQ=255)
                                                                                                                                ggACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAag                     >  1:4396425/1‑150 (MQ=255)
                                                                                                                                ggACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAAT‑GGCCCCCCACAATTTTACCCGGA‑AGTACATTATCCCCCCCCCAAGGGCAGTTATTCAGGCCTCACCCCCCCTGTGTTCCGTGGATGGGGTGGTCCGTTCCATGGAagt                    >  1:1798878/1‑149 (MQ=255)
                                                                                                                                 gACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAaga                    >  2:4005938/1‑150 (MQ=255)
                                                                                                                                 gACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAaga                    >  2:744583/1‑150 (MQ=255)
                                                                                                                                 gACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAaga                    >  2:3407626/1‑150 (MQ=255)
                                                                                                                                 gACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAaga                    >  1:4193776/1‑150 (MQ=255)
                                                                                                                                 gACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAaga                    >  2:2689867/1‑150 (MQ=255)
                                                                                                                                   ccTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGgt                                                      <  2:3574147/114‑1 (MQ=255)
                                                                                                                                   ccTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGgt                                                      >  1:3574147/1‑114 (MQ=255)
                                                                                                                                    cTTGTGCAGGTCCTGTGCGAGCGGGTAAGCCCGCGACTGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATACCCGACCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCTTTGACTGAGAGGGTCCGTTCAATGAAAGAGTg                 >  2:4348839/1‑150 (MQ=255)
                                                                                                                                    cTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGtt                                               >  2:4146071/1‑120 (MQ=255)
                                                                                                                                    cTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGtt                                               <  1:4146071/120‑1 (MQ=255)
                                                                                                                                     ttGTGCAGGTCCTGTGCGCGCGGGCATGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGTgg                <  1:4016366/150‑1 (MQ=255)
                                                                                                                                     ttGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGTgg                >  2:743248/1‑150 (MQ=255)
                                                                                                                                        ttcAGGTGGTGTGCGCGGGGCAAACCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGTGGGaa             <  1:2784536/148‑1 (MQ=255)
                                                                                                                                        tgCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGTGGGaa             <  2:4505081/150‑1 (MQ=255)
                                                                                                                                        tgCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGTGGGaa             <  2:4444158/150‑1 (MQ=255)
                                                                                                                                          cAGGTCCTGTGCGAGCGGGCAAGCCCGCGACTGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGTGGGAAGt           >  2:133156/1‑150 (MQ=255)
                                                                                                                                          cAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGGGGGAAGt           >  2:5389758/1‑150 (MQ=255)
                                                                                                                                             gTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGTGGGAAGTCAt        >  2:96857/1‑150 (MQ=255)
                                                                                                                                              tCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGTGGGa              >  1:2567601/1‑143 (MQ=255)
                                                                                                                                              tCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGTGGGa              <  2:2567601/143‑1 (MQ=255)
                                                                                                                                               ccTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTgaga                                       >  1:3324212/1‑117 (MQ=255)
                                                                                                                                               ccTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTgaga                                       <  2:3324212/117‑1 (MQ=255)
                                                                                                                                                 tgtgCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGACCATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGTGGGAAGTCAttttt    >  2:3042284/1‑150 (MQ=255)
                                                                                                                                                 ggtgCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGTGGGAAGTCAttttt    <  2:3458291/149‑1 (MQ=255)
                                                                                                                                                   tgCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGTGGGAAGTCATTTTTGc  >  1:3276146/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGC‑TTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTT‑CCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGTGGGAAGTCATTTTTGC  >  NC_002947/2065045‑2065340

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: ‑