Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,859,463 A→T 6.4% intergenic (‑55/+101) PP_1662 ← / ← PP_1663 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,859,4630AT 6.4% 237.1 / 4.4 78intergenic (‑55/+101)PP_1662/PP_1663hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (32/41);  new base T (2/3);  total (34/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

GTGGGACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCA  >  NC_002947/1859319‑1859610
                                                                                                                                                |                                                                                                                                                   
gTGGGACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGttt                                                                                                                                                 <  5:242325/149‑1 (MQ=255)
   ggACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACTTGTTTCAt                                                                                                                                              <  1:57790/149‑1 (MQ=255)
    gACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATg                                                                                                                                             <  3:241535/149‑1 (MQ=255)
       aTTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATg                                                                                                                                             >  5:204915/1‑146 (MQ=255)
              aTTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAg                                                                                                                                   <  3:126376/149‑1 (MQ=255)
               ttCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAgg                                                                                                                                  <  6:299451/149‑1 (MQ=255)
                      cgTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAggg                                                                                                                                 >  2:318974/1‑143 (MQ=255)
                      cgTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAggg                                                                                                                                 <  1:318974/143‑1 (MQ=255)
                        ttCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGcc                                                                                                                         <  6:127916/149‑1 (MQ=255)
                             ttACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCtgggt                                                                                                                    >  3:318175/1‑149 (MQ=255)
                              tACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCtgggt                                                                                                                    <  4:318175/148‑1 (MQ=255)
                               aCTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCtgggtg                                                                                                                   <  5:187559/148‑1 (MQ=255)
                                cTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCtgggtg                                                                                                                   >  7:138666/1‑147 (MQ=255)
                                cTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCtgggtg                                                                                                                   <  8:138666/147‑1 (MQ=255)
                                 ttGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGc                                                                                                                >  4:269103/1‑149 (MQ=255)
                                       aGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTc                                                                                                          >  7:170479/1‑149 (MQ=255)
                                         tCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCg                                                                                                         >  2:230425/1‑148 (MQ=255)
                                         tCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCg                                                                                                         >  5:48901/1‑148 (MQ=255)
                                           gtCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCgca                                                                                                       <  1:273716/147‑1 (MQ=255)
                                           aTCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCgca                                                                                                       <  8:157325/148‑1 (MQ=255)
                                           aTCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCgca                                                                                                       <  4:141300/148‑1 (MQ=255)
                                                tAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGttt                                                                                                                                                 >  6:351527/1‑101 (MQ=255)
                                                tAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGttt                                                                                                                                                 <  5:351527/101‑1 (MQ=255)
                                                tAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCgcagca                                                                                                    >  6:316992/1‑146 (MQ=255)
                                                tAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAg                                                                                                 <  6:204915/149‑1 (MQ=255)
                                                    cTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGaaaa                                                                                             <  1:230425/149‑1 (MQ=255)
                                                        aCTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGccc                                                                                         >  8:145827/1‑149 (MQ=255)
                                                         cttcttCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGcccc                                                                                        <  3:80174/149‑1 (MQ=255)
                                                          ttcttcAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCg                                                                                       >  2:288978/1‑149 (MQ=255)
                                                             ttcAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGAATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGTTGGGCTTTTTCGAAGCAAAGAAAAGCCCCGCAc                                                                                    >  7:302495/1‑149 (MQ=255)
                                                                aaGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCgca                                                                                                       >  5:37856/1‑127 (MQ=255)
                                                                aaGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCgca                                                                                                       <  6:37856/127‑1 (MQ=255)
                                                                aaGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTg                                                                                  <  6:158870/148‑1 (MQ=255)
                                                                  gATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACTTGTTTCAt                                                                                                                                              <  5:53405/86‑1 (MQ=255)
                                                                  gATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACTTGATTCAt                                                                                                                                              >  6:53405/1‑86 (MQ=255)
                                                                  gATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGc                                                                               >  4:163019/1‑149 (MQ=255)
                                                                    tGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCgg                                                                             <  1:288978/149‑1 (MQ=255)
                                                                    tGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCgg                                                                             <  6:48901/149‑1 (MQ=255)
                                                                    tGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCgg                                                                             <  8:170479/149‑1 (MQ=255)
                                                                    tGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCgg                                                                             <  5:60261/149‑1 (MQ=255)
                                                                     ggCTGCGACAGTGTGTCGCACGTCTTCTTTAGAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAACGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCggg                                                                            >  2:200403/1‑149 (MQ=255)
                                                                       gtGCGACACACTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCggg                                                                            >  1:282640/2‑147 (MQ=255)
                                                                       gtGCGACACACTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCggg                                                                            <  2:282640/146‑1 (MQ=255)
                                                                       cTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCCCAGTGCGGGGc                                                                          <  5:316992/149‑1 (MQ=255)
                                                                        tGCGACAGTGTGTCGCACGTCTTATTTATAGAGTTTTTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACGAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTATTCTCAGCAAAAAAAAACACCGAAATGGGCtggga                                                                          >  2:304391/1‑147 (MQ=255)
                                                                            agaGTGTGTCGCACGTGTTGTTTATAGAGTTTGTTAGCAGACTATGTAAGTAAAAGATGAATGAAACAATTTTAATGAAGCTCGAAGGGATGATGGCTGGGTGGGCTTTTTCGCATCTAAGAAATGCCGCGCACTGGGGGGTGCttaa                                                                      <  8:302495/146‑3 (MQ=255)
                                                                            aCAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTATGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAg                                                                     <  8:75206/149‑1 (MQ=255)
                                                                             cAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCCCAGTGCGGGGCttttct                                                                    <  3:269103/149‑5 (MQ=255)
                                                                                    tCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTgc                                                              <  1:200403/148‑1 (MQ=255)
                                                                                      gCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTgcagg                                                           >  2:88416/1‑149 (MQ=255)
                                                                                      gCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTgcag                                                            <  3:10418/148‑1 (MQ=255)
                                                                                          gTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGt                                                       <  3:163019/149‑1 (MQ=255)
                                                                                            cttcttTATAGAGTTTGTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCAATGGGCGGGGCTTCAGGTATTGCAGGCAGTgcg                                                    >  7:84468/1‑149 (MQ=255)
                                                                                             ttcttTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCgg                                                                             >  5:222324/1‑124 (MQ=255)
                                                                                             ttcttTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCgg                                                                             <  6:222324/124‑1 (MQ=255)
                                                                                             ttcttTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTgca                                                             <  3:335868/140‑1 (MQ=255)
                                                                                             ttcttTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTgca                                                             >  4:335868/1‑140 (MQ=255)
                                                                                                tttATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGa                                                 >  5:119442/1‑149 (MQ=255)
                                                                                                   ataGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTgc                                                              <  6:108607/133‑1 (MQ=255)
                                                                                                   ataGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTgc                                                              >  5:108607/1‑133 (MQ=255)
                                                                                                     agagTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAg                                                                     <  2:317102/124‑1 (MQ=255)
                                                                                                     agagTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAg                                                                     >  1:317102/1‑124 (MQ=255)
                                                                                                      gagTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAACAAAGAAAAGCCCCGCACTGGGCGGGGCTTCa                                                                      >  8:40637/1‑122 (MQ=255)
                                                                                                      gagTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAACAAAGAAAAGCCCCGCACTGGGCGGGGCTTCa                                                                      <  7:40637/122‑1 (MQ=255)
                                                                                                           tGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGc                                         <  4:173538/146‑1 (MQ=255)
                                                                                                           tGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGc                                         >  3:173538/1‑146 (MQ=255)
                                                                                                                gCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCAcc                                 <  1:88416/149‑1 (MQ=255)
                                                                                                                gCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCAc                                  <  6:119442/148‑1 (MQ=255)
                                                                                                                   gACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGc                              >  6:342633/1‑149 (MQ=255)
                                                                                                                     cTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCtt                            >  3:15289/1‑149 (MQ=255)
                                                                                                                                   gcttcATGAAACTTGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGccc                                       >  6:37709/5‑124 (MQ=255)
                                                                                                                                   gcttcATGAAACTTGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGccc                                       <  5:37709/120‑1 (MQ=255)
                                                                                                                                     ggCATGAAACGAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCa            <  1:304391/149‑1 (MQ=255)
                                                                                                                                          gagaCAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTg       >  5:80772/4‑149 (MQ=255)
                                                                                                                                            accAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGGACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTgca     <  4:15289/147‑1 (MQ=255)
                                                                                                                                            aaCAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTg       >  3:305114/1‑147 (MQ=255)
                                                                                                                                              cAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGc                              <  7:299772/122‑1 (MQ=255)
                                                                                                                                              cAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGc                              >  8:299772/1‑122 (MQ=255)
                                                                                                                                               aaGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTgcagca  >  6:85277/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GTGGGACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCA  >  NC_002947/1859319‑1859610

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: