Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,390,480 C→A 6.9% intergenic (+521/‑201) PP_4738 → / → PP_mr56 hypothetical protein/C4 ncRNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,390,4800CA 6.9% 229.1 / 4.6 72intergenic (+521/‑201)PP_4738/PP_mr56hypothetical protein/C4 ncRNA
Reads supporting (aligned to +/- strand):  ref base C (40/27);  new base A (2/3);  total (42/30)
Fisher's exact test for biased strand distribution p-value = 6.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

ACTCATCCAAAGCACCCTGCCCGTCGGGCGGCTGCGGATTAAATCAGTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTT  >  NC_002947/5390337‑5390629
                                                                                                                                               |                                                                                                                                                     
aCTCATCCAAAGCACCCTGCCCGTCGGGCGGCTGCGGGTTAAATCAGTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGACCAGCATATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGt                                                                                                                                                 >  1:979420/1‑150 (MQ=255)
   cATCCAAATCACCCTGCCCGTCGGGCGGCTGCGGATTAAATCAGTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGACGACGCCCGAGGGCGTCtt                                                                                                                                              <  1:346450/150‑1 (MQ=255)
   cATCCAAAGCACCCTGCCCGTCGGGCGGCTGCGGATTAAATTCGCGTGCCGGGGTTATCATTTAGTACCGACCGCCGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCtt                                                                                                                                              <  1:871533/150‑1 (MQ=255)
     tCCAAAGCACCCTGCCCGTCGGGCGGCTGCGGATTAAATCAGTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCtttt                                                                                                                                            >  2:505946/1‑150 (MQ=255)
     tCCAAAGCACCCTGCCCGTCGGGCGGCTGCGGATTAAATCAGTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCtttt                                                                                                                                            >  1:710807/1‑150 (MQ=255)
     tCCAAAGCACCCTGCCCGTCGGGCGGCTGCGGATTAAATCAGTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCtttt                                                                                                                                            >  2:1144485/1‑150 (MQ=255)
      ccAAAGCACCCTGCCCGTCGGGCGGCTGCGGATTAAATCAGTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCttttt                                                                                                                                           >  1:700955/1‑150 (MQ=255)
        aaaGCACCCTGCCCGTCGGGCGGCTGCGGATTAAATCAGTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCttttttt                                                                                                                                         <  1:10632/150‑1 (MQ=255)
             aCCCTGCCCGTCGGGCGGCTGCGGATTAAATCAGTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCGAGGGCGTCttc                                                                                                                                             <  2:753937/141‑2 (MQ=255)
             aCCCTGCCCGTCGGGCGGCTGCGGATTAAATCAGTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCGAGGGCGTCttc                                                                                                                                             >  1:753937/1‑140 (MQ=255)
               ccTGCCCGTCGGGCGGCTGCGGATTAAATCAGTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTg                                                                                                                                  >  2:676037/1‑150 (MQ=255)
               ccTGCCCGTCGGGCGGCTGCGGATTAAATCAGTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCGAGGGCGTCttc                                                                                                                                             >  2:440339/1‑138 (MQ=255)
               ccTGCCCGTCGGGCGGCTGCGGATTAAATCAGTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCGAGGGCGTCttc                                                                                                                                             <  1:440339/139‑2 (MQ=255)
                 tGCCCGTCGGGCGGCTGCGGATTAAATCAGTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGa                                                                                                                                >  1:189212/1‑150 (MQ=255)
                       tCGGGCGGCTGCGGATTAAATCAGTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTgg                                                                                                                          <  1:530870/150‑1 (MQ=255)
                           gcggcTGCGGATTAAATCAGTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGgcgc                                                                                                                      >  2:1063888/1‑150 (MQ=255)
                                tGCGGATTAAATCAGTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCgg                                                                                                                 >  2:295726/1‑150 (MQ=255)
                                    gATTAAATCAGTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCg                                                                                                             >  1:267670/1‑150 (MQ=255)
                                         aaTCAGTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGaaaaa                                                                                                        >  1:1049480/1‑150 (MQ=255)
                                          aTCAGTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAAt                                                                                                       >  2:263777/1‑150 (MQ=255)
                                           tCAGTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCgg                                                                                                                 <  2:798548/139‑1 (MQ=255)
                                           tCAGTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCgg                                                                                                                 >  1:798548/1‑139 (MQ=255)
                                              gTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGa                                                                                                   <  1:1144485/150‑1 (MQ=255)
                                                  acacGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTa                                                                                               >  2:527550/1‑150 (MQ=255)
                                                        cTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGa                                                                                                   >  1:615321/1‑140 (MQ=255)
                                                        cTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGa                                                                                                   <  2:615321/140‑1 (MQ=255)
                                                           aGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTc                                                                                      <  2:267670/150‑1 (MQ=255)
                                                                 tagtaCCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGaa                                                                                                                               >  2:1089460/1‑103 (MQ=255)
                                                                 tagtaCCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGaa                                                                                                                               <  1:1089460/103‑1 (MQ=255)
                                                                 tagtaCCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAg                                                                                <  1:297511/150‑1 (MQ=255)
                                                                 tagtaCCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAg                                                                                <  1:899090/150‑1 (MQ=255)
                                                                 tagtaCCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAg                                                                                <  2:710807/150‑1 (MQ=255)
                                                                  agtaCCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAgg                                                                               <  1:295726/150‑1 (MQ=255)
                                                                    taCCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGc                                                                             >  1:428279/1‑150 (MQ=255)
                                                                         aCAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAGATCGGATTTAACCCCTTTCGTGTAGGGCGTTTc                                                                        >  1:208721/1‑150 (MQ=255)
                                                                           aGCGGAGTACTGGCGGACGGGGGTTCAAATCCCGACGNGTCAGACGGNGATGGTGAACAATAGGCCATCGTTGGCTTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTTTAGGGCGttttgc                                                                      >  1:328604/1‑147 (MQ=255)
                                                                             cGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCaa                                                                    >  2:352995/1‑150 (MQ=255)
                                                                                 gTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATa                                                               <  1:505946/150‑1 (MQ=255)
                                                                                       gCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATActgcc                                                          >  2:344996/1‑150 (MQ=255)
                                                                                              ggggTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCt                                                   <  1:344996/150‑1 (MQ=255)
                                                                                               gggTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGaa                                                                                                           >  2:966809/1‑93 (MQ=255)
                                                                                               gggTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGaa                                                                                                           <  1:966809/93‑1 (MQ=255)
                                                                                                  ttCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGttt                                               >  1:953871/1‑150 (MQ=255)
                                                                                                      aaTCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCg                                           >  2:207429/1‑150 (MQ=255)
                                                                                                           ccccGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCg                                      >  2:751703/1‑150 (MQ=255)
                                                                                                                 cTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTgcgc                                                     <  1:431403/129‑1 (MQ=255)
                                                                                                                 cTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTgcgc                                                     >  2:431403/1‑129 (MQ=255)
                                                                                                                    caccaAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCAttttt                             >  1:1014425/1‑150 (MQ=255)
                                                                                                                     accaAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCAtttttt                            >  1:180682/1‑150 (MQ=255)
                                                                                                                      ccaAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCg                                                                                                                  >  1:947713/1‑63 (MQ=255)
                                                                                                                      ccaAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCg                                                                                                                  <  2:947713/63‑1 (MQ=255)
                                                                                                                      ccaAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTATTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTc                           >  1:997475/1‑150 (MQ=255)
                                                                                                                       caAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCa                          >  2:1074932/1‑150 (MQ=255)
                                                                                                                            ccAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGtt                     <  1:751703/150‑1 (MQ=255)
                                                                                                                            ccAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGtt                     <  1:263777/150‑1 (MQ=255)
                                                                                                                              aaaGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGAtt                                                                                                 <  2:1096958/72‑1 (MQ=255)
                                                                                                                              aaaGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGAtt                                                                                                 >  1:1096958/1‑72 (MQ=255)
                                                                                                                              aaaGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGtttt                   >  1:839417/1‑150 (MQ=255)
                                                                                                                              aaaGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGtttt                   >  1:202312/1‑150 (MQ=255)
                                                                                                                              aaaGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCACGGtttt                   >  2:434108/1‑150 (MQ=255)
                                                                                                                                    gaagaCGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGtctc             >  1:1100635/1‑150 (MQ=255)
                                                                                                                                       gaCGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAgg                                                                               <  2:421683/81‑1 (MQ=255)
                                                                                                                                       gaCGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAgg                                                                               >  1:421683/1‑81 (MQ=255)
                                                                                                                                       gaCGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGt          <  2:428279/150‑1 (MQ=255)
                                                                                                                                       gaCGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGt          <  1:1063888/150‑1 (MQ=255)
                                                                                                                                            ccTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTCACCCCTTTCTTGTAGGGCGTTTCCCAAAg                                                                  <  1:937830/89‑1 (MQ=255)
                                                                                                                                            ccTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAAc                                                                                             <  1:926164/62‑1 (MQ=255)
                                                                                                                                            ccTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAAc                                                                                             >  2:926164/1‑62 (MQ=255)
                                                                                                                                            ccTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAg                                                                  >  2:937830/1‑89 (MQ=255)
                                                                                                                                            ccTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTc     <  1:352995/150‑1 (MQ=255)
                                                                                                                                              tCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGt   >  1:786018/1‑150 (MQ=255)
                                                                                                                                              tCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGt   <  1:527550/150‑1 (MQ=255)
                                                                                                                                               cGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCCGAGTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGtt  <  1:676037/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
ACTCATCCAAAGCACCCTGCCCGTCGGGCGGCTGCGGATTAAATCAGTAGACACGGCTAAGCATGTAGTACCGACAGCGGAGTACTGGCGGACGGGGGTTCAAATCCCCCCGGCTCCACCAAATCCAAAGAAGAAGACGCCCTCGGGCGTCTTTTTTTGTGCCTGGAATTTGGGCGCTGCGGGTCGAAAAATCGGATTTAACCCCTTTCTTGTAGGGCGTTTCCCAAAGATACTGCCTGCGCCTGTTTTCCGTTGCGGTCATTTTTTCAGGGTTTTAGTCTCGGTGGGTCGTT  >  NC_002947/5390337‑5390629

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