Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,956,096 G→T 8.1% intergenic (‑67/+128) PP_3486 ← / ← PP_3487 cytochrome c/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,956,0960GT 8.1% 246.2 / 11.6 87intergenic (‑67/+128)PP_3486/PP_3487cytochrome c/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (31/49);  new base T (3/4);  total (34/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.12e-01

GAGACGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGG  >  NC_002947/3955947‑3956245
                                                                                                                                                     |                                                                                                                                                     
gagaCGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCGTGCGGCCCGGACGGGCAgcg                                                                                                                                                       >  2:2935964/1‑150 (MQ=255)
 agaCGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCGGTGGGGCGCGGAAGGGCTgggc                                                                                                                                                      >  1:5093818/1‑150 (MQ=255)
 agaCGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTgcgc                                                                                                                                                      <  1:2579410/150‑1 (MQ=255)
  gaCGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCgca                                                                                                                                                     <  1:1185643/150‑1 (MQ=255)
  gaCGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCgca                                                                                                                                                     <  2:1935435/150‑1 (MQ=255)
  gaCGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCgca                                                                                                                                                     <  1:4015663/150‑1 (MQ=255)
     gCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCgcagca                                                                                                                                                  <  2:2921162/150‑1 (MQ=255)
       ccTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCagcagc                                                                                                                                                <  2:2928542/150‑1 (MQ=255)
       ccTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCGTGCGGCGCGGAAGGGCAGCGCagcagc                                                                                                                                                >  2:1777993/1‑150 (MQ=255)
         tCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGccc                                                                                                                                              <  2:3626669/150‑1 (MQ=255)
          cgcgGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGACAGGGCTGCGCAGCAGcccc                                                                                                                                             <  2:2846536/150‑1 (MQ=255)
          cgcgGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGcccc                                                                                                                                             <  1:5090257/150‑1 (MQ=255)
            cgGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCTGCGCAGcccttt                                                                                                                                           <  1:2172670/150‑4 (MQ=255)
            cgGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCTGCGCAGcccttt                                                                                                                                           <  1:2172717/150‑4 (MQ=255)
            cgGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCgg                                                                                                                                           <  2:5529956/150‑1 (MQ=255)
            cgGTAAAGGCGTTTGCCGGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCgg                                                                                                                                           <  2:3026212/150‑1 (MQ=255)
             ggTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCTGCGCAGccctttc                                                                                                                                          <  2:1346496/150‑5 (MQ=255)
             ggTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGc                                                                                                                                          <  2:3777268/150‑1 (MQ=255)
              gTAAAGGCGTTTGCCTGCTTGTTCGTGGTTGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCa                                                                                                                                         <  1:718499/150‑1 (MQ=255)
              gTAAAGGCGTTTGCCTGCTTGTTCGTGGTTGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCa                                                                                                                                         <  1:718408/150‑1 (MQ=255)
              gTAAAGGCGTTTGCCTGCTTGTTCGTGGTTGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCa                                                                                                                                         <  1:1367146/150‑1 (MQ=255)
                 acggCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTg                                                                                                                                      <  2:2408440/148‑1 (MQ=255)
                 aaGGCGTGTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGACGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCTGCAGCCCCGGCAGTg                                                                                                                                      <  2:1664560/150‑1 (MQ=255)
                    gCGTTTGCCTGCTTGTTCGTGGTCGGTCTGGTCCTGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCGCAGCAGCCCCGGCAGTGGCa                                                                                                                                   <  1:2115000/150‑1 (MQ=255)
                    gCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCa                                                                                                                                   <  2:2895454/150‑1 (MQ=255)
                     cGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAg                                                                                                                                  <  1:5494948/150‑1 (MQ=255)
                      gTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCGGGGGTCGCGGCAGGGGGGCGCCGCCCCCCCCGCCGGGGCCGc                                                                                                                                 >  2:838754/1‑150 (MQ=255)
                          gcgtggttgttgGTGGTGGGTCTGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATCTCGCAAGCCCTGTGGGTGCGGCCTTGCTTCGCGTAAGGGCTGCGCTGCCGCCCCGGCAGTGGCAGCTCaa                                                                                                                             <  1:3626149/144‑1 (MQ=255)
                           ccTGCTTGTTCGTGGTCGGTCAGGTCCTGGGCCAGGGCTGAGGTCTAATTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAg                                                                                                                            <  1:806229/150‑1 (MQ=255)
                           ccTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAg                                                                                                                            <  1:1948244/150‑1 (MQ=255)
                             tGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGGCGTGCGTGGCGGAAGGGGTGGGCCCCAGCCCCCGCCAGGGCAGCGCAAgcg                                                                                                                          >  2:1157211/1‑150 (MQ=255)
                             tGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAgcg                                                                                                                          <  2:4484728/150‑1 (MQ=255)
                             tGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGATGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGAAGCTCAAgcg                                                                                                                          >  1:2944847/1‑150 (MQ=255)
                               cTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCa                                                                                                                        <  2:3716647/150‑1 (MQ=255)
                                 tGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGa                                                                                                                      <  1:215598/150‑1 (MQ=255)
                                  gTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCGTGGGTCGCGCGAGGGGTGCGCAGCAGCCCCGGCAGGGGCAGGGCAAgccccgag                                                                                                                     >  2:4199415/1‑144 (MQ=255)
                                   gtCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGaaa                                                                                                                    <  2:87291/149‑1 (MQ=255)
                                        ggtcggtcAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGGCGTGCGGGGCGGAAGGGCTGCGCACCACCCCCCGCCGGGGCCGCACAAACCCCGAAAGAggg                                                                                                               >  2:713070/1‑150 (MQ=255)
                                        ggtcggtcAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAggg                                                                                                               <  2:2944847/150‑1 (MQ=255)
                                        ggtcggtcAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAggg                                                                                                               <  2:2640027/150‑1 (MQ=255)
                                        ggtcggtcAGGTCCAGGGCCAGGGCGGAGGTCGAAGTTGCCATGCCAAGCACGCAGGCAAGCCAGATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAggg                                                                                                               <  2:2640086/150‑1 (MQ=255)
                                          tcggtcAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGa                                                                                                             <  2:2291871/150‑1 (MQ=255)
                                                 ggTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTgc                                                                                                      <  2:2640462/150‑1 (MQ=255)
                                                 ggTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCGGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTgc                                                                                                      <  2:5221321/150‑1 (MQ=255)
                                                    ccagggccagggCTGAGGTCGAAGTTGCCATGCCTAGCACGCTGGTAAGCCATCTCGCAAGCCCTGGGGGAGCGGCCGTGCGTCGCGAATGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTgcgcg                                                                                                   <  1:105084/150‑1 (MQ=39)
                                                     cagggccagggcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGGCGTTCGGCGCCGAAGGGCTGCGCCCCACCCCCCGCCCGGGCCGCACAAGAGCCGCAAGAGGGGAGCGCAgcgcgc                                                                                                  >  2:813228/1‑150 (MQ=255)
                                                     cagggccagggcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTgcgcgc                                                                                                  <  1:150759/150‑1 (MQ=255)
                                                      agggccagggcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCt                                                                                                 <  2:2957285/150‑1 (MQ=255)
                                                      agggccagggcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCt                                                                                                 <  2:2233120/150‑1 (MQ=255)
                                                      agggccagggcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCAGGCGGCGCGGGAGGGCGGCGCACCAGCCCCCGCCGGGGCAGCACAACCCCCGCAAGAGGGGGGGGGGgcgcgcg                                                                                                 >  2:1352777/1‑149 (MQ=255)
                                                       gggccagggcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCGAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTc                                                                                                <  1:5116198/150‑1 (MQ=255)
                                                       gggccagggcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTc                                                                                                <  1:776900/150‑1 (MQ=255)
                                                         gccagggcTGAGGTCGCAGTTGCCAGGCCACGCACGCAGGTAAGCCAGAGCGCCCGCCCGGGGGGCGCGGCCTGGCGTCGCGAAAGGGCGGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCGCAGCGCTccc                                                                                              <  2:3396716/150‑2 (MQ=255)
                                                          ccagggcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGAGGAGCGCGCAGCGCTcccc                                                                                             <  1:2734242/150‑3 (MQ=255)
                                                            agggcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAACGGGCGGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTc                                                                                           <  1:3883477/150‑1 (MQ=255)
                                                            agggcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTc                                                                                           <  1:2128758/150‑1 (MQ=255)
                                                             gggcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTcg                                                                                          <  1:2949819/150‑1 (MQ=255)
                                                              ggcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTcgc                                                                                         <  2:2497643/150‑1 (MQ=255)
                                                              ggcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTcgc                                                                                         <  2:1076290/150‑1 (MQ=255)
                                                              ggcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTcgc                                                                                         <  1:3063823/150‑1 (MQ=255)
                                                                 tGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGAc                                                                                      <  1:363619/150‑1 (MQ=255)
                                                                 tGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGAc                                                                                      <  1:3178791/150‑1 (MQ=255)
                                                                   aGGTCGAAGTTGCCATGCCAAGCACGCTGGTAAGCCAGATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCTGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGc                                                                                    <  2:3376143/150‑1 (MQ=255)
                                                                     gTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCGTGCGGCGCGAAAGGGCGGCGCCCCCGCCCCCGCAGGTGCCGCACAAGCGCAGAAAGAGGGGAGCGCTGCGCGCGCCGTTCCCGCCGCaa                                                                                  >  2:1801499/1‑150 (MQ=255)
                                                                     gTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCGGGCGGCGCGGAAGGGGTGCGCCGCAGCCCCCGCCGGGGGCGCGCCAGCGCCGAAAGAGGGGTGCGCTGCGCGCGCCGATCGCGACGCac                                                                                  >  2:1325542/1‑149 (MQ=255)
                                                                       cGACGTTGCCATGCCAAGCACGCAGGTACGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAgg                                                                                <  1:1011943/150‑1 (MQ=255)
                                                                       cGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGGCGTGCGGCGCGCGCGGGGTGGGCAGCCGCCCCGGCCCGGGCAGCGCCAGCGCCGCAAGAGGGGGGCGCGGCGCGCCCCCTTCGCGGCGCACgg                                                                                >  2:4930136/1‑150 (MQ=255)
                                                                       cGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAgg                                                                                <  2:1728513/150‑1 (MQ=255)
                                                                           gTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCGTTCGTGGCGAAAGGGCTGCGTATCAGCCCCCGCAGTGGCCGCGCCACCGCACCAAGAGGGGGGCGCAGCCCCCACCGTTAGCCGCGCAACGcccc                                                                            >  2:4726628/1‑148 (MQ=255)
                                                                              gccatgccaAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCGTGCGGCGCGCCACGGCGGCGCCGCCCCCCCCGCCCGGGCCGCGCCCGCGCCCCAAGGGGGGCGCGCGGCGCGCGCCCGTCGCGGCGCACGGCCGCGcc                                                                         >  2:4147209/1‑150 (MQ=255)
                                                                               ccatgccaAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTAGCGGAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCGCAAGCCCAGAAAGAGGGGGGCGCGGCGCGCGCCTTTCGCGACGCAGGGCCGCGccc                                                                        >  2:4148813/1‑150 (MQ=38)
                                                                                      aaGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGGCGTGCGGCGCGGAAGGGGGGCGCAGCAGCCCCCGCAGGGGGCGCGCAACCGCAGGAAGAGGGGGGCGCTGCGCGCGCCCTTGGCGCCGCCAGGCCCCGccccccacac                                                                 >  2:4485313/1‑143 (MQ=255)
                                                                                      aaGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTTGGTGGGGCAAGGGGTGCGCAGCAGCCCCCGCCCTGGCCGCTCCACCGCAGCAAGGGGGGGGCGCGGCGCGCGCCGTTCGCGCCGCCAGGGCGCGcccccccacg                                                                 >  2:4485372/1‑143 (MQ=255)
                                                                                          aCGCTGTTACGCAATATCGCAAGCCCTGGGGGCGCGGCCTTGCGTCGCCAATGGGGTGCGCAGCAGCCCCGGCCGTGTCAGCTCACGCGCATAATGTGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGa                                                             <  2:4148306/150‑1 (MQ=25)
                                                                                                       atatCGCAAGCCCTGTGGGAGCGGCCTTGCGGCGCGAAAGGGCTGCGCCGCCCCCCCGGCACTGGCCGCTCAAGCGCAGAACGAGGGGGGCGCGGCGCGCTCCGTTCGCGAAGCAAGGCCGCTCCCCCAAAAGCGGAAGCCACAGCgggc                                                >  2:3207276/1‑148 (MQ=255)
                                                                                                                              gCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATa                         >  2:3360110/1‑150 (MQ=255)
                                                                                                                               ccTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATAc                        >  1:2717367/1‑150 (MQ=255)
                                                                                                                               ccTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATAc                        >  2:423613/1‑150 (MQ=255)
                                                                                                                               ccTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATAc                        >  2:1479582/1‑150 (MQ=255)
                                                                                                                                cTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACt                       >  2:1432831/1‑150 (MQ=255)
                                                                                                                                     gTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCaccg                  >  1:2388235/1‑150 (MQ=255)
                                                                                                                                      tCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCaccga                 >  1:2154335/1‑150 (MQ=255)
                                                                                                                                        gcgAACGGGCTGCGCAGCACCCCCGGCAGTGGCAGCTCACGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCaccgacc               <  2:5093818/150‑1 (MQ=255)
                                                                                                                                        gcgAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCaccgacc               >  1:59922/1‑150 (MQ=255)
                                                                                                                                          gAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCgg             >  1:2615935/1‑150 (MQ=255)
                                                                                                                                          gAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCgg             >  2:2332632/1‑150 (MQ=255)
                                                                                                                                          gAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCgg             >  1:1328060/1‑150 (MQ=255)
                                                                                                                                            cggggCTGCTGCGCAGCCCTGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGa                           <  1:1079289/132‑1 (MQ=39)
                                                                                                                                            cggggCTGCTGCGCAGCCCTGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGa                           >  2:1079289/3‑134 (MQ=39)
                                                                                                                                             ggggCTGCTGCGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCa          >  2:603775/2‑150 (MQ=255)
                                                                                                                                              gggCTGCTGCGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAg         >  1:4879894/1‑150 (MQ=255)
                                                                                                                                              gggCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAg         >  2:397309/1‑150 (MQ=255)
                                                                                                                                              gggCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAg         >  1:5112437/1‑150 (MQ=255)
                                                                                                                                                  tGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGc     >  1:4173443/1‑150 (MQ=255)
                                                                                                                                                   gcgcAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCt    >  1:2088325/1‑150 (MQ=255)
                                                                                                                                                     gcagccGCCCGGGCTGTGGCTGCTCACGCGCTTAACGTGGGTCGGGGTGCGCGCTCCTTTCGCGCCGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTgg  <  2:2127309/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAGACGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGG  >  NC_002947/3955947‑3956245

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