Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,953,069 G→T 8.3% 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,953,0690GT 8.3% 243.4 / 11.9 85intergenic (‑67/+128)PP_3486/PP_3487cytochrome c/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (34/44);  new base T (2/5);  total (36/49)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.60e-01

GAGACGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGC  >  NC_002947/3952920‑3953215
                                                                                                                                                     |                                                                                                                                                  
gagaCGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTgcg                                                                                                                                                    <  2:3715024/150‑1 (MQ=255)
gagaCGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTgcg                                                                                                                                                    <  2:5319239/150‑1 (MQ=255)
  gaCGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCgca                                                                                                                                                  <  1:712777/150‑1 (MQ=255)
  gaCGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCgca                                                                                                                                                  <  1:463780/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:1608804/1‑150 (MQ=255)
                            cTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAgc                                                                                                                        >  2:1608791/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:1215681/150‑1 (MQ=255)
                                        ggtcggtcAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAggg                                                                                                            <  2:82037/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                                                                                         <  1:3122359/150‑1 (MQ=255)
                                                           cagggcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCttt                                                                                         <  2:2784159/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                                                                                                                       <  2:3151110/93‑1 (MQ=255)
                                                                                      aaGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAgcg                                                                                                                       >  1:3151110/1‑93 (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                                                                                                                         >  1:1504438/1‑75 (MQ=255)
                                                                                                      cATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAg                                                                                                                         <  2:1504438/75‑1 (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)
                                                                                                                                                     |                                                                                                                                                  
GAGACGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGC  >  NC_002947/3952920‑3953215

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: