Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 58,437 T→A 7.8% intergenic (‑348/+340) 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,4370TA 7.8% 180.6 / 5.1 64intergenic (‑348/+340)PP_5424/PP_0050hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (30/29);  new base A (2/3);  total (32/32)
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.03e-01

CCTCCGTGAGGCCGGCCATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCTTTCCAGGCCAA  >  NC_002947/58290‑58584
                                                                                                                                                     |                                                                                                                                                   
ccTCCGTGAGGCCGGCCATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTg                                                                                                                                                    >  5:54388/1‑149 (MQ=255)
actacGTGAGGCCGGCCATCGTCATGCATCAGATCTTGATTGACACGCGTGACGCTTTGGATAGCCTTGAAATTGAGGGCTGCAACTGCGCATTAGTCG‑‑GGCATGGGTACCAGAATGAGCCCCACCAATTTTGGCAGCAGAGGAAGATg                                                                                                                                                    >  8:321982/5‑149 (MQ=255)
 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                                                                                                                                        >  6:104155/1‑129 (MQ=255)
                             cAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAACCAGCTTCCTCTgctg                                                                                                                                        <  5:104155/132‑4 (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                                                                                                                                                  <  5:285964/105‑1 (MQ=255)
                                              gcgTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGt                                                                                                                                                  >  6:285964/1‑105 (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)
                                                                ccTTGCAATTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCt                                                                                                                                                     <  7:118853/84‑1 (MQ=255)
                                                                ccTTGCAATTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCt                                                                                                                                                     >  8:118853/1‑84 (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                                                        <  2:105226/149‑1 (MQ=255)
                                                                                            tttGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTActcc                                                        <  4:139809/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)
                                                                                                                                                     |                                                                                                                                                   
CCTCCGTGAGGCCGGCCATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCG‑‑GGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCTTTCCAGGCCAA  >  NC_002947/58290‑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: