Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,557,993 A→C 5.6% S55A (TCG→GCG)  PP_4045 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,557,9930AC 5.6% 202.9 / 2.0 67S55A (TCG→GCG) PP_4045hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (31/32);  new base C (2/2);  total (33/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 = 9.66e-01

GCCCATTCCCTGTGGGGGCTGAACCGCTCCCACTTCGCTTAATCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAGGGTGACCAAGCAGAACGATGGCAG  >  NC_002947/4557851‑4558142
                                                                                                                                              |                                                                                                                                                     
gCCCATTCCCTGTGGGGGCTGAACCGCTCCCACTTCGCTTAATCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCt                                                                                                                                                 >  2:171285/1‑149 (MQ=255)
 cccATTCCCTGTGGGGGCTGAACCGCTCCCACTTCGCTTAATCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTg                                                                                                                                                >  3:231434/1‑149 (MQ=255)
  ccATTCCCTGTGGGGGCTGAACCGCTCCCACTTCGCTTAATCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTgg                                                                                                                                               <  8:127058/149‑1 (MQ=255)
  ccATTCCCTGTGGGGGCTGAACCGCTCCCACTTCGCTTAATCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTgg                                                                                                                                               <  6:175398/149‑1 (MQ=255)
         cTGTGGGGGCTGAACCGCTCCCACTTCGCTTAATCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCACCAGCGGCTCGGCGCTgatgatgga                                                                                                                                        <  2:170759/149‑4 (MQ=255)
              ggggCTGAACCGCTCCCACTTCGCTTAATCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGccc                                                                                                                                   <  4:263634/149‑1 (MQ=255)
                          cTCCCACTTCGCTTAATCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGcc                                                                                                                       >  4:309310/1‑149 (MQ=255)
                           tCCCACTTCGCTTAATCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCt                                                                                                                      >  6:115098/1‑149 (MQ=255)
                               aCTTCGCTTAATCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCg                                                                                                                  >  3:140202/1‑149 (MQ=255)
                                    gCTTAATCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCaa                                                                                                              <  4:140202/148‑1 (MQ=255)
                                    gCTTAATCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAAc                                                                                                             >  3:77958/1‑149 (MQ=255)
                                      ttAATCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAAccc                                                                                                           >  8:265668/1‑149 (MQ=255)
                                        aaTCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATcgccgc                                                                                                                                     <  2:311212/121‑1 (MQ=255)
                                        aaTCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCACCAGCGGCTCGGCGCtgatga                                                                                                                                            <  7:108945/114‑1 (MQ=255)
                                        aaTCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCACCAGCGGCTCGGCGCtgatga                                                                                                                                            >  8:108945/1‑114 (MQ=255)
                                        aaTCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGAAGCATCGACTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATcgccgc                                                                                                                                     >  1:311212/1‑121 (MQ=255)
                                         aTCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTgg                                                                                                        <  1:270513/149‑1 (MQ=255)
                                                    aacTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTg                                                                                             >  6:54750/1‑149 (MQ=255)
                                                    aacTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTg                                                                                             >  6:49355/1‑149 (MQ=255)
                                                      cTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGcc                                                                                           <  5:115098/149‑1 (MQ=255)
                                                       tGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCTACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCg                                                                                          >  5:52856/1‑149 (MQ=255)
                                                         ttATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGtt                                                                                        <  2:238828/149‑1 (MQ=255)
                                                         atgtCGCGATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGtt                                                                                        <  7:293133/146‑1 (MQ=255)
                                                            tCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGgcc                                                                                                   <  4:338537/135‑1 (MQ=255)
                                                            tCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGgcc                                                                                                   >  3:338537/1‑135 (MQ=255)
                                                            tCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGa                                                                                      <  5:179148/148‑1 (MQ=255)
                                                            tCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGa                                                                                      <  1:171285/148‑1 (MQ=255)
                                                            tCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACACTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCACCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGAc                                                                                     >  2:238761/1‑149 (MQ=255)
                                                              gCCATTTTCCAGCACGCCGCCATGACTGGTAGCATCGCCAACCCTAGAAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCTCCGCCCCACAGCAGAGCCTGTCGCCAACCCTGGCGGCCGCCATGCCgat                                                                                        >  8:293133/1‑142 (MQ=255)
                                                               ccATTTTCCAGGACGCCGCCATGACTGGTCGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGCGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACgat                                                                                  >  7:317701/1‑148 (MQ=255)
                                                                    ttCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAAc                                                                              <  4:231434/148‑1 (MQ=255)
                                                                    ttCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAAc                                                                              <  3:309310/148‑1 (MQ=255)
                                                                       cAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAgg                                                                           <  8:182820/148‑1 (MQ=255)
                                                                       cAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGc                                                                          >  6:225653/1‑149 (MQ=255)
                                                                       cAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGc                                                                          <  5:225653/149‑1 (MQ=255)
                                                                        aGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACg                                                                                    <  6:148043/138‑1 (MQ=255)
                                                                        aGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACg                                                                                    >  5:148043/1‑138 (MQ=255)
                                                                        aGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCt                                                                         >  2:232834/1‑149 (MQ=255)
                                                                           aCGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGa                                                                      <  7:265668/149‑1 (MQ=255)
                                                                           aCGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTACCGTTGACGAAGAACAGGCTGGa                                                                      <  6:299206/149‑1 (MQ=255)
                                                                                  cATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCgctg                                                               >  2:218888/1‑149 (MQ=255)
                                                                                    tGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCgctgct                                                             <  4:77958/149‑1 (MQ=255)
                                                                                     gACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGtt                                                                                        >  7:233085/1‑121 (MQ=255)
                                                                                     gACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGtt                                                                                        <  8:233085/121‑1 (MQ=255)
                                                                                       cTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGAc                                                          >  4:14527/1‑149 (MQ=255)
                                                                                          gTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGTCGAAGAACAGGCTGGAGCCGCTGCTGACGCt                                                       >  2:60776/1‑149 (MQ=255)
                                                                                          gTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCt                                                       >  5:234190/1‑149 (MQ=255)
                                                                                          gTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCt                                                       >  1:306985/1‑149 (MQ=255)
                                                                                              cATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCTACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCtggt                                                    <  6:52856/148‑1 (MQ=255)
                                                                                                 cGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCTAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTAGAGCCGCTGCTGACGCTGGTgggg                                                >  2:13662/1‑149 (MQ=255)
                                                                                                  gCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTgggg                                                <  3:14527/148‑1 (MQ=255)
                                                                                                            gCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAAt                                       <  1:238761/147‑1 (MQ=255)
                                                                                                              aaGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATggg                                    >  5:339595/1‑148 (MQ=255)
                                                                                                                  ggCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATg                               <  5:54750/149‑1 (MQ=255)
                                                                                                                  ggCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATg                               <  5:49355/149‑1 (MQ=255)
                                                                                                                       tCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGa                           <  1:232834/148‑1 (MQ=255)
                                                                                                                          ttCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAgg                         <  3:97877/147‑1 (MQ=255)
                                                                                                                          ttCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAgg                         >  4:97877/1‑147 (MQ=255)
                                                                                                                               ccatcATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCt                                                                                                          <  2:114399/61‑1 (MQ=255)
                                                                                                                               ccatcATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCt                                                                                                          >  1:114399/1‑61 (MQ=255)
                                                                                                                                catcatCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAGGGTGAcc                   >  4:140358/1‑147 (MQ=255)
                                                                                                                                  tcatcaGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAGGGTGACCAAg                >  8:143835/1‑148 (MQ=255)
                                                                                                                                      caGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGc                                            >  6:367154/1‑116 (MQ=255)
                                                                                                                                      caGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGc                                            <  5:367154/116‑1 (MQ=255)
                                                                                                                                        gcgcCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTgggg                                                >  3:326239/1‑110 (MQ=255)
                                                                                                                                        gcgcCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTgggg                                                <  4:326239/110‑1 (MQ=255)
                                                                                                                                            cGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAGGGTGACCAAGCAGAACGATgg     <  2:306985/149‑1 (MQ=255)
                                                                                                                                              aGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAGGGTGACCAAGCAGAACGATggca   <  1:218888/149‑1 (MQ=255)
                                                                                                                                              aGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGT‑GGGCCGCACAATGGGCAATGGTGAGGGTGACCAAGCAGAACGATggcag  >  6:191482/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                     
GCCCATTCCCTGTGGGGGCTGAACCGCTCCCACTTCGCTTAATCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAGGGTGACCAAGCAGAACGATGGCAG  >  NC_002947/4557851‑4558142

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: ‑