Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 58,439 G→C 8.4% intergenic (‑350/+338) 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,4390GC 8.4% 175.1 / 4.3 59intergenic (‑350/+338)PP_5424/PP_0050hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (27/27);  new base C (2/3);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

CTCCGTGAGGCCGGCCATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCTTTCCAGGCCAA  >  NC_002947/58291‑58584
                                                                                                                                                      |                                                                                                                                                 
cTCCGTGAGGCCGGCCATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGGATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTgg                                                                                                                                                   >  4:134554/1‑149 (MQ=255)
cTCCGTGAGGCCGGCCATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTgg                                                                                                                                                   <  2:43676/149‑1 (MQ=255)
cTCCGTGAGGCCGGCCATCGGCATGCATCAGATCTTGATTGACAGGCGTGCCGCTTTGGCTAGCCGTGCAATTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTgg                                                                                                                                                   <  4:120040/149‑1 (MQ=255)
 tCCGTGAGGCCGGCCATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCGCCGCAAGTTTGGCAGCAGAGGCAGCTGGt                                                                                                                                                  >  8:170174/1‑149 (MQ=255)
        ggccggccATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCctct                                                                                                                                            >  5:161601/1‑148 (MQ=255)
         gcaggccATAGTCATGCATCAGATTTTGATTGACACGCGTGCCTCTTTGGATAGCCTTGCAATTGCGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGa                                                                                                                                          <  7:170174/146‑1 (MQ=255)
          ccggccATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGaa                                                                                                                                         <  4:221355/149‑1 (MQ=255)
           cggccATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAACCAGCTTCCTCTgctg                                                                                                                                        <  5:330921/149‑4 (MQ=255)
             gccATCGTCATTCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGCTGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAag                                                                                                                                        <  1:111681/149‑1 (MQ=255)
             gccATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAagta                                                                                                                                      <  6:221731/149‑1 (MQ=255)
                   gTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGAcg                                                                                                                                >  3:49975/1‑149 (MQ=255)
                          aTCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCa                                                                                                                         >  3:139809/1‑149 (MQ=255)
                          aTCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCa                                                                                                                         >  5:373817/1‑149 (MQ=255)
                           tCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCAt                                                                                                                        >  2:60702/1‑149 (MQ=255)
                            cAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAACCAGCTTCCTCTgctg                                                                                                                                        <  5:104155/132‑4 (MQ=255)
                            cAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAACCAGCTTCCTCTgctg                                                                                                                                        >  6:104155/1‑129 (MQ=255)
                             aGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTg                                                                                                                      >  3:361197/1‑149 (MQ=255)
                                tCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGc                                                                                                                     >  6:208144/1‑147 (MQ=255)
                                    gattgaCACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTc                                                                                                                <  5:352055/148‑1 (MQ=255)
                                    gattgaCACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAACCAGCTTCCTCTgct                                                                                                                                         <  7:98218/123‑3 (MQ=255)
                                    gattgaCACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAACCAGCTTCCTCTgct                                                                                                                                         >  8:98218/1‑121 (MQ=255)
                                     attgaCACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCg                                                                                                               <  3:58997/148‑1 (MQ=255)
                                        gaCACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGAt                                                                                                            <  6:373817/148‑1 (MQ=255)
                                           acGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGa                                                                                                                                  >  2:360847/1‑123 (MQ=255)
                                             gcgTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATTAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGaa                                                                                                      <  7:321982/149‑1 (MQ=255)
                                             gcgTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGt                                                                                                                                                  >  6:285964/1‑105 (MQ=255)
                                             gcgTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGt                                                                                                                                                  <  5:285964/105‑1 (MQ=255)
                                                 gCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAAccc                                                                                                   >  6:367780/1‑148 (MQ=255)
                                                   cGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAg                                                                                                <  7:292542/149‑1 (MQ=255)
                                                   cGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGGTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAg                                                                                                <  6:164346/149‑1 (MQ=255)
                                                     cTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGt                                                                                              >  1:105226/1‑149 (MQ=255)
                                                        tGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTc                                                                                             <  4:49975/147‑1 (MQ=255)
                                                        tGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTc                                                                                             <  1:183331/147‑1 (MQ=255)
                                                          gATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCa                                                                                         <  3:58013/149‑1 (MQ=255)
                                                              gCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGAc                                                                                      <  8:126453/148‑1 (MQ=255)
                                                                 ttGCAATTGATGCCTTCAACTGCGCATTTGTCG‑‑GGCCTGGGGACCAGCATGAGCCCCCGCAATTATGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGTTCGCCAGACTAgg                                                                                  <  2:156318/149‑1 (MQ=255)
                                                                     aaTTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGc                                                                               >  5:113260/1‑148 (MQ=255)
                                                                     aaTTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTAAAGTAGTGACGGCGGGCGTTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGc                                                                               >  3:264341/1‑148 (MQ=255)
                                                                      aTTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCt                                                                              >  7:105800/1‑148 (MQ=255)
                                                                         gAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGTGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTaaa                                                                           <  6:115118/148‑1 (MQ=255)
                                                                         gAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCTGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTaaa                                                                           <  1:290733/148‑1 (MQ=255)
                                                                         gAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTaaa                                                                           <  6:54388/148‑1 (MQ=255)
                                                                           gggCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGa                                                                        <  4:361197/149‑1 (MQ=255)
                                                                               tGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGc                                                                    >  3:299299/1‑149 (MQ=255)
                                                                                  aaCTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGACCCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGc                                                                 <  6:266959/149‑1 (MQ=255)
                                                                                     tGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGac                                                              >  1:338588/1‑149 (MQ=255)
                                                                                       cgcATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGact                                                             >  4:83837/1‑148 (MQ=255)
                                                                                          aTTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTActc                                                         <  5:367780/149‑1 (MQ=255)
                                                                                           tttGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTActcc                                                        <  4:139809/149‑1 (MQ=255)
                                                                                           tttGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTActcc                                                        <  2:105226/149‑1 (MQ=255)
                                                                                              gTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAag                                                                                                                                        >  8:294127/1‑66 (MQ=255)
                                                                                                 g‑‑ggCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATcc                                                  <  3:83837/149‑1 (MQ=255)
                                                                                                    ggCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATccc                                                 >  3:242043/1‑149 (MQ=255)
                                                                                                     gCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGAAGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCt                                                >  6:85274/1‑149 (MQ=255)
                                                                                                            tACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGa                                          >  5:262612/1‑148 (MQ=255)
                                                                                                               cAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAg                                      >  5:62865/1‑149 (MQ=255)
                                                                                                                    tGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCAcg                                 <  6:161601/149‑1 (MQ=255)
                                                                                                                          cccGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGCCTACTCCTCATCCCTTACCGAAAAGTCACGCgtcgg                           >  6:32123/1‑149 (MQ=255)
                                                                                                                          ccagCAATTTTGGCAGAAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCACTTACCGAAAAGTCACGCttctg                           >  2:197582/4‑147 (MQ=255)
                                                                                                                                      gcagcaGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGcc               >  2:163947/1‑149 (MQ=255)
                                                                                                                                              ggAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAATCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCTTTCCAg       >  6:62656/1‑149 (MQ=255)
                                                                                                                                                   cTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCTTTCCAGGCCaa  <  4:299299/149‑1 (MQ=255)
                                                                                                                                                     ggTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCTTTCCAGGc     <  6:189725/144‑1 (MQ=255)
                                                                                                                                                     ggTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCTTTCCAGGc     >  5:189725/1‑144 (MQ=255)
                                                                                                                                                      gTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGt                       <  6:358588/125‑1 (MQ=255)
                                                                                                                                                      gTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGt                       >  5:358588/1‑125 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
CTCCGTGAGGCCGGCCATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCTTTCCAGGCCAA  >  NC_002947/58291‑58584

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