Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 719,424 C→T 6.6% A165T (GCA→ACA)  PP_0613 ← amidase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947719,4240CT 6.6% 180.4 / 2.5 60A165T (GCA→ACA) PP_0613amidase family protein
Reads supporting (aligned to +/- strand):  ref base C (24/32);  new base T (2/2);  total (26/34)
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.25e-01

CGGTCATGGCGCTGGCCGGCCATTGCGGCACATAGCCGAAGGTCGGTTTGATGCCGAAGGTGCCGGTGAAGGCGCAAGGGATGCGAATCGAACCACCGGCATCG‑‑CTGCCCTGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGAACCACCGGCAGTCAGCCTGGTGTCCCACGGGTTGCGAGTGATGCCGTACAGCGGGTTGTCGGTCACCCCTTTCCAGCCGAACTCTGGCGTGGTGGTCTTGCCCACCAGCACTGCACCGGCTTCACGCAT  >  NC_002947/719280‑719571
                                                                                                                                                  |                                                                                                                                                   
cGGTCATGGCGCTGGCCGGCCATTGCGGCACATAGCCGAAGGTCGGTTTGATGCCGAAGGTGCCGGTGAAGGCGCAAGGGATGCGAATCGAACCACCGGCATCG‑‑CTGCCCTGATGCAGCACCCCCAGGTCCAGCGCGGCGGCTGCCgca                                                                                                                                                 <  2:145879/149‑1 (MQ=255)
 ggTCATGGCGCTGGCCGGCCATTGCGGCACATAGCCGAAGGTCGGTTTGATGCCGAAGGTGCCGGTGAAGGCGCAAGGGATGCGAATCGCACCACCGGCATCG‑‑CTGCCCTGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGCCgcag                                                                                                                                                >  6:318166/1‑149 (MQ=255)
       ggCGCTGGCCGGCCATTGCGGCACATAGCCGAAGGTCGGTTTGATGCCGAAGGTGCCGGTGAAGGCGCAAGGGATGCGAATCGAACCACCGGCATCG‑‑CTGCCCTGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGCCGCAGcaccac                                                                                                                                          >  4:233465/1‑149 (MQ=255)
        gcgcTGGCCGGCCATTGCGGCACATAGCCGAAGGTCGGTTTGATGCCGAAGGTGCCGGTGAAGGCGCAAGGGATGCGAATCGAACCACCGGCATCG‑‑CTGCCCTGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGCCGCAGcaccac                                                                                                                                          <  5:232627/148‑1 (MQ=255)
              gccggccATTGCGGCACATAGCCGAAGGTCGGTTTGATGCCGAAGGTGCCGGTGAAGGCGCAAGGGATGCGAATCGAACCACCGGCATCG‑‑CTGCCCTGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACg                                                                                                                                    <  1:301075/148‑1 (MQ=255)
               ccggccATTGCGGCACATAGCCGAAGGTCGGTTTGATGCCGAAGGTGCCGGTGAAGGCGCAAGGGATGCGAATCGAACCACCGGCATCG‑‑CTGCCCTGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGaa                                                                                                                                  >  8:230379/1‑149 (MQ=255)
                cggccATTGCGGCACATAGCCGAAGGTCGGTTTGATGCCGAAGGTGCCGGTGAAGGCGCAAGGGATGCGAATCGAACCACCGGCATCG‑‑CTGCCCTGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGAac                                                                                                                                 <  8:149032/149‑1 (MQ=255)
                cggccATTGCGGCACATAGCCGAAGGTCGGTTTGATGCCGAAGGTGCCGGTGAAGGCGCAAGGGATGCGAATCGAACCACCGGCATCG‑‑CTGCCCTGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGAac                                                                                                                                 <  8:144938/149‑1 (MQ=255)
                 ggccATTGCGGCACATAGCCGAAGGTCGGTTTGATGCCGAAGGTGCCGGTGAAGGCGCAAGGGATGCGAATCGAACCACCGGCATCG‑‑CTGCCCTGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGAac                                                                                                                                 >  7:174515/1‑148 (MQ=255)
                   ccATTGCGGCACATAGCCGAAGGTCGGTTTGATGCCGAAGGTGCCGGTGAAGGCGCAAGGGATGCGAATCGAACCACCGGCATCG‑‑CTGCCCTGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGAaccac                                                                                                                              <  3:266401/149‑1 (MQ=255)
                           gCACATAGCAGAAGGTCGGTTTGATGCCGAAGGTGCCGGTGAAGGCGCAAGGGATGCGAATCGAACCACCGGCATCG‑‑CTGCCCTGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGAACCACCGGCAGTc                                                                                                                      <  2:237670/149‑1 (MQ=255)
                               aTAGCCGAAGGTCGGTTTGATGCCGAAGGTGCCGGTGAAGGCGCAAGGGATGCGAATCGAACCACCGGCATCG‑‑CTGCCCTGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGCACCACCGGCAGTCAGcc                                                                                                                  >  8:96935/1‑149 (MQ=255)
                               aTAGCCGAAGGTCGGTTTGATGCCGAAGGTGCCGGTGAAGGCGCAAGGGATGCGAATCGAACCACCGGCATCG‑‑CTGCCCTGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGAACCACCGGCAGTCAGcc                                                                                                                  <  1:40435/149‑1 (MQ=255)
                                tAGCCGAAGGTCGGTTTGATGCCGAAGGTGCCGGTGAAGGCGCAAGGGATGCGAATCGAACCACCGGCATCG‑‑CTGCCCTGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGAACCACCGGCAGTCAGCCt                                                                                                                 <  7:230379/149‑1 (MQ=255)
                                  tccGAAGGTCGGTTTGATTCCGAAGGTGCCGGTGAATGCGCAAGGGATGCGAATCGAATCAACGGCATCGGTCTTTTCTGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGAACCACCGGCAGTCAGcc                                                                                                                  <  1:256344/147‑1 (MQ=255)
                                  gCCGAAGGTCGGTTTGATGCCGAAGGTGCCGGTGAAGGCGCAAGGGATGCGAATCGAACCACCGGCATCG‑‑CTGCCCTGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGAACCACCGGCAGTCAGCCTgg                                                                                                               <  3:233465/149‑1 (MQ=255)
                                  gCCGAAGGTCGGTTTGATGCCGAAGGTGCCGGTGAAGGCGAAAGGGATGCGAATCGAACCACCGGCATCG‑‑CTGCCCTGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGAACCACCGGCAGTCAGCCTg                                                                                                                <  7:296207/148‑1 (MQ=255)
                                    cGAAGGTCGGTTTGATGCCGAAGGTGCCGGTGAAGGCGCAAGGGATGCGAATCGAACCACCGGCATCG‑‑CTGCCCTGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGAACCACCGGCAGTCAGCCTGgtg                                                                                                             >  6:214112/1‑149 (MQ=255)
                                      aaGGTCGGTTTGATGCCGAAGGTGCCGGTGAAGGCGCAAGGGATGCGAATCGAACCACCGGCATCG‑‑CTGCCCTGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGAACCACCGGCAGTCAGCCTGGTGTc                                                                                                           <  8:174515/149‑1 (MQ=255)
                                              tttGATGCCGAAGGTGCCGGTGAAGGCGCAAGGGAAGCGAATCGAACCACCGGCATCG‑‑CTGCCCTGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGAACCACCGGCAGTCAGCCTGGTGTCCCACGGGt                                                                                                   >  1:207098/1‑149 (MQ=255)
                                                 gATGCCGAAGGTGCCGGTGAAGGCGCAAGGGATGCGAATCGAACCACCGGCATCG‑‑CTGCCCTGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGAACCACCGGCAGTCAGCCTGGTGTCCCACGGGTTGc                                                                                                >  5:323812/1‑149 (MQ=255)
                                                     ccGAAGGTGCCGGTGAAGGCGCAAGGGATGCGAATCGAACCACCGGCATCG‑‑CTGCCCTGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGAACCACCGGCAGTCAGcc                                                                                                                  <  2:31414/127‑1 (MQ=255)
                                                     ccGAAGGTGCCGGTGAAGGCGCAAGGGATGCGAATCGAACCACCGGCATCG‑‑CTGCCCTGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGAACCACCGGCAGTCAGcc                                                                                                                  <  4:249359/127‑1 (MQ=255)
                                                     ccGAAGGTGCCGGTGAAGGCGCAAGGGATGCGAATCGAACCACCGGCATCG‑‑CTGCCCTGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGAACCACCGGCAGTCAGcc                                                                                                                  >  3:249359/1‑127 (MQ=255)
                                                     ccGAAGGTGCCGGTGAAGGCGCAAGGGATGCGAATCGAACCACCGGCATCG‑‑CTGCCCTGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGAACCACCGGCAGTCAGcc                                                                                                                  >  1:31414/1‑127 (MQ=255)
                                                       gAAGGTGCCGGTGAAGGCGCAAGGGATGCGAATCGAACCACCGGCATCG‑‑CTGCCCTGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGTCGCAGCACCACCAGACGAACCACCGGCAGTCAGCCTGGTGTCCCACGGGTTGCGAGTGa                                                                                          >  6:76134/1‑149 (MQ=255)
                                                       gAAGGTGCCGGTGAAGGCGCAAGGGATGCGAATCGAACCACCGGCATCG‑‑CTGCCCTGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGTCGCAGCACCACCAGACGAACCACCGGCAGTCAGCCTGGTGTCCCACGGGTTGCGAGTGa                                                                                          >  6:80642/1‑149 (MQ=255)
                                                           gTGCCGGTGAAGGCGCAAGGGATGCGAATCGAACCACCGGCATCG‑‑CTGCCCTGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGAACCACCGGCAGTCAGCCTGGTGTCCCACGGGTTGCGAGTGATGcc                                                                                      <  7:191684/149‑1 (MQ=255)
                                                                ggTGAAGGCGCAAGGGATGCGAATCGAACCACCGGCATCG‑‑CTGCCCTGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGAACCACCGGCAGTCAGCCTGGTGTCCCACGGGTTGCGAGTGATGCCGTa                                                                                   >  4:181495/1‑147 (MQ=255)
                                                                    aaGGCGCAAGGGATGCGAATCGAACCACCGGCATCG‑‑CTGCCCTGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGCCGCCGCACCACCAGACGAACCACCGGCAGTCAGCCTGGTGTCCCACGGGTTGCGAGTGATGCCGTACAGCgg                                                                             >  2:264692/1‑149 (MQ=255)
                                                                     aGGCGCAAGGGATGCGAATCGAACCACCGGCATCG‑‑CTGCCCTGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGAACCACCGGCAGTCAGCCTGGTGTCCCACGGGTTGCGAGTGATGCCGTACAGCggg                                                                            <  7:96935/149‑1 (MQ=255)
                                                                       gcgcAAGGGATGCGAATCGAACCACCGGCATCG‑‑CTGCCCTGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCATACGAACCACCGGCAGTCACCCACGAGTCCCACACGTTGCGAGTGATCCCGTACAGCGGGtt                                                                          >  8:283846/1‑149 (MQ=255)
                                                                       gcgcAAGGGATGCGAATCGAACCACCGGCATCG‑‑CTGCCCTGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGAACCACCGGCAGTCAGCCTGGTGTCCCACGGGTTGCGAGTGATGCCGTACAGCggg                                                                            >  1:13001/1‑147 (MQ=255)
                                                                       gcgcAAGGGATGCGAATCGAACCACCGGCATCG‑‑CTGCCCTGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGAACCACCGGCAGTCAGCCTGGTGTCCCACGGGTTGCGAGTGATGCCGTACAGCggg                                                                            <  2:13001/147‑1 (MQ=255)
                                                                                       tCGAACCACCGGCATCG‑‑CTGCCCTGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGAACCACCGGCAGTCAGCCTGGTGTCCCACGGGTTGCGAGTGATGCCGTACAGCGGGTTGTCGGTCACCCCTTTc                                                          <  6:323812/149‑1 (MQ=255)
                                                                                       tCGAACCACCGGCATCG‑‑CTGCCCTGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGAACCACCGGCAGTCAGCCTGGTGTCCCACGGGTTGCGAGTGATGCCGTACAGCGGGTTGTCGGTCACCCCTTTc                                                          <  4:108852/149‑1 (MQ=255)
                                                                                       tCGAACCACCGGCATCG‑‑CTGCCCTGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGAACCACCGGCAGTCAGCCTGGTGTCCCACGGGTTGCGAGTGATGCCGTACAGCGGGTTGTCGGTCACCCCTTTc                                                          >  1:248481/1‑149 (MQ=255)
                                                                                            ccaccGGCATCG‑‑CTGCCCTGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGAACCACCGGCAGTCAGCCTGGTGTCCCACGGGTTGCGAGTGATGCCGTACAGCGGGTTGTCGGTCACCCCTTTCCAGcc                                                     >  4:272350/1‑149 (MQ=255)
                                                                                            ccaccGGCATCG‑‑CTGCCCTGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGAACCACCGGCAGTCAGCCTGGTGTCCCACGGGTTGCGAGTGATGCCGTACAGCGGGTTGTCGGTCACCCCTTTCCAGcc                                                     >  2:146456/1‑149 (MQ=255)
                                                                                                   cATCG‑‑CTGCCCTGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGAACCACCGGCAGTCAGCCTGGTGTCCCACGGGTTGCGAGTGATGCCGTACAGCGGGTTGTCGGt                                                                    <  4:129142/127‑1 (MQ=255)
                                                                                                   cATCG‑‑CTGCCCAGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGAACCACCGGCAGTCAGCCTGGTGTCCCACGGGTTGCGAGTGATGCCGTACAGCGGGTTGTCGGTCACCCCTTTCCAGCCGAActct                                              >  5:4051/1‑149 (MQ=255)
                                                                                                             cccTGATGCAGCACCCCCAGGTTCCGCGCGGCGGCTGCCGCAGCACCACCAGACGAACCACCGGCAGTCAGCCTGGTGTCCCACGGGTTGCGAGTGATGCCGTACAGCGGGTTGTCGGTCACCGCTTTCCAGCCGAACTCTGGCgtggt                                      >  7:258965/1‑149 (MQ=255)
                                                                                                             cccTGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGAACCACCGGCAGTCAGCCTGGTGTCCCACGGGTTGCGAGTGATGCCGTACAGCGGGTTGTCGGTCACCCCTTTCCAGCCGAACTCTGGCgtg                                        <  5:318166/147‑1 (MQ=255)
                                                                                                                tGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGAACCACCGGCAGTCAGCCTGGTGTCCCACGGGTTGCGAGTGATGCCGTACAGCggg                                                                            >  1:282917/1‑108 (MQ=255)
                                                                                                                tGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGAACCACCGGCAGTCAGCCTGGTGTCCCACGGGTTGCGAGTGATGCCGTACAGCggg                                                                            <  2:282917/108‑1 (MQ=255)
                                                                                                                   tGCAGCACCCCCAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGAACCACCGGCAGTCAGCCTGGTGTCCCACGGGTTGCGAGTGATGCCGTACAGCGGGTTGTCGGTCACCCCTTTCCAGCCGAACTCTGGCGTGGTGGTCtt                                <  1:146456/149‑1 (MQ=255)
                                                                                                                    gcagcaCCCCCAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGAACCACCGGCAGTCAGCCTGGTGTCCCACGGGTTGCGAGTGATGCCGTACAGCGGGTTGTCGGTCACCCCTTTCCAGCCGAACTCTGGCGTGGTGGTCTTg                               >  2:93164/1‑149 (MQ=255)
                                                                                                                         aCCCCCAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGAACCACCGGCAGTCAGCCTGGTGTCCCACGGGTTGCGAGTGATGCCGTACAGCGGGTTGTCGGTCACCCCTTTCCAGCCGAACTCTGGCGTGGTGGTCTTGCccac                          <  6:5252/149‑1 (MQ=255)
                                                                                                                           ccccAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGAACCACCGGCAGTCAGCCTGGTGTCCCACGGGTTGCGAGTGATGCCGTACAGCGGGTTGTCGGTCACCCCTTTCCAGCCGAACTCTGGCGTGGTGGTCTTGCccacca                        <  3:181495/149‑1 (MQ=255)
                                                                                                                            cccAGGTTCAGCGCGGCGGCTGTCGCAGCACCACCAGACGAACCACCGGCAGTCAGCCTGGTGTCCCACGGGTTGCGAGTGATGCCGTACAGCGGGTTGTCGGTCACCCCTTTCCAGCCGAACTCTGGCGTGGTGGTCTTGCCCACCAg                       <  5:80642/149‑1 (MQ=255)
                                                                                                                            cccAGGTTCAGCGCGGCGGCTGTCGCAGCACCACCAGACGAACCACCGGCAGTCAGCCTGGTGTCCCACGGGTTGCGAGTGATGCCGTACAGCGGGTTGTCGGTCACCCCTTTCCAGCCGAACTCTGGCGTGGTGGTCTTGCCCACCAg                       <  5:76134/149‑1 (MQ=255)
                                                                                                                            cccAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGAACCACCGGCAGTCAGCGTGGTGTCCCACGGGTTGCGAGTGATGCCGTACAGCGGGTTGTCGGTCACCCCTTTCCAGCCGAACTCTGGCGTGGTGGTCTTGCccac                          >  1:170036/1‑146 (MQ=255)
                                                                                                                            cccAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGAACCACCGGCAGTCAGCCTGGTGTCCCACGGGTTGCGAGTGATGCCGTACAGCGGGTTGTCGGTCACCCCTTTCCAGCCGAACTCTGGCGTGGTGGTCTTGCccacca                        <  2:170036/148‑1 (MQ=255)
                                                                                                                             ccAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGAACCACCGGCAGTCAGCCTGGTGTCCCACGGGTTGCGAGTGATGCCGTACAGCGGGTTGTCGGTCACCCCTTTCCAGCCGAACTCTGGCGTggtggt                                   <  2:21119/136‑1 (MQ=255)
                                                                                                                             ccAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGAACCACCGGCAGTCAGCCTGGTGTCCCACGGGTTGCGAGTGATGCCGTACAGCGGGTTGTCGGTCACCCCTTTCCAGCCGAACTCTGGCGTggtggt                                   >  1:21119/1‑136 (MQ=255)
                                                                                                                                ggTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGAACCACCGGCAGTCAGCCTGGTGTCCCACGGGTTGCGAGTGATGCCGTACAGCGGGTTGTCGGTCACCCCTTTCCAGCCGAACTCTGGCGTGGTGGGCAAGAccaccacgcca                   <  3:272350/149‑6 (MQ=255)
                                                                                                                                        gcggcgGCTGCCGCAGCACCACCAGACGAACCACCGGCAGTCAGCCTGGTGTCCCACGGGTTGCGAGTGATGCCGTACAGCGGGTTGTCGGTCACCCCTTTCCAGCCGAACTCTGGCGTGGTGGTCTTGCCCACCAGCACTGCACCGGc           <  1:93164/149‑1 (MQ=255)
                                                                                                                                          ggcggcTGCCGCAGCACCACCAGACGAACCACCGGCAGTCAGCCTGGTGTCCCACGGGTTGCGAGTGATGCCGTACAGCGGGTTGTCGGTCACCCCTTTCCAGCCGAACTCTGGCGTGGTGGTCTTGCCCACCAGCACTGCACCGGCtt         <  2:207098/149‑1 (MQ=255)
                                                                                                                                               cTGCCGCAGCACCACCAGACGAACCACCGGCAGTCAGCCTGGTGTCCCACGGGTTGCGAGTGATGCCGTACAGCGGGTTGTCGGTCACCCCTTTCCAGCCGAACTCTGGCGTGGTGGTCTTGCCCACCAGCACTGCACCGGCTTCACGc    <  2:248481/149‑1 (MQ=255)
                                                                                                                                                tGCCGCAGCACCACCAGACGAACCACCGGCAGTCAGCCTGGTGTCCCACGGGTTGCGAGTGATGCCGTACAGCGGGTTGTCGGTCACCCCTTTCCAGCCGAACTCTGGCGTGGTGGTCTTGCCCACCAGCACTGCACCGGCTTCACGCa   >  2:278708/1‑149 (MQ=255)
                                                                                                                                                  ccGCAGCACCACCAGACGAACCACCGGCAGTCAGCCTGGTGTCCCACGGGTTGCGAGTGATGCCGTACAGCGGGTTGTCGGTCACCCCTTTCCAGCCGAACTCTGGCGTGGTGGTCTTGCCCACCAGCACTGCACCGGCTTCACGCAt  >  8:162585/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CGGTCATGGCGCTGGCCGGCCATTGCGGCACATAGCCGAAGGTCGGTTTGATGCCGAAGGTGCCGGTGAAGGCGCAAGGGATGCGAATCGAACCACCGGCATCG‑‑CTGCCCTGATGCAGCACCCCCAGGTTCAGCGCGGCGGCTGCCGCAGCACCACCAGACGAACCACCGGCAGTCAGCCTGGTGTCCCACGGGTTGCGAGTGATGCCGTACAGCGGGTTGTCGGTCACCCCTTTCCAGCCGAACTCTGGCGTGGTGGTCTTGCCCACCAGCACTGCACCGGCTTCACGCAT  >  NC_002947/719280‑719571

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: