Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,956,098 A→C 15.1% 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,956,0980AC15.1% 184.0 / 19.3 86intergenic (‑69/+126)PP_3486/PP_3487cytochrome c/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (24/46);  new base C (10/4);  total (34/52)
Fisher's exact test for biased strand distribution p-value = 1.56e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGT  >  NC_002947/3955949‑3956247
                                                                                                                                                     |                                                                                                                                                     
gaCGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCgca                                                                                                                                                       <  2:1935435/150‑1 (MQ=255)
gaCGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCgca                                                                                                                                                       <  1:1185643/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:2172717/150‑4 (MQ=255)
          cgGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCTGCGCAGcccttt                                                                                                                                             <  1:2172670/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:1367146/150‑1 (MQ=255)
            gTAAAGGCGTTTGCCTGCTTGTTCGTGGTTGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCa                                                                                                                                           <  1:718499/150‑1 (MQ=255)
            gTAAAGGCGTTTGCCTGCTTGTTCGTGGTTGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCa                                                                                                                                           <  1:718408/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:2640027/150‑1 (MQ=255)
                                      ggtcggtcAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAggg                                                                                                                 <  2:2944847/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:1076290/150‑1 (MQ=255)
                                                            ggcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTcgc                                                                                           <  2:2497643/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                          >  2:1479582/1‑150 (MQ=255)
                                                                                                                             ccTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATAc                          >  1:2717367/1‑150 (MQ=255)
                                                                                                                             ccTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATAc                          >  2:423613/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                             >  2:1079289/3‑134 (MQ=39)
                                                                                                                                          cggggCTGCTGCGCAGCCCTGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGa                             <  1:1079289/132‑1 (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)
                                                                                                                                                     agcagcCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGt  >  2:3280527/1‑150 (MQ=255)
                                                                                                                                                     agcagcCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGt  >  2:4735557/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GACGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGT  >  NC_002947/3955949‑3956247

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: