Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 344,032 A→G 12.3% intergenic (+121/+68) PP_5440 → / ← PP_t72 hypothetical protein/tRNA‑Phe

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947344,0320AG12.3% 157.2 / 12.8 65intergenic (+121/+68)PP_5440/PP_t72hypothetical protein/tRNA‑Phe
Reads supporting (aligned to +/- strand):  ref base A (30/27);  new base G (6/2);  total (36/29)
Fisher's exact test for biased strand distribution p-value = 2.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.05e-01

GACGAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGAC  >  NC_002947/343903‑344178
                                                                                                                                   |                                                                                                                                                   
gACGAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGt                                                                                                                                 >  3:53366/1‑149 (MQ=255)
     gTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCtt                                                                                                                            >  2:55964/1‑149 (MQ=255)
     gTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCtt                                                                                                                            <  8:219998/149‑1 (MQ=255)
       aaGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCTATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCttg                                                                                                                           <  3:244152/148‑1 (MQ=255)
        aGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAGAAACCCGCCTAGTGCGGGTTTTttgtt                                                                                                                         <  7:121354/149‑1 (MQ=39)
         ggTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCttgttg                                                                                                                        <  5:268351/149‑1 (MQ=255)
         ggTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCtt                                                                                                                            >  1:317980/1‑145 (MQ=255)
         ggTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCtt                                                                                                                            <  2:317980/145‑1 (MQ=255)
          gTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTtgttgt                                                                                                                       >  8:224781/1‑149 (MQ=255)
           tttGGCCCGCAGTCTGGTCCTGGCTGCTTGTTTCTCGATTGCATGGCTAGCCTTCAGGTGGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGCCGGGTTTCTTGTTGTc                                                                                                                      <  5:254007/149‑1 (MQ=255)
                 ccGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCg                                                                                                                >  2:277386/1‑149 (MQ=255)
                           cTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAGAAACCCGCCTAGTGCGGGTTTTTtgttgt                                                                                                                       >  2:39523/1‑132 (MQ=38)
                           cTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAGAAACCCGCCTAGTGCGGGTTTTTtgttgt                                                                                                                       <  1:39523/132‑1 (MQ=38)
                           cTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCttgt                                                                                                                          >  6:154351/1‑129 (MQ=255)
                           cTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCttgt                                                                                                                          <  5:154351/129‑1 (MQ=255)
                            tCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCttgt                                                                                                                          <  3:318293/128‑1 (MQ=255)
                            tCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCttgt                                                                                                                          >  4:318293/1‑128 (MQ=255)
                            tCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTg                                                                                                      >  5:189876/1‑148 (MQ=255)
                              aTGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCttgttg                                                                                                                        >  1:169280/1‑128 (MQ=255)
                              aTGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCttgttg                                                                                                                        <  2:169280/128‑1 (MQ=255)
                                 gCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAacca                                                                                                <  6:189876/149‑1 (MQ=255)
                                  cTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAaccac                                                                                               >  1:35085/1‑149 (MQ=255)
                                     cGGGTTTCTCGATTGCATGGCTTGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAacccc                                                                                               >  6:319652/1‑144 (MQ=255)
                                         tttCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGc                                                                                        <  1:277386/149‑1 (MQ=255)
                                             tcGATTGCATGGCTAGCCTTCAGGTTGCAGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCtt                                                                                      >  1:193811/1‑148 (MQ=255)
                                             tcGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTg                                                                                     >  7:296200/1‑148 (MQ=255)
                                               gATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCt                                                                                       <  1:188753/144‑1 (MQ=255)
                                               gATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCt                                                                                       >  2:188753/1‑144 (MQ=255)
                                               gATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGaa                                                                                   <  8:296200/148‑1 (MQ=255)
                                                   gCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAggtg                                                                               <  2:35085/148‑1 (MQ=255)
                                                    cATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAggtgg                                                                              <  4:53366/148‑1 (MQ=255)
                                                     aTGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGgtggtg                                                                            >  2:115548/1‑149 (MQ=255)
                                                     aTGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGgtggtg                                                                            >  3:113440/1‑149 (MQ=255)
                                                       ggCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGc                                                                           >  5:197759/1‑148 (MQ=255)
                                                       cgatttccTTCAGGTTGCGTTGTGATGTCTCGTGTTATTCGACGTCGCCTGAAGTTAAAAATCACAGAAAACAAAAATCGCGCAATTGGAGGGTTTCTTGTTGTCGCTTCGGGTATCTTTGCAACGACCAGAATGTTGTAGGTGGTGc                                                                           <  8:296132/142‑1 (MQ=255)
                                                           aGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCa                                                                        >  7:7975/1‑147 (MQ=255)
                                                              cTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTTTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGc                                                                           >  7:296132/1‑141 (MQ=255)
                                                                  aGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCTGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCTACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGa                                                                <  4:313597/148‑1 (MQ=255)
                                                                  aGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGa                                                                <  2:95788/148‑1 (MQ=255)
                                                                   ggTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGt                                                              <  6:197759/149‑1 (MQ=255)
                                                                    gTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTc                                                             >  4:320461/1‑149 (MQ=255)
                                                                              tgAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAGAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCgg                                                                                                               >  1:55696/1‑91 (MQ=255)
                                                                              tgAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCgg                                                                                                               <  2:55696/90‑1 (MQ=255)
                                                                                   tctcGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGacac                                                >  8:277107/1‑148 (MQ=255)
                                                                                       gtgtTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCg                                                            >  5:197290/1‑132 (MQ=255)
                                                                                       gtgtTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCg                                                            <  6:197290/132‑1 (MQ=255)
                                                                                            cttcGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGAtttt                                       <  3:320461/148‑1 (MQ=255)
                                                                                             ttcGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGa                                                   <  4:29777/135‑1 (MQ=255)
                                                                                             ttcGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTc                                      >  7:79213/1‑148 (MQ=255)
                                                                                               cGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAAt                                   >  7:17178/1‑149 (MQ=255)
                                                                                                gACGCCGCCTTAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATc                                  <  8:60485/149‑1 (MQ=255)
                                                                                                    ccgccTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCctct                              <  7:224781/149‑1 (MQ=255)
                                                                                                           aaGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTAc                       <  1:115548/149‑1 (MQ=255)
                                                                                                              ttAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGa                    >  2:32522/1‑149 (MQ=255)
                                                                                                                aaaaaTCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTAc                       >  6:243535/1‑144 (MQ=255)
                                                                                                                aaaaaTCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTAc                       <  5:243535/144‑1 (MQ=255)
                                                                                                                      tCACAGACAACAAAAAACCCGCACTAGGCGGTTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCATAGACGGAGTCGAACCGCCGACACGAGTATTTTCAATCCTCTGCTCTACCGACTGAGCt             <  1:90508/149‑1 (MQ=255)
                                                                                                                       cacaGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAGCCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTa            >  7:219777/1‑149 (MQ=255)
                                                                                                                           gACAACAAGAAACCCGCCTAGTGCGGGTTTTTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTg        >  5:152300/1‑149 (MQ=39)
                                                                                                                            acaacaAGAAACCCGCCTAGTGCGGGTTTTTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTgg       >  6:107856/1‑149 (MQ=39)
                                                                                                                            acaacaAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTgg       >  5:141304/1‑149 (MQ=255)
                                                                                                                             caacaaGAAACCCGCCTAGTGCGGGTTTTTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTggg      >  5:290421/1‑149 (MQ=39)
                                                                                                                              aacaaAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCAGCTCTACCGACTGAGCTAACTGGGc     >  7:51434/1‑149 (MQ=255)
                                                                                                                                caaGAAACCCGCCTAGTGCGGGTTTTTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGcga   >  7:14483/1‑149 (MQ=255)
                                                                                                                                caaGAAACCCGCCTAGTGCGGGTTTTTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGcga   >  7:242953/1‑149 (MQ=255)
                                                                                                                                  aaaaaCCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGcgac  >  3:287441/1‑148 (MQ=255)
                                                                                                                                   |                                                                                                                                                   
GACGAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGAC  >  NC_002947/343903‑344178

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