Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,557,984 T→C 6.3% M58V (ATG→GTG)  PP_4045 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,557,9840TC 6.3% 185.6 / 2.0 63M58V (ATG→GTG) PP_4045hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (29/30);  new base C (2/2);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

TCAGCCTCTGATCTTGCCCATTCCCTGTGGGGGCTGAACCGCTCCCACTTCGCTTAATCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAGGGTGACCAAG  >  NC_002947/4557836‑4558128
                                                                                                                                                    |                                                                                                                                                
tCAGCCTCTGATCTTGCCCATTCCCTGTGGGGGCTGAACCGCTCCCACTTCGCTTAATCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCcatcat                                                                                                                                                  <  5:331989/149‑1 (MQ=255)
    ccTCTGATCTTGCCCATTCCCTGTGGGGGCTGAACCGCTCCCACTTCGCTTAATCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCcatcat                                                                                                                                                  <  5:159707/145‑1 (MQ=255)
    ccTCTGATCTTGCCCATTCCCTGTGGGGGCTGAACCGCTCCCACTTCGCTTAATCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCcatcat                                                                                                                                                  >  6:159707/1‑145 (MQ=255)
    ccTCTGATCTTGCCCATTCCCTGTGGGGGCTGAACCGCTCCCACTTCGCTTAATCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAgc                                                                                                                                              <  7:123259/149‑1 (MQ=255)
     ctctGATCTTGCCCATTCCCTGTGGGGGCTGAACCGCTCCCACTTCGCTTAATCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAgcg                                                                                                                                             >  1:238828/1‑149 (MQ=255)
      tctGATCTTGCCCATTCCCTGTGGGGGCTGAACCGCTCCCACTTCGCTTAATCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCACCAgcgg                                                                                                                                            >  1:170759/1‑148 (MQ=255)
               gCCCATTCCCTGTGGGGGCTGAACCGCTCCCACTTCGCTTAATCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCt                                                                                                                                   >  2:171285/1‑149 (MQ=255)
                cccATTCCCTGTGGGGGCTGAACCGCTCCCACTTCGCTTAATCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTg                                                                                                                                  >  3:231434/1‑149 (MQ=255)
                 ccATTCCCTGTGGGGGCTGAACCGCTCCCACTTCGCTTAATCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTgg                                                                                                                                 <  6:175398/149‑1 (MQ=255)
                 ccATTCCCTGTGGGGGCTGAACCGCTCCCACTTCGCTTAATCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTgg                                                                                                                                 <  8:127058/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:49355/1‑149 (MQ=255)
                                                                   aacTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTg                                                                               >  6:54750/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                                                                                     >  3:338537/1‑135 (MQ=255)
                                                                           tCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGgcc                                                                                     <  4:338537/135‑1 (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                                                            <  5:225653/149‑1 (MQ=255)
                                                                                      cAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGc                                                            >  6:225653/1‑149 (MQ=255)
                                                                                       aGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACg                                                                      >  5:148043/1‑138 (MQ=255)
                                                                                       aGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACg                                                                      <  6:148043/138‑1 (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                                                                          <  8:233085/121‑1 (MQ=255)
                                                                                                    gACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGtt                                                                          >  7:233085/1‑121 (MQ=255)
                                                                                                      cTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGAc                                            >  4:14527/1‑149 (MQ=255)
                                                                                                         gTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGTCGAAGAACAGGCTGGAGCCGCTGCTGACGCt                                         >  2:60776/1‑149 (MQ=255)
                                                                                                         gTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCt                                         >  1:306985/1‑149 (MQ=255)
                                                                                                         gTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCt                                         >  5:234190/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           >  4:97877/1‑147 (MQ=255)
                                                                                                                                         ttCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAgg           <  3:97877/147‑1 (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)
                                                                                                                                                    |                                                                                                                                                
TCAGCCTCTGATCTTGCCCATTCCCTGTGGGGGCTGAACCGCTCCCACTTCGCTTAATCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAGGGTGACCAAG  >  NC_002947/4557836‑4558128

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