Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,847,875 A→T 5.2% intergenic (+25/‑272) 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,847,8750AT 5.2% 416.6 / 4.9 131intergenic (+25/‑272)ccoP‑II/PP_4259cbb3‑type cytochrome c oxidase subunit/iron‑sulfur cluster‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (34/89);  new base T (3/4);  total (38/93)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 135 total aligned reads displayed.

CGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGT  >  NC_002947/4847728‑4848023
                                                                                                                                                   |                                                                                                                                                     
cGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTaaa                                                                                                                                                     >  2:2578865/1‑150 (MQ=255)
cGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTaaa                                                                                                                                                     <  1:953585/150‑1 (MQ=255)
cGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTaaa                                                                                                                                                     <  1:953950/150‑1 (MQ=255)
 gTTATGGCCTCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGAAACGACAAGGTGCACTTGCTGGCGTCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAAc                                                                                                                                                    <  2:1467007/150‑1 (MQ=255)
    aTGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCg                                                                                                                                                 <  1:4795307/150‑1 (MQ=255)
     tGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGc                                                                                                                                                <  1:3299378/150‑1 (MQ=255)
       gccgccAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTc                                                                                                                                              <  1:5158691/150‑1 (MQ=255)
        ccgccAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCCACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTTTACCCGcgaa                                                                                                                                             <  2:5330984/150‑4 (MQ=255)
        ccgccAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTTTACCCGcgaa                                                                                                                                             <  2:1662329/150‑4 (MQ=255)
        ccgccAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTTTACCCGcgaa                                                                                                                                             <  2:5330933/150‑4 (MQ=255)
        ccgccAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTcc                                                                                                                                             <  2:4199086/150‑1 (MQ=255)
          gccAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTTTACCCGcgaaga                                                                                                                                           <  2:4462167/150‑6 (MQ=255)
             aGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTacac                                                                                                                                        <  2:5348397/150‑1 (MQ=255)
             aGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTacac                                                                                                                                        <  2:2524800/150‑1 (MQ=255)
              ggggCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTAcaca                                                                                                                                       <  1:3107545/150‑1 (MQ=255)
              ggggCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTAcaca                                                                                                                                       <  2:5610998/150‑1 (MQ=255)
              ggggCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTAcaca                                                                                                                                       <  2:2970069/150‑1 (MQ=255)
              ggggCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCCCTCCTAaaca                                                                                                                                       >  1:3490095/1‑150 (MQ=255)
               gggCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAg                                                                                                                                      <  1:183363/150‑1 (MQ=255)
                 gCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGAt                                                                                                                                   <  1:1899957/150‑1 (MQ=255)
                  cAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATc                                                                                                                                  <  1:4001130/150‑1 (MQ=255)
                   aGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATc                                                                                                                                  <  1:3281327/150‑1 (MQ=255)
                    gATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATcc                                                                                                                                 <  2:3847041/150‑1 (MQ=255)
                    gATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCt                                                                                                                                <  2:1142794/150‑1 (MQ=255)
                    gATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCt                                                                                                                                <  2:678928/150‑1 (MQ=255)
                      tGCCGGCCCAGGAACAACTGCAGGGCACCGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACGGATCCtgt                                                                                                                              <  2:328504/150‑3 (MQ=255)
                      tGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCtt                                                                                                                              <  1:2576093/150‑1 (MQ=255)
                      tGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCtt                                                                                                                              <  1:2245550/150‑1 (MQ=255)
                      tGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCtgt                                                                                                                              <  2:795331/150‑3 (MQ=255)
                       gCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCtgtg                                                                                                                             <  2:1528566/150‑4 (MQ=255)
                         cGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTa                                                                                                                           <  1:2111081/150‑1 (MQ=255)
                         cGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTa                                                                                                                           <  2:135715/150‑1 (MQ=255)
                          ggCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAg                                                                                                                          <  2:174681/150‑1 (MQ=255)
                          ggCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAg                                                                                                                          <  2:4432874/150‑1 (MQ=255)
                            cccAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGa                                                                                                                        <  1:5311322/150‑1 (MQ=255)
                            cccAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGa                                                                                                                        <  2:1108091/150‑1 (MQ=255)
                            cccAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGa                                                                                                                        <  1:4139720/150‑1 (MQ=255)
                            cccAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGa                                                                                                                        <  2:4971985/150‑1 (MQ=255)
                             ccAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTacac                                                                                                                                        >  1:1645346/1‑134 (MQ=255)
                             ccAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTacac                                                                                                                                        <  2:1645346/134‑1 (MQ=255)
                             ccAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAg                                                                                                                       <  2:3307626/150‑1 (MQ=255)
                             ccAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCGACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTCCACAGGATCCTGGTAGGAg                                                                                                                       <  2:3307523/150‑1 (MQ=255)
                              cAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGc                                                                                                                      <  2:3964495/150‑1 (MQ=255)
                              cAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGc                                                                                                                      <  1:287926/150‑1 (MQ=255)
                              cAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCCCTCCAACACAGGATCCGTGTAGGAGc                                                                                                                      >  1:695918/1‑150 (MQ=255)
                                ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAGGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCgg                                                                                                                    <  2:1922047/150‑1 (MQ=255)
                                ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCgg                                                                                                                    <  1:3958773/150‑1 (MQ=255)
                                ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCgg                                                                                                                    <  1:3837343/150‑1 (MQ=255)
                                ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCgg                                                                                                                    <  1:430886/150‑1 (MQ=255)
                                ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCgg                                                                                                                    <  1:2691616/150‑1 (MQ=255)
                                ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCgg                                                                                                                    <  1:2692046/150‑1 (MQ=255)
                                ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCgg                                                                                                                    <  1:3958686/150‑1 (MQ=255)
                                ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCgg                                                                                                                    <  2:1092295/150‑1 (MQ=255)
                                ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCgg                                                                                                                    <  2:1389885/150‑1 (MQ=255)
                                ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCgg                                                                                                                    <  2:1389942/150‑1 (MQ=255)
                                ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCgg                                                                                                                    <  2:1390454/150‑1 (MQ=255)
                                 gAACAACTGCAGGGCAACGACACGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCggg                                                                                                                   <  2:593194/150‑1 (MQ=255)
                                 gAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCggg                                                                                                                   <  2:593166/150‑1 (MQ=255)
                                 gAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCggg                                                                                                                   <  1:3122076/150‑1 (MQ=255)
                                 gAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCggg                                                                                                                   <  2:1903336/150‑1 (MQ=255)
                                  aacaacTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCGGGt                                                                                                                  <  1:4295266/150‑1 (MQ=255)
                                  aacaacTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCGGGt                                                                                                                  <  1:2160470/150‑1 (MQ=255)
                                  aacaacGGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCGTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTACGGCAGGCCTGGCCTCCTCGCGGGTACACCCGCTCCTACCCAGGATCCTGGTAGGAGCGGGt                                                                                                                  <  2:1238402/150‑1 (MQ=255)
                                     aacTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCGGGTTTa                                                                                                               <  1:1084043/150‑1 (MQ=255)
                                      acTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCGGGATTAc                                                                                                              <  1:4063897/150‑1 (MQ=255)
                                           aGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTAGCGGGGTAAAACGCGCCTACACAGGTTGCGTGTTTGAGGGGGTTTAACcgcg                                                                                                         >  1:605801/1‑150 (MQ=255)
                                               cAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAaga                                                                                                     <  2:4159959/150‑1 (MQ=255)
                                                aaCGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCgg                                                                                                                    >  2:4488037/1‑134 (MQ=255)
                                                aaCGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAagag                                                                                                    >  2:4137522/1‑150 (MQ=255)
                                                aaCGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTGGTAGGAGCGGGTTTACCCGCGAagag                                                                                                    >  2:3171635/1‑150 (MQ=255)
                                                    aCAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAg                                                                                                                                      <  1:4552865/113‑1 (MQ=255)
                                                          tGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTg                                                                                          >  2:3330221/1‑150 (MQ=255)
                                                          tGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTCGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGATTACCCGCGAAGAGGCCGGTACTg                                                                                          >  2:4587137/1‑150 (MQ=255)
                                                                 gCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACTGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGc                                                                                   >  2:3712643/1‑150 (MQ=255)
                                                                  gtGGCGGCTTACGTCTACCGCCTGTCGCACCTGGCTGAGCCGGCCAAGGCTGCGTACGGCCGGCCTGGCATCTGCGCGGGTAAACCCGCTCCTACACCGGATCCTGGTAGGTGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCa                                                                                  <  2:5015574/149‑1 (MQ=37)
                                                                   tGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCAt                                                                                 >  2:2505075/1‑150 (MQ=255)
                                                                   tGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAAACCGCTCCGACACAGGATACGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCAt                                                                                 >  2:2504864/1‑150 (MQ=255)
                                                                                cTACTGCGTGTCGCACCCGGCTGAGCCGGCCACGGCTGTGTACGGCCGGCCGGGCCTCTTCGGGGGAAACCCCGCTCCTACACTGTATCCTTGTCGGCGGGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCaa                                                                    <  1:4137522/150‑1 (MQ=25)
                                                                                cTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGTGTAAGGCTGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGTATCAGTGTCGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCaa                                                                    <  1:5133490/150‑1 (MQ=38)
                                                                                  tcAGCCTGGTGCGACAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTTAAACCGCTCCTAAACAGGATGCTTGTTGGTGCGGGTTTTACCGCGAAGAGGGCGGTACTGGGACAGCCTTTATAGCAGGGAAGt                                                                  >  1:1257245/2‑150 (MQ=255)
                                                                                     gCCTGTCGCACCCGGCTGAGCCGGCAACGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGTTAAACCGGCTCCTACACCGTATCCTGGTAGGAGGGGTTTTACCGGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGAc                                                               <  2:2050271/150‑1 (MQ=255)
                                                                                              aCCAGGCTGAGCCGGCCAAGGCTGAGTAATGAAGGCCTGGCCTCTTCGCGGGTAAAAACGCTACGACAAATGATCCTTGTTTTAGCGGGTTTACCAGCGAATAGGCCCGTAATGGCCCAGCAATTATAGCAGGCAAGTGACCTGGATCaa                                                      >  1:4551862/1‑150 (MQ=255)
                                                                                              aCCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCaa                                                      >  1:1847085/1‑150 (MQ=255)
                                                                                              aCCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCAACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCaa                                                      >  2:3963758/1‑150 (MQ=255)
                                                                                                   gCTGAGCCGGCAAAGGCTGAGTACGGCCGGCCGGGCCTCTTCGCGGTAAACCCCGCTCCTACACCGGAGCCTGGTCGGAGGGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGac                                                 <  2:1975585/150‑1 (MQ=37)
                                                                                                             ccAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATa                                       <  2:4622117/150‑1 (MQ=255)
                                                                                                                aGGCTGTGTAAGGCTGGCCTGGCCTCTTGGCGGGTAAACCCGCTCATACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAaca                                    <  2:2960644/150‑1 (MQ=255)
                                                                                                                aGGCTGAGTACGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGGAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAaca                                    <  2:5024961/150‑1 (MQ=255)
                                                                                                                            gcaggcCTGGCCTCTTCGCGGGTAAACCCGCTCCTACCCAGGATCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCa                        <  2:4591701/150‑1 (MQ=255)
                                                                                                                              aggcCTGGCCTCTTCGCGGGTAAAACCGCTCCGACCCACGATCCTTGGAGGAGCGGGTTTACCCGCGAAGGGGCCGGTGCTGGGCCACCAGTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCCTAACACGTTTCATAACACg                      >  2:2597200/1‑150 (MQ=25)
                                                                                                                              aggcCTGGCCTCTTCGCGGGTAAAACCGCGCCGACAAAAGATCCGTGGTGGAGCGGGGTTACCCGCGCAGAGGCCGGTACGGGCACCCCCGTTCTCGCCGGCCAGTGTACTGGATCAATTGACACCCGCCCGAACACGGTTCATACCAcc                      >  2:2596869/1‑149 (MQ=19)
                                                                                                                                    gccctCTTCTCGGGTAAACCCGCTCCCACCCCGGTTCCTTTTAGGAGCTGTTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGGAGGCTAGTGACCTGGATGAATTGACAGCCGCCATAACACGTTTCATACCACGAACCCa                <  2:123031/148‑1 (MQ=255)
                                                                                                                                           gagcgGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCCCAGCATTTCTCGCAGGCCAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGAc         >  2:2551970/3‑150 (MQ=255)
                                                                                                                                           gagcgGGTAAACCAGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGAAACCCGCCATAACACGTTTCATAACACGAACCCAACGCGAc         >  2:2552168/3‑150 (MQ=255)
                                                                                                                                             gcgGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTa       >  1:5551538/1‑150 (MQ=255)
                                                                                                                                              cgGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTaa      >  2:2907053/1‑150 (MQ=255)
                                                                                                                                                ggTAACCCCGCTCCTACACCGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCACCGCGACTAAAg    <  2:360862/150‑1 (MQ=255)
                                                                                                                                                 ttAAACCCGCTCCTACACTGTTTCGTTGTTGGTGGGGTTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAgg   <  1:4450685/149‑1 (MQ=255)
                                                                                                                                                 gTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAgg   >  1:367897/1‑150 (MQ=255)
                                                                                                                                                  tAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGt  >  2:3198559/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGT  >  NC_002947/4847728‑4848023

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: