Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 58,747 T→A 10.4% intergenic (‑658/+30) PP_5424 ← / ← PP_0050 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294758,7470TA10.4% 218.3 / 12.4 84intergenic (‑658/+30)PP_5424/PP_0050hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (39/35);  new base A (5/4);  total (44/40)
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.46e-01

TGAGTAGGGATTTCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAGGTGGCGC  >  NC_002947/58602‑58888
                                                                                                                                                   |                                                                                                                                              
tGAGTAGGGATTTCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCgcagca                                                                                                                                             <  2:189532/149‑2 (MQ=255)
  aGTAGGGATTTCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGctgctg                                                                                                                                            <  1:154747/148‑1 (MQ=255)
     aGGGATTTCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCgg                                                                                                                                         <  3:227658/148‑1 (MQ=255)
      gggATTTCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGc                                                                                                                                       >  2:318179/1‑149 (MQ=255)
       gtttttCACGGATTTTATGTAAAATGAACCCATTAATTTAAAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGGCACTAAAAAGCCCGCAGCAGCGGGCtt                                                                                                                                     <  5:192262/146‑1 (MQ=255)
       ggATTTCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCt                                                                                                                                      >  2:224585/1‑149 (MQ=255)
        gATTTCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGGCACTAAAAAGCCCGCAGCAGCGGGCtt                                                                                                                                     <  2:187396/149‑1 (MQ=255)
          tttCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCtttt                                                                                                                                   <  2:271233/149‑1 (MQ=255)
             ctcGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAg                                                                                                                                >  5:345505/1‑149 (MQ=255)
             ctcGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCAGCAGCGGGCTTTTTa                                                                                                                                 <  3:161120/148‑1 (MQ=255)
              tcGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCtttt                                                                                                                                   <  7:278309/145‑1 (MQ=255)
              tcGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTCTAATTCGCCAGACACTAAAAAGACCGCTGCTGCGGCCtttt                                                                                                                                   >  8:278309/1‑145 (MQ=255)
               gtgCTTTTCACAGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTg                                                                                                                              >  1:133426/3‑149 (MQ=255)
                ggCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTa                                                                                                                                 >  4:205412/1‑145 (MQ=255)
                ggCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATCCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTa                                                                                                                                 <  3:205412/145‑1 (MQ=255)
                        acagAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCaaa                                                                                                                      >  5:333877/4‑148 (MQ=255)
                          tgAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACaa                                                                                                                   >  3:76043/1‑149 (MQ=255)
                           gAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAg                                                                                                                                >  6:346371/1‑135 (MQ=255)
                           gAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAg                                                                                                                                <  5:346371/135‑1 (MQ=255)
                            aaaaGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCggg                                                                                                                                        >  5:162535/1‑126 (MQ=255)
                            aaaaGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCggg                                                                                                                                        <  6:162535/126‑1 (MQ=255)
                            aaaaGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAAt                                                                                                                  <  1:318179/148‑1 (MQ=255)
                              aaGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAg                                                                                                               <  7:116935/149‑1 (MQ=255)
                                gCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCTCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTa                                                                                                             <  1:60221/149‑1 (MQ=255)
                                     ccTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCgg                                                                                                        >  8:280738/1‑149 (MQ=255)
                                      cTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTTGTGCCTGCAAACAATCAGTAGCCGGt                                                                                                       <  1:9598/149‑1 (MQ=255)
                                      cTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGt                                                                                                       <  6:345505/149‑1 (MQ=255)
                                            ttCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAAt                                                                                                 >  5:340733/1‑149 (MQ=255)
                                              cacaAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACTA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGt                                                                                               >  6:340842/1‑149 (MQ=255)
                                              cacaAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGt                                                                                               >  5:73752/1‑149 (MQ=255)
                                               acaAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGt                                                                                               >  5:321386/1‑148 (MQ=255)
                                                 aaaaaaaaTGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAAGGGGAGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAAAAATCAGTAGCCGGTGAGAATGt                                                                                               >  1:280068/1‑148 (MQ=255)
                                                 aaaaaaaaTGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCa                                                                                            <  4:56570/149‑1 (MQ=255)
                                                 aaaaaaaaTGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCa                                                                                            >  5:220285/1‑149 (MQ=255)
                                                  aaaaaaaTGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCAc                                                                                           <  7:280738/149‑1 (MQ=255)
                                                   aaaaaaTGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATAAGTAGCCGGTGAGAATGATCAc                                                                                           >  6:241142/1‑148 (MQ=255)
                                                           gCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACAcc                                                                                  >  8:91956/1‑149 (MQ=255)
                                                              aaGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTa                                                                                                                                 >  8:150926/1‑99 (MQ=255)
                                                              aaGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTa                                                                                                                                 <  7:150926/99‑1 (MQ=255)
                                                               aGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTg                                                                              <  6:321386/149‑1 (MQ=255)
                                                                gATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTgg                                                                             <  6:340733/149‑1 (MQ=255)
                                                                   tCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTa                                                                                                                                 <  6:296716/94‑1 (MQ=255)
                                                                     cAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGTCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGTCAATg                                                                        <  5:241142/149‑1 (MQ=255)
                                                                     cAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAAt                                                                         <  2:133426/148‑1 (MQ=255)
                                                                     cAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATg                                                                        >  1:69884/1‑149 (MQ=255)
                                                                       aTTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACTA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGc                                                                      <  5:340842/149‑1 (MQ=255)
                                                                       aTTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGc                                                                      <  2:142941/149‑1 (MQ=255)
                                                                           cACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGccac                                                                   <  2:239673/148‑1 (MQ=255)
                                                                               gggTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGaa                                                              <  6:220285/149‑1 (MQ=255)
                                                                                     ggCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGccac                                                                   <  2:331915/138‑1 (MQ=255)
                                                                                     ggCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGccac                                                                   >  1:331915/1‑138 (MQ=255)
                                                                                         tATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGc                                                    >  2:345786/1‑148 (MQ=255)
                                                                                            cACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTa                                                                                                             >  8:233981/1‑89 (MQ=255)
                                                                                            cACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTa                                                                                                             <  7:233981/89‑1 (MQ=255)
                                                                                             aCCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTc                                               <  6:333877/149‑1 (MQ=255)
                                                                                               cTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATAAGTAGCCGTTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTca                                              >  6:116620/1‑148 (MQ=255)
                                                                                                      gtgAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAg                                       <  4:76043/148‑1 (MQ=255)
                                                                                                        gAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGccac                                                                   <  6:239502/119‑1 (MQ=255)
                                                                                                        gAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGccac                                                                   >  5:239502/1‑119 (MQ=255)
                                                                                                         aaCCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTga                                                                                                     <  8:71881/84‑1 (MQ=255)
                                                                                                         aaCCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTga                                                                                                     >  7:71881/1‑84 (MQ=255)
                                                                                                           ccTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAAc                                  >  1:307123/1‑148 (MQ=255)
                                                                                                           ccTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGAGCCTGCAAACAATCAGTAGCCGGTGACAATGTTCACGATCACACACCTGGCAATGGCCACCAGAACAAAGTCGCCATTCAAATACAGCCAAc                                  >  5:144234/1‑148 (MQ=255)
                                                                                                                    ggTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCaa                                   <  7:199578/140‑1 (MQ=255)
                                                                                                                    ggTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCaa                                   >  8:199578/1‑140 (MQ=255)
                                                                                                                      tgtAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCAGCCATTCACATACAGCCAACGATGATCa                          <  1:21317/148‑1 (MQ=255)
                                                                                                                      tgtAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGAATTTTAGTGCCTCCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCCCTCTCAATGGCCACCAGAACATACTCGCCATTCACATACAGCCAACGATGATAACg                        >  8:42429/1‑149 (MQ=255)
                                                                                                                        tAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGa                       <  7:91956/148‑1 (MQ=255)
                                                                                                                         aaTTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAgg                     >  1:119891/1‑149 (MQ=255)
                                                                                                                              gCCCGACCCTAACACCCCCGCGGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGgc                 <  2:280068/148‑1 (MQ=255)
                                                                                                                              gCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGgcg                <  2:69884/149‑1 (MQ=255)
                                                                                                                                cAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCa              >  6:53450/1‑149 (MQ=255)
                                                                                                                                   gcacTAAAAAGCCCGCAGCAGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGaa                                                              >  6:67669/2‑99 (MQ=255)
                                                                                                                                   gcacTAAAAAGCCCGCAGCAGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGaa                                                              <  5:67669/98‑1 (MQ=255)
                                                                                                                                     acTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAAGGTTCACGATCACACCGATGGCAATGGCCACCAGAACATAGTCGGCATTCACATACAGCCAACGATGATCACGAGGTGGAGCATACAg         >  3:18251/1‑149 (MQ=255)
                                                                                                                                      cTAAAAAGCCCGCAGCAGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAgg        >  6:323705/1‑149 (MQ=255)
                                                                                                                                       tAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAGGt       >  4:160262/1‑149 (MQ=255)
                                                                                                                                       tAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAGGt       >  3:334995/1‑149 (MQ=255)
                                                                                                                                       tAAAAAGCCCGCAGCAGCGGGCTTTTTAGTGTCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAAAATAGTCGCCAGTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAGGt       >  1:247227/1‑149 (MQ=39)
                                                                                                                                        aaaaaGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGGATACAGGTg      >  7:277938/1‑149 (MQ=255)
                                                                                                                                         aaaaGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAg         <  8:197712/145‑1 (MQ=255)
                                                                                                                                         aaaaGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAg         >  7:197712/1‑145 (MQ=255)
                                                                                                                                         aaaaGCCCGCAGCAGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAGGTgg     >  3:20443/1‑149 (MQ=255)
                                                                                                                                         aaaaGCCCGCAGCAGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAGGTgg     >  6:49115/1‑149 (MQ=255)
                                                                                                                                            aGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAGGTGgcgc  <  4:334995/149‑1 (MQ=255)
                                                                                                                                              cccGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTca                                              >  6:70191/1‑103 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
TGAGTAGGGATTTCTCGGCTTTTCTGTGAAAAGCACCCCTTCATTTCACAAAAAAAATGGCGAAGATTCCAATTTCACGGGGTCTGGCCTATCACCTATTTTGTGAACCTACAA‑‑GGTGTAATTCGCCAGACACTAAAAAGCCCGCTGCTGCGGGCTTTTTAGTGCCTGCAAACAATCAGTAGCCGGTGAGAATGTTCACGATCACACCGCTGGCAATGGCCACCAGAACATAGTCGCCATTCACATACAGCCAACGATGATCACGAGGTGGCGCATACAGGTGGCGC  >  NC_002947/58602‑58888

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