Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,172,004 G→C 6.4% R78P (CGG→CCG)  PP_2784 → short‑chain dehydrogenase/reductase family oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,172,0040GC 6.4% 296.9 / 9.1 95R78P (CGG→CCG) PP_2784short‑chain dehydrogenase/reductase family oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (45/42);  new base C (3/3);  total (49/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

GAGCACCTTTTGATCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTTGTGCCTGCCACTGATGGCGGATG  >  NC_002947/3171856‑3172153
                                                                                                                                                     |                                                                                                                                                     
gAGCACCTTTTGATCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCgg                                                                                                                                                      >  2:2346714/1‑150 (MQ=255)
 aGCACCTTTTGATCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCggg                                                                                                                                                     <  1:847633/150‑1 (MQ=255)
 aGCACCTTTTGATCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCggg                                                                                                                                                     >  2:1847870/1‑150 (MQ=255)
 aGCACCTTTTGATCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCggg                                                                                                                                                     <  1:1720170/150‑1 (MQ=255)
  gCACCTTTTGATCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGt                                                                                                                                                    <  1:2346714/150‑1 (MQ=255)
    aCCTTTTGATCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCg                                                                                                                                                  >  1:1867837/1‑150 (MQ=255)
     ccTTTTGATCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCgg                                                                                                                                                 >  2:2085755/1‑150 (MQ=255)
       ttttGATCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCt                                                                                                                                               >  2:533279/1‑150 (MQ=255)
        tttGATCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCtt                                                                                                                                              <  1:89390/150‑1 (MQ=255)
        tttGATCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCCGGTCGGcca                                                                                                                                              >  2:1444079/1‑148 (MQ=255)
          tGATCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCCGGTCGGccagg                                                                                                                                            <  1:1444079/150‑5 (MQ=255)
            aTCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTc                                                                                                                                          >  2:2047429/1‑150 (MQ=255)
                  cttgATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTAGCCGACCGGGTCGGCTTTGTCGAtgtg                                                                                                                                    <  2:1548919/150‑1 (MQ=255)
                     gATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTg                                                                                                                                 <  1:2047429/150‑1 (MQ=255)
                     gATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCAAACCCGCTACTGCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTg                                                                                                                                 <  1:1376632/150‑1 (MQ=255)
                               tGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAgcg                                                                                                                       >  1:966486/1‑150 (MQ=255)
                               tGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAgcg                                                                                                                       >  2:493028/1‑150 (MQ=255)
                                  cAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCacaag                                                                                                                         <  2:1010979/145‑5 (MQ=255)
                                  cAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAg                                                                                                                         >  1:1010979/1‑145 (MQ=255)
                                  cAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCg                                                                                                                    <  1:1847870/150‑1 (MQ=255)
                                   aaCCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATTAAGCGAACATCGCCGGTGCCGCCAGCCTGCAAGTTGTGCCCTAGGATCGGACGAGCCGGGGCGGCTTTGTCGATGTGCTGGTCAACAGCGCCgg                                                                                                                   >  1:2341872/1‑150 (MQ=255)
                                   aaCCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTg                                                                                                                                            >  1:1809609/1‑125 (MQ=255)
                                   aaCCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTg                                                                                                                                            <  2:1809609/125‑1 (MQ=255)
                                   aaCCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCgg                                                                                                                   >  2:2237241/1‑150 (MQ=255)
                                   aaCCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCgg                                                                                                                   <  1:2237241/150‑1 (MQ=255)
                                   aaCCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCCATGACTTCAGCCTGCACATTATGTTCGGGCTTCTGTTCGAGGTGTGCGGCTTTGTCGATGTGCTGGTCAACAGCGCCgg                                                                                                                   >  1:2341923/1‑150 (MQ=255)
                                    aCCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCggg                                                                                                                  <  2:1534160/150‑1 (MQ=255)
                                      cTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCggggg                                                                                                                >  2:1209689/1‑149 (MQ=255)
                                      cTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGt                                                                                                                >  2:790297/1‑150 (MQ=255)
                                      cTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGt                                                                                                                <  1:493028/150‑1 (MQ=255)
                                       ttGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTg                                                                                                                                 >  2:377931/1‑132 (MQ=255)
                                       ttGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTg                                                                                                                                 <  1:377931/132‑1 (MQ=255)
                                         gCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTg                                                                                                             <  1:1007546/150‑1 (MQ=255)
                                         gCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTg                                                                                                             <  1:956403/150‑1 (MQ=255)
                                            ggggTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGcga                                                                                                          >  2:763195/1‑150 (MQ=255)
                                             gggtgggtCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAAGATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGcgac                                                                                                         >  2:1311092/1‑150 (MQ=255)
                                               gtgggtCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCttt                                                                                                                                             <  2:614937/112‑1 (MQ=255)
                                               gtgggtCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCttt                                                                                                                                             >  1:614937/1‑112 (MQ=255)
                                               gtgggtCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGcgacga                                                                                                       >  2:1929566/1‑150 (MQ=255)
                                               gtgggtCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGcgacga                                                                                                       <  1:110717/150‑1 (MQ=255)
                                                                aaTTTGCCGCTCGTATCGCAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGAcc                                                                                      <  1:1929566/150‑1 (MQ=255)
                                                                  tttGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTc                                                                                    <  1:2062397/150‑1 (MQ=255)
                                                                         cTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACt                                                                             >  2:2257873/1‑150 (MQ=255)
                                                                          tCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCggg                                                                                                                  <  1:930350/112‑1 (MQ=255)
                                                                          tCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCggg                                                                                                                  >  2:930350/1‑112 (MQ=255)
                                                                          tCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTAGACGACCGGGTCGGCTTTGTTTATTTGGTGGTCTTCTTCGTTGGTTTCTTCGACGAGAACGAACCCGAAGACCTCGACAACTg                                                                            >  2:272571/1‑150 (MQ=255)
                                                                            gTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGc                                                                          >  2:1858757/1‑150 (MQ=255)
                                                                             tATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGca                                                                         >  2:1591694/1‑150 (MQ=255)
                                                                                        aTCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTTGCCGACCGGGTCGGCTGTGTGGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTAAttt                                                              >  1:1066901/1‑150 (MQ=255)
                                                                                        aTCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCAttt                                                              <  2:428387/150‑1 (MQ=255)
                                                                                        aTCAAGCGAACATCGCCGACCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCAttt                                                              <  1:533279/150‑1 (MQ=255)
                                                                                         tCAAGCGAACATCGCCGATCTGGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGGTGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGTGGTCTGCGCCGAGCCCGAACGCGAAGACCTCGACAACTGGCACTAGGTCATTTc                                                             >  2:48552/1‑150 (MQ=255)
                                                                                           aaGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAgaacgaac                                                                                               <  2:1829725/114‑1 (MQ=255)
                                                                                           aaGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAgaacgaac                                                                                               >  1:1829725/1‑114 (MQ=255)
                                                                                           aaGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGa                                                           <  1:1858757/150‑1 (MQ=255)
                                                                                           aaGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGa                                                           <  2:733035/150‑1 (MQ=255)
                                                                                           aaGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGAATTGGCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGa                                                           <  1:763195/150‑1 (MQ=255)
                                                                                            aGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGAt                                                          >  1:2278561/1‑150 (MQ=255)
                                                                                                aaCATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAAc                                                      <  1:2085755/150‑1 (MQ=255)
                                                                                                aaCATCGCCGATCGCGCCAGCCTGCAACTTTTGGTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCAGGGTCTGCGACTAGAACCAGCCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAAc                                                      <  1:48552/150‑1 (MQ=255)
                                                                                                      gCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAAc                                                >  2:1482873/1‑150 (MQ=255)
                                                                                                        caaTCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAAc                                                                                                                           <  1:2229161/71‑1 (MQ=255)
                                                                                                        caaTCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAAc                                                                                                                           >  2:2229161/3‑73 (MQ=255)
                                                                                                        cGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGcc                                                                                                                     <  2:610480/79‑1 (MQ=255)
                                                                                                        cGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGcc                                                                                                                     >  1:610480/1‑79 (MQ=255)
                                                                                                               gCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGcgacga                                                                                                       <  2:541256/86‑1 (MQ=255)
                                                                                                               gCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGcgacga                                                                                                       >  1:541256/1‑86 (MQ=255)
                                                                                                                ccAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGaa                                                                                          >  2:2140683/1‑98 (MQ=255)
                                                                                                                ccAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGaa                                                                                          <  1:2140683/98‑1 (MQ=255)
                                                                                                                 cAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCggg                                                                                                                  >  2:1495327/1‑73 (MQ=255)
                                                                                                                 cAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCggg                                                                                                                  <  1:1495327/73‑1 (MQ=255)
                                                                                                                 cAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTa                                     <  1:2257873/150‑1 (MQ=255)
                                                                                                                     cTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTa                                     <  2:1085633/146‑1 (MQ=255)
                                                                                                                     cTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTa                                     >  1:1085633/1‑146 (MQ=255)
                                                                                                                       gCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGcc                                                                                                                     >  1:582128/1‑64 (MQ=255)
                                                                                                                       gCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGcc                                                                                                                     <  2:582128/64‑1 (MQ=255)
                                                                                                                        cAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCAcc                              <  1:1482873/150‑1 (MQ=255)
                                                                                                                        cAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCAAGACGTTGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCAcc                              >  2:676381/1‑150 (MQ=255)
                                                                                                                        cAACTTCTGTTCAAGGGCCGGGCCGAACAAGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCAcc                              <  1:1209689/150‑1 (MQ=255)
                                                                                                                         aaCTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGaa                                                                                          >  1:422103/1‑89 (MQ=255)
                                                                                                                         aaCTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGaa                                                                                          <  2:422103/89‑1 (MQ=255)
                                                                                                                         aaCTTCTGTTCAAGGGCCTGGCCGACAGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCt                             >  2:2170211/1‑150 (MQ=255)
                                                                                                                            ttCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGt                          >  2:2250386/1‑150 (MQ=255)
                                                                                                                               tGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTtgt                       <  1:1591694/150‑1 (MQ=255)
                                                                                                                                gTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTtgtg                      >  1:609638/1‑150 (MQ=255)
                                                                                                                                gTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTtgtg                      >  1:72747/1‑150 (MQ=255)
                                                                                                                                    aaGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTTGtgcctg                  >  1:329857/1‑150 (MQ=255)
                                                                                                                                    aaGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTTGtgcctg                  >  1:2031207/1‑150 (MQ=255)
                                                                                                                                        gCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCAACGAGAACGAACCAGCCGAAGGGGACAACTGGGACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTTGTGCCTGCCAc              >  1:1495940/1‑150 (MQ=255)
                                                                                                                                        acaaagCCGACCCGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTtgt                       <  1:112394/136‑1 (MQ=255)
                                                                                                                                        acaaagCCGACCCGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTtgt                       >  2:112394/6‑141 (MQ=255)
                                                                                                                                          cTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTAc                                    <  1:1272725/126‑1 (MQ=255)
                                                                                                                                          cTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTAc                                    >  2:1272725/1‑126 (MQ=255)
                                                                                                                                          aaagCCGACCCGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGTGTCTGCGACGAGAACGAAcc                                                                                              <  1:20892/65‑1 (MQ=255)
                                                                                                                                          aaagCCGACCCGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGTGTCTGCGACGAGAACGAAcc                                                                                              >  2:20892/4‑68 (MQ=255)
                                                                                                                                            ggCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTTGTGCCTGCCACTGAt          <  1:669450/150‑1 (MQ=255)
                                                                                                                                              ccGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTTGTGCCTGCCACTGATgg        <  1:790297/150‑1 (MQ=255)
                                                                                                                                               cGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTTGTGCCTGCCACTGATGGc       <  2:966486/150‑1 (MQ=255)
                                                                                                                                                gACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTTGTGCCTGCCACTGATGGCg      >  2:16263/1‑150 (MQ=255)
                                                                                                                                                     gggTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTTGTGCCTGCCACTGATGGCGGAtg  <  2:2031207/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAGCACCTTTTGATCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTTGTGCCTGCCACTGATGGCGGATG  >  NC_002947/3171856‑3172153

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: