Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,163,542 T→C 6.1% intergenic (‑37/‑145) pssA ← / → PP_3665 CDP‑diacylglycerol‑‑serine O‑phosphatidyltransferase/AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,163,5420TC 6.1% 240.3 / 4.2 82intergenic (‑37/‑145)pssA/PP_3665CDP‑diacylglycerol‑‑serine O‑phosphatidyltransferase/AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (36/41);  new base C (2/3);  total (38/44)
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.97e-01

GATGCGCTCGAACAGGCAGGCGCGGTAAGCCCTGGCATCGGGCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGCCTCAGCTTGCACAGGAACAGCGCAATGGT  >  NC_002947/4163395‑4163691
                                                                                                                                                   |                                                                                                                                                       
gATGCGCTCGAACAGGCAGGCGCGGTAAGCCCTGGCATCGGGCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCcctgc                                                                                                                                                       <  1:792727/150‑1 (MQ=255)
    cgcTCGAACAGGCAGGCGCGGTAAGCCCTGGCATCGGGCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCctgcctgc                                                                                                                                                   <  1:978413/150‑1 (MQ=255)
       tCGAACAGGCAGGCGCGGTAAGCCCTGGCATCGGGCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCcgg                                                                                                                                                >  2:1238933/1‑150 (MQ=255)
             aggcaggcGCGGTAAGCCCTGGCATCGGGCAGCCCCTGGATCGCTTCGGGGGCCAGCTCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCCGCCGGCAGGCAGGGc                                                                                                                                          <  1:1454174/150‑1 (MQ=255)
             aggcaggcGCGGTAAGCCCTGGCATCGGGCAGCACCTGGATCGCTTCGGGAGCGCTGGGGTAGCGTGGCAAGCCCGACAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCCGCCGGCAGGCAGGGc                                                                                                                                          <  2:731324/150‑1 (MQ=255)
               gcaggcGCGGTAAGCCCTGGCATCGGGCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCggg                                                                                                                                            <  1:1392885/146‑1 (MQ=255)
               gcaggcGCGGTAAGCCCTGGCATCGGGCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCggg                                                                                                                                            >  2:1392885/1‑146 (MQ=255)
                        gTAAGCCCTGGCATCGGGCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGcc                                                                                                                               <  2:180850/150‑1 (MQ=255)
                        gTAAGCCCTGGCATCGGGCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGcc                                                                                                                               <  2:177332/150‑1 (MQ=255)
                            gCCCTGGCATCGGGCAGCACCTGGAGCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGTGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTg                                                                                                                           <  2:561897/150‑1 (MQ=255)
                             cccTGGCATCGGGCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGc                                                                                                                          >  1:121742/1‑150 (MQ=255)
                             cccTGGCAGCGGGCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGCTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGc                                                                                                                          >  2:1538479/1‑150 (MQ=255)
                               cTGGCATCGGGCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCAt                                                                                                                        >  1:23907/1‑150 (MQ=255)
                                    aTCGGGCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCt                                                                                                                                   >  1:1118526/1‑134 (MQ=255)
                                    aTCGGGCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCt                                                                                                                                   <  2:1118526/134‑1 (MQ=255)
                                    aTCGGGCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCggg                                                                                                                   <  1:597170/150‑1 (MQ=255)
                                     tCGGGCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCcctgc                                                                                                                                                       >  2:610383/1‑113 (MQ=255)
                                     tCGGGCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCcctgc                                                                                                                                                       <  1:610383/113‑1 (MQ=255)
                                     tCGGGCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTgg                                                                                                                                 <  1:30847/135‑1 (MQ=255)
                                     tCGGGCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTgg                                                                                                                                 >  2:30847/1‑135 (MQ=255)
                                     tCGGGCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCgggg                                                                                                                  >  1:996926/1‑150 (MQ=255)
                                      cGGGCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGgt                                                                                                                 <  2:342216/150‑1 (MQ=255)
                                      cGGGCAGCACCTGGATCGCTGCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTCCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGgt                                                                                                                 <  2:5923/150‑1 (MQ=255)
                                       gggCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGgtg                                                                                                                >  1:1008667/1‑150 (MQ=255)
                                            gcaCCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGa                                                                                                         <  2:1131933/150‑1 (MQ=255)
                                            gcaCCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGa                                                                                                         <  2:1131940/150‑1 (MQ=255)
                                              aCCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATc                                                                                                       <  1:1213834/150‑1 (MQ=255)
                                              aCCTGGATCGCTTCGGGGGCCAGGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACCGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATc                                                                                                       <  2:301698/150‑1 (MQ=255)
                                                    aTCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCg                                                                                                 >  2:1925947/1‑150 (MQ=255)
                                                        cTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCaa                                                                                             >  1:664268/1‑150 (MQ=255)
                                                         ttCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAAc                                                                                            >  1:1637057/1‑150 (MQ=255)
                                                          ttgggggCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACg                                                                                           >  1:878923/3‑150 (MQ=255)
                                                              gggCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGAccc                                                                                       <  2:1008667/150‑1 (MQ=255)
                                                                 ccATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATa                                                                                    >  2:109431/1‑150 (MQ=255)
                                                                 ccATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATa                                                                                    <  1:109431/150‑1 (MQ=255)
                                                                  taTGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATaa                                                                                   <  2:1126378/149‑1 (MQ=255)
                                                                  cATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATaa                                                                                   >  1:455908/1‑150 (MQ=255)
                                                                  cATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATaa                                                                                   <  2:1510498/150‑1 (MQ=255)
                                                                    tGGCGTAGCCTGGTAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATaaca                                                                                 >  2:1759234/1‑150 (MQ=255)
                                                                    tGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGCGTATGATCGACCCGGCAACGACCCATaaca                                                                                 >  2:1325141/1‑150 (MQ=255)
                                                                      gCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATaacaac                                                                               <  1:722523/150‑1 (MQ=255)
                                                                        gTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAAcaacaa                                                                             <  1:838861/150‑1 (MQ=255)
                                                                         tAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAg                                                                            <  2:1362886/150‑1 (MQ=255)
                                                                         tAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATCCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAg                                                                            <  1:233844/150‑1 (MQ=255)
                                                                           gCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGaa                                                                          >  1:1042740/1‑150 (MQ=255)
                                                                           gCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGaa                                                                          <  2:121742/150‑1 (MQ=255)
                                                                              tGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCa                                                                       >  2:1724974/1‑150 (MQ=255)
                                                                              tGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCa                                                                       >  2:749040/1‑150 (MQ=255)
                                                                                     cccGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGaa                                                                          <  2:1447724/140‑1 (MQ=255)
                                                                                     cccGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGaa                                                                          >  1:1447724/1‑140 (MQ=255)
                                                                                     cccGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCCGCCGGCAGGCAGGGc                                                                                                                                          <  2:840186/78‑1 (MQ=38)
                                                                                     cccGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCCGCCGGCAGGCAGGGc                                                                                                                                          >  1:840186/1‑78 (MQ=38)
                                                                                         ccAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGGCGTCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGctct                                                            >  1:1079228/1‑150 (MQ=255)
                                                                                           aGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGc                                                          >  1:1659844/1‑150 (MQ=255)
                                                                                            gCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCt                                                         <  2:1659844/150‑1 (MQ=255)
                                                                                               tGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTtc                                                      <  1:1238933/150‑1 (MQ=255)
                                                                                                ggTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTtct                                                     <  2:23907/150‑1 (MQ=255)
                                                                                                  tttCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCAG‑‑ACGTCTGAACTCCAGTCAAGCTAGCCGGCAACGACCCATAACAACAAGAACCACGCAAGGCTCTGCTTTCTCa                                                   >  1:431995/1‑150 (MQ=255)
                                                                                                      ggggTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATaacaac                                                                               >  1:73182/1‑118 (MQ=255)
                                                                                                      ggggTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATaacaac                                                                               <  2:73182/118‑1 (MQ=255)
                                                                                                        ggtgaaCTGATGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGccc                                             >  2:1013342/6‑150 (MQ=255)
                                                                                                          ttCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCt                                           <  1:1282475/150‑1 (MQ=255)
                                                                                                            cACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGa                                         <  2:621076/150‑1 (MQ=255)
                                                                                                            cACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGa                                         <  1:1925947/150‑1 (MQ=255)
                                                                                                             aCTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGAt                                                                                                        <  2:1313704/86‑1 (MQ=255)
                                                                                                             aCTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGAt                                                                                                        >  1:1313704/1‑86 (MQ=255)
                                                                                                               tGACGCTCACCTTCTCTTTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTtct                                                     <  2:240667/135‑1 (MQ=255)
                                                                                                               tGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTtct                                                     >  1:240667/1‑135 (MQ=255)
                                                                                                                   gCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGaccaacca                                  >  2:248527/1‑150 (MQ=255)
                                                                                                                     tCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAgg                                >  1:1212876/1‑150 (MQ=255)
                                                                                                                      cACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCGTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGc                               <  2:996926/150‑1 (MQ=255)
                                                                                                                         tgttcctGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGCCTc                            >  2:6820/6‑150 (MQ=255)
                                                                                                                         cTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGCCTc                            <  1:248527/150‑1 (MQ=255)
                                                                                                                              ctGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCAGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGCCTCAGCtt                       >  1:645937/1‑150 (MQ=255)
                                                                                                                               tgtgCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGCCTCAGCTTg                      <  2:1042740/150‑1 (MQ=255)
                                                                                                                                    aaTCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTtct                                                     <  1:1186654/114‑1 (MQ=255)
                                                                                                                                    aaTCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTtct                                                     >  2:1186654/1‑114 (MQ=255)
                                                                                                                                    aaTCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGCCTCAGCTTGCACAg                 >  1:1840610/1‑150 (MQ=255)
                                                                                                                                       cACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGCCTCAGCTTGCACAGGaa              <  2:1840610/150‑1 (MQ=255)
                                                                                                                                         cTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGCCTCAGCTTGCACAGGAACa            <  1:6820/150‑1 (MQ=255)
                                                                                                                                          tGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGCCt                             <  2:482345/132‑1 (MQ=255)
                                                                                                                                          tGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGCCt                             >  1:482345/1‑132 (MQ=255)
                                                                                                                                          tGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGCCTCAGCTTGCACAGGAACAg           <  2:878923/150‑1 (MQ=255)
                                                                                                                                            tCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGaccaacca                                  >  1:1753293/1‑125 (MQ=255)
                                                                                                                                            tCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGaccaacca                                  <  2:1753293/125‑1 (MQ=255)
                                                                                                                                              ggCCCCGCCGGCAGGCAGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGCCTCAGCTTGCACAGGAACAGCGCa       >  2:1251229/1‑150 (MQ=255)
                                                                                                                                                 cctgcctgCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGCCTCAGCTTGCACAGGAACAGCGCAATg    <  2:431995/150‑1 (MQ=255)
                                                                                                                                                   tgcctgccGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCCACGACCCATAACACCAAGAAGGAGGGGAGGGTCTACTAGCACACGACCCCTACCCACCCAGCCCCCAGCTTGCACAGGCACACCGCCACGGt  >  1:697904/1‑150 (MQ=255)
                                                                                                                                                   tgcctgccGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGCGTATGAACGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACGTCCCCTGACCAACCAGGCCTCAGCTTGCACAGGAACAGCGCAATGGt  >  2:1298215/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                       
GATGCGCTCGAACAGGCAGGCGCGGTAAGCCCTGGCATCGGGCAGCACCTGGATCGCTTCGGGGGCCATGGCGTAGCCTGGCAAGCCCGCCAGCATGGTTTCGGGGTTCACTGACGCTCACCTTCTCTGTGCAATCACTGTCGGCCCTGCCTGCCGGCGGGGCCATTCCTGGCCGCTGCATCCGGGGTGTATGATCGACCCGGCAACGACCCATAACAACAAGAACCACGCGAGGCTCTGCTTTCTCACTGCCCCTGACCAACCAGGCCTCAGCTTGCACAGGAACAGCGCAATGGT  >  NC_002947/4163395‑4163691

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