Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,844,863 Δ1 bp 26.0% intergenic (+40/‑257) ccoP‑II → / → PP_4259 cbb3‑type cytochrome c oxidase subunit/iron‑sulfur cluster‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,844,8630C.26.0% 333.2 / 135.7 131intergenic (+40/‑257)ccoP‑II/PP_4259cbb3‑type cytochrome c oxidase subunit/iron‑sulfur cluster‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (47/50);  new base . (16/18);  total (63/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 141 total aligned reads displayed.

AGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCA  >  NC_002947/4844714‑4845011
                                                                                                                                                      |                                                                                                                                                      
aGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTacac                                                                                                                                                        >  2:4503916/1‑150 (MQ=255)
 ggggCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTAcaca                                                                                                                                                       <  1:2986992/150‑1 (MQ=255)
       gATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATcc                                                                                                                                                 <  1:2992001/150‑1 (MQ=255)
         tGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCtt                                                                                                                                              <  1:1580252/150‑1 (MQ=255)
         tGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCtt                                                                                                                                              <  2:3666037/150‑1 (MQ=255)
         tGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCtt                                                                                                                                              <  2:3130815/150‑1 (MQ=255)
         tGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCtgt                                                                                                                                              <  1:2938977/150‑3 (MQ=255)
            cGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTa                                                                                                                                           <  1:2583585/150‑1 (MQ=255)
             ggCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAg                                                                                                                                          <  1:1552806/150‑1 (MQ=255)
             ggCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAg                                                                                                                                          <  1:28175/150‑1 (MQ=255)
             ggCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAg                                                                                                                                          <  2:4283429/150‑1 (MQ=255)
             cgCCCAGGAACAACTGCAGGGCAACGACCAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTACGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAg                                                                                                                                          <  2:1338211/149‑1 (MQ=38)
              gCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGCCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAgg                                                                                                                                         <  2:4446517/150‑1 (MQ=37)
               cccAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGa                                                                                                                                        >  1:4495269/1‑150 (MQ=255)
               cccAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGa                                                                                                                                        <  1:31292/150‑1 (MQ=255)
               cccAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGa                                                                                                                                        <  1:230378/150‑1 (MQ=255)
                ccAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAg                                                                                                                                       <  1:1755031/150‑1 (MQ=255)
                 cAGGAACAACTGCAGGGCATCGACACGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGTAGc                                                                                                                                      <  2:3009011/150‑1 (MQ=38)
                 cAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGc                                                                                                                                      <  1:2604673/150‑1 (MQ=255)
                   tgAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCgg                                                                                                                                    <  1:2047933/149‑1 (MQ=255)
                   ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCgg                                                                                                                                    <  1:367324/150‑1 (MQ=255)
                   ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCgg                                                                                                                                    <  2:3314095/150‑1 (MQ=255)
                   ggAACAACTGCAGGGCAACGACAAGGTGCACGTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCgg                                                                                                                                    <  1:3332764/150‑1 (MQ=255)
                    gAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCggg                                                                                                                                   <  1:1098867/150‑1 (MQ=255)
                    gAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCggg                                                                                                                                   <  1:591126/150‑1 (MQ=255)
                    gAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCggg                                                                                                                                   <  1:1470720/150‑1 (MQ=255)
                    gAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCggg                                                                                                                                   <  1:1098773/150‑1 (MQ=255)
                    gAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCGTGTAGGAGCggg                                                                                                                                   <  1:3898879/150‑1 (MQ=255)
                    gAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCggg                                                                                                                                   <  1:1531114/150‑1 (MQ=255)
                     aacaacTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGt                                                                                                                                  <  2:3056698/150‑1 (MQ=255)
                     aacaacTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGt                                                                                                                                  <  2:775244/150‑1 (MQ=255)
                        aacTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCGGGTAGGAGCGGGTTTa                                                                                                                               >  2:3719461/1‑150 (MQ=255)
                         acTGCAGGGCAACGACAAGGTGCACTTGTTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTAc                                                                                                                             <  1:3777421/150‑1 (MQ=255)
                         acTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTAc                                                                                                                             <  1:305213/150‑1 (MQ=255)
                           tGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTAccc                                                                                                                           >  2:1864677/1‑150 (MQ=255)
                             cAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGTGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCcgc                                                                                                                         <  1:2455747/150‑1 (MQ=255)
                                   aaCGCCACGGTGCACTGGCTGGCGGCTTACGTCTACCGCCTGGCGCACCCGGCTGCGCCGGCCAAGGCTGCGTAAGCCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTCCACGGATCCTGTGTAGGAGCGGGTTTACCCGCGAagag                                                                                                                   <  2:2743558/150‑1 (MQ=21)
                                      caCAAGGTGCACTTGCTGGCGGCTTCCGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTACGGCAGGCCTGGCCACTTCGCGGGTAACCCCGCCCCTACACCGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGcc                                                                                                                <  2:1378092/149‑1 (MQ=38)
                                       ccatggTGCACTTGCTGGCGGCTTACGTCTACCGCCTGTCGCACCCGGCTGAGCCGGCCACGGCGGCGTAAGGCAGGCCGGGCCTCTTCGCGGGTAAACCCGCTCCTACACCGGATCCTTGTTGGAGGGGGTTTACCCGCGAAGAGGCCg                                                                                                               <  2:3167391/146‑1 (MQ=25)
                                                 cTTGGTGGCGGCTTACGTCTACCGCCTGTCGCACCTGGCTGCGCCGGCCAAGGGTGAGTACGGCTGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACCGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGcac                                                                                                     <  2:3937669/150‑1 (MQ=25)
                                                     cTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACCGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCa                                                                                                 >  2:675550/1‑150 (MQ=255)
                                                        gcggcTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCACCTACACACGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAATGGGCCCGCAttt                                                                                              >  1:2867593/1‑150 (MQ=37)
                                                          ggcTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCTGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACCCAGGATCCTTGTCGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTtct                                                                                            <  2:374787/150‑1 (MQ=255)
                                                             ttACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGc                                                                                         <  2:4067581/150‑1 (MQ=255)
                                                              tACGTCTACCGCCTGTCGCACCCGGCTGAGCCGGCCACGGCTTATTAAGGCCGGCCTGGCCTCTTCGGGGTTAAACCCGCTCCTACACAGGATCCTGGTCGGAGGGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCa                                                                                        <  2:1653988/150‑1 (MQ=37)
                                                               aCGTCTACAGCCGGTCGCACCTGGCTGAGCCGGCCACGGCTGAGTAAGGCAGGCCTGGCCTCTGCGCGGGTAAACCCGCGCCTCCACCCGATCCTTGTAGGAGCTGGTTAACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAg                                                                                       <  2:3316981/150‑1 (MQ=38)
                                                                 gTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACCCAGGATCCGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGc                                                                                     >  2:1370749/1‑150 (MQ=255)
                                                                 gTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGc                                                                                     >  2:3647837/1‑150 (MQ=255)
                                                                     tcAGCCTGGTGCGACAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCCCGCCTACCCAGGATCCGTGTAGGAGCGGGTTTAACCCCCAAGAGGCCCGGGCTGGCCCCCCAATTAACGCAGGCaact                                                                                 >  1:4368074/2‑148 (MQ=19)
                                                                       aGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCGACACAGGATCCGTGTGGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGGACAGCATTTCTCGCAGGCAAGTGa                                                                               >  2:2045028/1‑150 (MQ=255)
                                                                           tGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACCCAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGCCCTg                                                                           >  2:2641226/1‑150 (MQ=255)
                                                                              cGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGAt                                                                        >  2:1488406/1‑150 (MQ=255)
                                                                                 aCCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACCCAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCCAGTGACCGGGATCaa                                                                     >  2:668906/1‑150 (MQ=255)
                                                                                  ccAGGCGGAGCCGGCCACGGCGGAGTAAGGCCGGCCTGGCCTCTTCGCGGGAAAACCCGCTCCTACACAGGATCCTTGTTGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAAt                                                                    <  2:4222639/150‑1 (MQ=37)
                                                                                   cAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAACCCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAAtt                                                                   <  2:1460076/150‑1 (MQ=255)
                                                                                           gCCGGCCACGGCTGAGTAAGGCTGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCGTGTAGGAGGGGTTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCg                                                           <  2:470012/150‑1 (MQ=38)
                                                                                                 catggCTTAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATaa                                                     <  2:4369071/147‑1 (MQ=255)
                                                                                                    ggCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAacac                                                  >  2:4107017/1‑150 (MQ=255)
                                                                                                       tGAGTACGGCTGGCGTGGCCTGTTCGCGGGTAAACCCGCTCCTACACCGGATCCTGGTTGGCGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGtt                                               <  2:4364608/150‑1 (MQ=37)
                                                                                                       ggtgTACGGCGGGCCGGGCCTGTTCGGGGGAACCCCCGCTCCTACACTGGATCGTTGGCGGAGGGGGTTTCCCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGtt                                               <  1:3266645/147‑1 (MQ=255)
                                                                                                             cggcaggcCTGGCCTCTTCGCGGGTAACCCCGCTCCAACACAGGATCCTTGTAGGAGGGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATAc                                         <  1:4083479/149‑1 (MQ=255)
                                                                                                             aggcaggcCTGGCCTCTTCGCGGGTTAAACCGCACCAACACAGGAAGCATGTAGGAGCCGGGTTAACCGCCAAGAGGGCCGGACTGGCACAGCCATTCTCGCCGGCACGTGACCGGGGTCAATTGTCACCCCCCCGGAACCGTTTCATAc                                         >  2:1913423/1‑150 (MQ=255)
                                                                                                               gctggcgtggcctGTTGGGGGTAAACCCCGCCCCCACACTGTTTCTTTGTTGGTGGGGTTTACCCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCa                                       <  1:902225/147‑1 (MQ=255)
                                                                                                               gccggcctggcctCTTCGGGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGTGGGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCa                                       <  1:2196980/147‑1 (MQ=255)
                                                                                                                           tctTCGCGGGTAAACCCCCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCACAGCAATTCTCGCAGGCAAGTGACCTGGATCAAATGACACCCGCCAAAACACGTTTCATAACACGAAACCAAcgc                           >  2:755969/1‑150 (MQ=37)
                                                                                                                                gcgGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTa                      >  1:1227744/1‑150 (MQ=255)
                                                                                                                                    gTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAgg                  <  2:1111667/150‑1 (MQ=255)
                                                                                                                                    gTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAgg                  >  1:3919481/1‑150 (MQ=255)
                                                                                                                                    gTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAgg                  >  2:419490/1‑150 (MQ=255)
                                                                                                                                       aaCCCGCACCTACACAGGATCCTTGTCGGAGGGGTTTAACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCg               <  1:692983/150‑1 (MQ=255)
                                                                                                                                          cccGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCa             >  2:4002221/1‑150 (MQ=255)
                                                                                                                                          cccGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCa             >  1:3367981/1‑150 (MQ=255)
                                                                                                                                          cccGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCa             >  1:1594796/1‑150 (MQ=255)
                                                                                                                                           ccGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAAtt                                                                   <  1:4112998/95‑1 (MQ=255)
                                                                                                                                           ccGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAg            >  1:941041/1‑150 (MQ=255)
                                                                                                                                           ccGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAg            >  1:1951683/1‑150 (MQ=255)
                                                                                                                                           ccGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAg            >  2:1194387/1‑150 (MQ=255)
                                                                                                                                           ccGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAg            >  2:3497221/1‑150 (MQ=255)
                                                                                                                                           ccGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAg            >  2:1320232/1‑150 (MQ=255)
                                                                                                                                           ccGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAg            >  1:277099/1‑150 (MQ=255)
                                                                                                                                           ccGCTCCTACAAGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAg            >  1:732770/1‑150 (MQ=255)
                                                                                                                                           ccGCTCCTACAAGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAg            >  1:1965946/1‑150 (MQ=255)
                                                                                                                                           ccGCTCCTACAAGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAg            >  2:1587192/1‑150 (MQ=255)
                                                                                                                                           ccGCTCCTACAAGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAg            >  1:2280503/1‑150 (MQ=255)
                                                                                                                                           ccGCTCCTACAAGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAg            >  1:2956584/1‑150 (MQ=255)
                                                                                                                                           ccGCTCCTACAAGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAg            >  1:1965942/1‑150 (MQ=255)
                                                                                                                                            cGCTCCTACAAGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGt           >  2:2267062/1‑150 (MQ=255)
                                                                                                                                            cGCTCCTACAAGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGt           >  2:2531386/1‑150 (MQ=255)
                                                                                                                                             gCTCCTACAAGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTg          >  2:1255103/1‑150 (MQ=255)
                                                                                                                                             gCTCCTACAAGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTg          >  1:2387635/1‑150 (MQ=255)
                                                                                                                                               tCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCg        >  1:3907178/1‑150 (MQ=255)
                                                                                                                                                 cTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGaa      >  2:2376520/1‑150 (MQ=255)
                                                                                                                                                 cTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGaa      >  2:172524/1‑150 (MQ=255)
                                                                                                                                                 cTACACACGATACTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCAAAGCATTTCACTCAGGCAAGTGACCTGGATCAATTGACACCAGCCATAACACGTATCATACAACGAACCCAACGCGACTAAAGGTCGCAGTGCGaa      >  2:2142440/1‑150 (MQ=38)
                                                                                                                                                  tACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAAc     >  1:982780/1‑150 (MQ=255)
                                                                                                                                                   acacAGGATCGTTTTTGGAGCGGTTTTACCCGCGATGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAAcc    <  2:3833094/150‑1 (MQ=255)
                                                                                                                                                   acacAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAAcc    >  1:1986364/1‑150 (MQ=255)
                                                                                                                                                     acaGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCa  >  1:2100259/1‑150 (MQ=255)
                                                                                                                                                     acaGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCa  >  2:2988431/1‑150 (MQ=255)
                                                                                                                                                     acaGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCa  >  1:3793464/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
AGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCA  >  NC_002947/4844714‑4845011

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