Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,953,071 A→C 17.6% intergenic (‑69/+126) 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,953,0710AC17.6% 170.9 / 17.7 85intergenic (‑69/+126)PP_3486/PP_3487cytochrome c/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (26/41);  new base C (8/8);  total (36/49)
Fisher's exact test for biased strand distribution p-value = 5.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGT  >  NC_002947/3952922‑3953220
                                                                                                                                                     |                                                                                                                                                     
gaCGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCgca                                                                                                                                                       <  1:463780/150‑1 (MQ=255)
gaCGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCgca                                                                                                                                                       <  1:712777/150‑1 (MQ=255)
 aCGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCgcag                                                                                                                                                      <  1:4055281/150‑1 (MQ=255)
      cTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCTGCGCAGcc                                                                                                                                                 <  2:742896/150‑1 (MQ=255)
      cTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCTGCGCAGcc                                                                                                                                                 <  2:742801/150‑1 (MQ=255)
        cgcgGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGcccc                                                                                                                                               <  1:575030/150‑1 (MQ=255)
          cgGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCTGCGCAGcccttt                                                                                                                                             <  1:1098841/150‑4 (MQ=255)
           ggTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAGGGGCTGCTGCGCAGccctttc                                                                                                                                            <  2:4186749/150‑5 (MQ=255)
           ggTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCTGCGCAGccctttc                                                                                                                                            <  2:3282117/150‑5 (MQ=255)
             tAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGGCGCGAAAGGGCAGCGCAGCAGCCCCCGCac                                                                                                                                          >  2:4527776/1‑149 (MQ=255)
                aGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTgg                                                                                                                                       <  2:1759454/150‑1 (MQ=255)
                aGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTgg                                                                                                                                       <  2:1759397/150‑1 (MQ=255)
                     tttGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCt                                                                                                                                  <  1:4226385/150‑1 (MQ=255)
                         ccTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAg                                                                                                                              >  1:1215681/1‑150 (MQ=255)
                          cTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGGCTTGCCTCGCGCAAGGGCAGCGCAGCAGCCCCCGCACAGGCACCGCAAgc                                                                                                                             >  2:2211224/1‑150 (MQ=255)
                          cTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAgc                                                                                                                             >  2:1608791/1‑150 (MQ=255)
                          cTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAgc                                                                                                                             >  2:1608804/1‑150 (MQ=255)
                              ttGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCATGCGTCGCGGAACGGCTTCGCCCCACCCCCCCCCCTGGCAGCACAAGCGCac                                                                                                                         >  2:739116/1‑149 (MQ=255)
                               tGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGa                                                                                                                        <  2:3097343/150‑1 (MQ=255)
                               tGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGa                                                                                                                        >  2:2763267/1‑150 (MQ=255)
                                    gTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGGCGTGCGGCGCGGCAGGGGTGCGCCGCCGCCCCCGCCGTGGCCGCGCAAGCGCAGAAagag                                                                                                                   >  2:3560719/1‑150 (MQ=255)
                                    gTGGGCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTGGCCATGCCATGCACGCAGGTAACCCATAGCGCAAGCCCGGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAagag                                                                                                                   <  1:5381497/150‑1 (MQ=39)
                                     tGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGGGGGGCGGGGCGGGAGGGGGGCGCGGCCGCCCCGGCCGGGGCGGCTCGAGCGCCGaccgggg                                                                                                                  >  1:3171754/1‑144 (MQ=255)
                                     tGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAgg                                                                                                                  <  2:3342856/150‑1 (MQ=255)
                                     tGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAgg                                                                                                                  <  2:3931873/150‑1 (MQ=255)
                                      ggtcggtcAGGTCCTGGGCCTGGGCTGTGGTCGAAGTTGCCTTGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAggg                                                                                                                 <  1:4164196/150‑1 (MQ=255)
                                      ggtcggtcAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCGGTGCGGCGCGGGGGGGCTGCGCCGCCGCCCCCGCAGGGGGCGCGCAAGCGCCCGAAGAggg                                                                                                                 >  2:2656809/1‑150 (MQ=255)
                                      ggtcggtcAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAggg                                                                                                                 <  2:82037/150‑1 (MQ=255)
                                      ggtcggtcAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAggg                                                                                                                 <  2:1215681/150‑1 (MQ=255)
                                       gtcggtcAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAgggg                                                                                                                <  2:663940/150‑1 (MQ=255)
                                        tcggtcAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGa                                                                                                               <  1:2763267/150‑1 (MQ=255)
                                            tcAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCgcag                                                                                                                                                      <  2:2699013/107‑1 (MQ=255)
                                            tcAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCgcag                                                                                                                                                      >  1:2699013/1‑107 (MQ=255)
                                            tcAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCGTGCGGGGCGCAAGGGCGGCGCGCCCGCCCCCGCCGGGGCAGCGCCAGCCCAGACAGGGGGgcgcgc                                                                                                           >  2:2575899/1‑150 (MQ=255)
                                                  ccagggccagggCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTgcgcg                                                                                                     <  1:2915702/150‑1 (MQ=255)
                                                   ctgggccagggcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTgcgcgc                                                                                                    <  1:1490874/148‑1 (MQ=255)
                                                   cagggcaagggcTGAGTTCTAAGTTGCAATGCAAAGCACGCAGGAAAGCAATATCGAAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTgcgcgc                                                                                                    <  1:3390898/150‑1 (MQ=38)
                                                     gggccagggcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTc                                                                                                  <  1:1608791/150‑1 (MQ=255)
                                                     gggccagggcTGAGGTCGAAGTTGCCATGCAAAGCACGCAGGTAAGCAATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTc                                                                                                  <  1:1608804/150‑1 (MQ=255)
                                                        ccagggcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCgcag                                                                                                                                                      <  1:3441489/95‑1 (MQ=255)
                                                        ccagggcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCgcag                                                                                                                                                      >  2:3441489/1‑95 (MQ=255)
                                                        ccagggcGGAGGTCGAAGTGGCCATGCCCAGCACGCAGGTAAGCCATCTCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCGCAGCGCTcccc                                                                                               <  2:4286412/150‑3 (MQ=38)
                                                         cagggcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCttt                                                                                              <  2:2784159/150‑1 (MQ=255)
                                                         cagggcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCttt                                                                                              <  1:3122359/150‑1 (MQ=255)
                                                          agggcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCCTGCGGCGCGCAAGGGCGGCGCCCCCGCCCCCGCCCGGGCCCCGCAAGCGCCGCACGCGGGGGGCGCGGCGCGCGcccggc                                                                                             >  2:1773851/1‑146 (MQ=255)
                                                          agggcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGGGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTc                                                                                             <  1:1042135/150‑1 (MQ=255)
                                                           gggcTGTGGTCGAAGTTGCCATGCCAAGCACGCTGGTAAGCCTTATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAACGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTcg                                                                                            <  1:5008714/150‑1 (MQ=255)
                                                           gggcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTcg                                                                                            <  2:3181465/150‑1 (MQ=255)
                                                            ggcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTcgc                                                                                           <  1:830450/150‑1 (MQ=255)
                                                               tGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTCTCGCGAc                                                                                        <  1:709394/150‑1 (MQ=255)
                                                                 aGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGGCTTGCGGCGCGGAAGGGGGGCGCCGCAGCCCCCGCAGTGGCAGCAAAAGCGCAGGAACAGGGGCGCGCGGCGCGCGGCGTTCGCGgcgc                                                                                      >  2:1743706/1‑150 (MQ=255)
                                                                 aGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGc                                                                                      <  2:291176/150‑1 (MQ=255)
                                                                   gTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCaa                                                                                    <  2:5383296/150‑1 (MQ=255)
                                                                      gAAGTTGCCTTGCCAAGCACGCAGTTAAGCAATATCGCAAGCCCGGTGGTAGCGGCCTTGCGTCGCGAACGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGc                                                                                 <  1:2467879/150‑1 (MQ=255)
                                                                       ttgTTGCCTTGCCAAGCACGCTGGTATGACATATCGCAAGCCCTGTGGGTGCGGCCTTGCGTCGCGACTGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCATAAAGAGGGGCGCGCTGCGCGCTCCTTTCGCGACGCAAGGcc                                                                                <  1:756236/148‑1 (MQ=38)
                                                                       aaGTTGCCTTGCCAAGCCCGCAGGTAAGCCATATCGCAAGCCCGGGGGGAGCGGCCTTGCGTCGCGAACGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGcc                                                                                <  1:756386/150‑1 (MQ=39)
                                                                              catgccaAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGGCGTTGGGGGCGGAAGGGCTGCGCAGCAGCCCCCGCAGTGGCAGCGCCAGCGCCGAAAGAGGGGGGCGCTGCGCGCTCCGGTCGCGACGCCAGGGCGCGcccc                                                                         >  2:824023/1‑149 (MQ=255)
                                                                                 gccaAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGGCGTGCGGGGCGCAAAGGCGGCGCCGCAGCCCCCGCCCGGGGCGCGCAAGCCCCCGAACAGGGGGGCGGGGCGCGCACCCGTGGCGGCCCCACGCCGCGCCcacac                                                                      >  2:4318866/1‑149 (MQ=255)
                                                                                    aaGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAgcg                                                                                                                            >  1:3151110/1‑93 (MQ=255)
                                                                                    aaGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAgcg                                                                                                                            <  2:3151110/93‑1 (MQ=255)
                                                                                    aaGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTc                                                                                             >  1:1121236/1‑124 (MQ=255)
                                                                                    aaGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTc                                                                                             <  2:1121236/124‑1 (MQ=255)
                                                                                      gCACGCAGGTACGCCTTATCGCACGCCGTGTGGGAGCGGCTTGGGGTCGCGAACGGGCTGCGCAGCAGCCCCGGCTGTGGCAGCACACGCGAATAATGTGGGGTGCGCTGCGCGCTCATTTGGCGCCGCAAGGCCGCTCCCACAAAAGCg                                                                 <  2:4915382/150‑1 (MQ=255)
                                                                                          gCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCGTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCCGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCATTCGCGACTCAAGGCCGCACCCACAAAAGCCGATg                                                             >  2:1653128/1‑150 (MQ=39)
                                                                                            aGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAagag                                                                                                                   >  1:1052395/1‑94 (MQ=255)
                                                                                            aGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAagag                                                                                                                   <  2:1052395/94‑1 (MQ=255)
                                                                                                    cATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAg                                                                                                                              <  2:1504438/75‑1 (MQ=255)
                                                                                                    cATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAg                                                                                                                              >  1:1504438/1‑75 (MQ=255)
                                                                                                          gCAAGCCCTGTGGGAGCGGGCCTTCGTCGCGCAAGGGCGGCGCAGCCGCCCCCGCACGGGCAGCGCCAGCGCCGCCAGCGGGGCGCGCGGCGCCCGCCGTTCGCGGCGCCAGCGCGCGGCCCCCAAAGCGGCGGCCACAGAGGcccaata                                             >  2:3548713/1‑145 (MQ=255)
                                                                                                          gCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCt                                                                                                          <  2:311677/89‑1 (MQ=255)
                                                                                                          gCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCt                                                                                                          >  1:311677/1‑89 (MQ=255)
                                                                                                          gCAAGCCCTGGGGGCGCGGCCTTGCGTCGCTATCGGGCTGCGCAGCAGCCCGGGCTTTGGCTGAACAAGCGCTTATAGTGGGGTGCGCTGCGCTATCATTTCGCGACGCACGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTc                                             <  2:2725016/150‑1 (MQ=255)
                                                                                                                          tgggcTTGCGTCGGGCTCGGGCTGCGCCGCCCCCCGGGCCGGGGCAGATCAAGCGCATAACGTGGGGAGCGCTGCGCGCTCATTGCGCGCCGCACGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGa                             <  2:1173343/146‑1 (MQ=255)
                                                                                                                            gCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATa                           >  1:1173343/1‑150 (MQ=255)
                                                                                                                             ccTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATAc                          >  2:2740739/1‑150 (MQ=255)
                                                                                                                                 gCGTCGCTAAAGGGCTGCGCCGCCGCCCCGGCCGTGGCAGCCCACGCGCATAACGAGGGGAGCGCTGCGCGCTCATTGCGCGACGCATGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCac                      <  2:4951468/150‑1 (MQ=255)
                                                                                                                                    tCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCaccga                   >  1:998975/1‑150 (MQ=255)
                                                                                                                                     cgcgAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCaccgac                  >  1:1423702/1‑150 (MQ=255)
                                                                                                                                     cgcgAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCaccgac                  >  1:178355/1‑150 (MQ=255)
                                                                                                                                        gAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCgg               >  2:1354385/1‑150 (MQ=255)
                                                                                                                                        gAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCgg               >  2:4961031/1‑150 (MQ=255)
                                                                                                                                         ctctttcTGCGCAGCAGCCCCGGCAGTGGCAGCTCACGCGCAGACAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGc              >  2:2595165/7‑150 (MQ=255)
                                                                                                                                         ccggggCTGCTGCGCAGCCCTGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGc              >  2:5409959/4‑150 (MQ=39)
                                                                                                                                           ggggCTGCTGCGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCa            >  1:1309947/2‑150 (MQ=255)
                                                                                                                                           aGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCa            <  1:4961031/150‑1 (MQ=255)
                                                                                                                                             ggCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGc          >  1:4733924/1‑150 (MQ=255)
                                                                                                                                              gCTGCGCCGCCGCCCCGGCTGTGGCTGCTCATGCGCCTAATGGGGGGTGGGCTGGGCGCTCTTTTGGCGACGCACGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCg         <  2:2673248/150‑1 (MQ=255)
                                                                                                                                              gCTGCGCAGCAGCCCCGGCTGTGGCTGCTCACGCGCAGAATGTGGGGAGCGCTGCGCGCTCATTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCg         <  2:5397897/150‑1 (MQ=255)
                                                                                                                                              gCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCg         >  1:2642760/1‑150 (MQ=255)
                                                                                                                                                tGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACagag                                                     >  1:3915808/1‑104 (MQ=255)
                                                                                                                                                tGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACagag                                                     <  2:3915808/104‑1 (MQ=255)
                                                                                                                                                tGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGc       >  2:169594/1‑150 (MQ=255)
                                                                                                                                                     agcagcCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGt  >  2:5074466/1‑150 (MQ=255)
                                                                                                                                                     agcagcCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGt  >  2:505998/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GACGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGT  >  NC_002947/3952922‑3953220

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: