Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,622,921 C→G 5.9% A146P (GCC→CCC)  PP_1422 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,622,9210CG 5.9% 218.6 / 3.2 68A146P (GCC→CCC) PP_1422hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (32/32);  new base G (2/2);  total (34/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.71e-01

CTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTCGATCACCAGGCCGATGGC  >  NC_002947/1622775‑1623068
                                                                                                                                                  |                                                                                                                                                   
cTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCt                                                                                                                                                   >  6:362311/1‑149 (MQ=255)
   gAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGt                                                                                                                                                >  8:53869/1‑149 (MQ=255)
   gAGGATCATGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTg                                                                                                                                                  <  5:331663/147‑1 (MQ=255)
            gTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGGCCACCTTGGGTGGGCTGGTcaa                                                                                                                                             >  6:358202/1‑141 (MQ=255)
            gTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGGCCACCTTGGGTGGGCTGGTcaa                                                                                                                                             <  5:358202/143‑3 (MQ=255)
             tCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTga                                                                                                                                       <  1:194366/148‑1 (MQ=255)
               ccccAGCACCGGCTGGAGTTGTTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGACATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGagag                                                                                                                                    <  3:57260/149‑1 (MQ=255)
               ccccAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGt                                                                                                                                           >  7:288854/1‑142 (MQ=255)
               ccccAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGt                                                                                                                                           <  8:288854/142‑1 (MQ=255)
               ccccAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTgaga                                                                                                                                     <  6:214708/148‑1 (MQ=255)
               ccccAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGagag                                                                                                                                    <  2:187278/149‑1 (MQ=255)
                cccAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAgg                                                                                                                                   <  8:49726/149‑1 (MQ=255)
                cccAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAgg                                                                                                                                   >  7:49726/1‑149 (MQ=255)
                    gCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGgtgt                                                                                                                                >  6:104792/1‑148 (MQ=255)
                          gCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGa                                                                                                                         >  2:334929/1‑149 (MQ=255)
                             ggAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACgg                                                                                                                      >  5:261415/1‑149 (MQ=255)
                                 ttGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGa                                                                                                                  <  5:362311/149‑1 (MQ=255)
                                 ttGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGa                                                                                                                  <  4:149094/149‑1 (MQ=255)
                                   gCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCa                                                                                                                    <  6:192763/145‑1 (MQ=255)
                                   gCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCa                                                                                                                    >  5:192763/1‑145 (MQ=255)
                                   gCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTCACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAATCGACGGCAGAtt                                                                                                                >  7:59443/1‑149 (MQ=255)
                                     tCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTgaga                                                                                                                                     >  8:189063/1‑126 (MQ=255)
                                     tCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTgaga                                                                                                                                     <  7:189063/126‑1 (MQ=255)
                                     tCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTc                                                                                                              >  6:10558/1‑149 (MQ=255)
                                       gATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCg                                                                                                             <  5:350724/148‑1 (MQ=255)
                                        aTCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGTGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTa                                                                                                           <  7:16538/149‑1 (MQ=255)
                                        aTCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTa                                                                                                           <  5:104792/149‑1 (MQ=255)
                                         tCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGt                                                                                                                                           >  5:1945/1‑116 (MQ=255)
                                         tCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGt                                                                                                                                           <  6:1945/116‑1 (MQ=255)
                                           acaTAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCt                                                                                                        <  5:10558/149‑1 (MQ=255)
                                                      gCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCa                                                                                             >  4:191337/1‑149 (MQ=255)
                                                      gCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCa                                                                                             >  5:351193/1‑149 (MQ=255)
                                                        aaaCCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACt                                                                                           >  8:205523/1‑149 (MQ=255)
                                                             gggTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAg                                                                                      <  6:351193/149‑1 (MQ=255)
                                                                    acaacaGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTg                                                                               >  4:146077/1‑149 (MQ=255)
                                                                      aacaGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGAt                                                                                                     <  1:102514/125‑1 (MQ=255)
                                                                      aacaGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGAt                                                                                                     >  2:98523/1‑125 (MQ=255)
                                                                      aacaGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGAt                                                                                                     <  1:98523/125‑1 (MQ=255)
                                                                      aacaGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGAt                                                                                                     >  2:102514/1‑125 (MQ=255)
                                                                      aacaGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGc                                                                             >  4:286549/1‑149 (MQ=255)
                                                                      aacaGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGc                                                                             >  8:196437/1‑149 (MQ=255)
                                                                        caGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGagag                                                                                                                                    <  8:272923/92‑1 (MQ=255)
                                                                        caGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGagag                                                                                                                                    >  7:272923/1‑92 (MQ=255)
                                                                         aGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCa                                                                          <  1:334929/149‑1 (MQ=255)
                                                                          gTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGcc                                                                            >  1:174590/1‑146 (MQ=255)
                                                                           ttCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGcc                                                                            >  6:115934/1‑145 (MQ=255)
                                                                                    ctctTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCa                                                               <  3:191337/149‑1 (MQ=255)
                                                                                       ttCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGtt                                                            >  5:130761/1‑149 (MQ=255)
                                                                                         cGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTCACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGttt                                                           <  8:59443/148‑1 (MQ=255)
                                                                                              tAGATGAGTAGCGTCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTc                                                     <  6:261415/149‑1 (MQ=255)
                                                                                                      tACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTg                                                                               <  8:56868/115‑1 (MQ=255)
                                                                                                      tACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTg                                                                               >  7:56868/1‑115 (MQ=255)
                                                                                                         cGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGc                                           >  4:77354/1‑148 (MQ=255)
                                                                                                           gCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTg                                                        >  7:266200/1‑133 (MQ=255)
                                                                                                           gCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTg                                                        <  8:266200/133‑1 (MQ=255)
                                                                                                              gcgcgcCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCgg                                     <  3:77354/149‑1 (MQ=255)
                                                                                                               cgcgcCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCggg                                    <  7:205523/149‑1 (MQ=255)
                                                                                                                   ccAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATg                                <  6:130761/149‑1 (MQ=255)
                                                                                                                       tCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGaa                             >  2:46850/1‑148 (MQ=255)
                                                                                                                        cTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACTGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGaaaa                           >  1:123515/1‑149 (MQ=255)
                                                                                                                           cccgACCAGGCCACCTTGGGTGGGCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCAt                                 >  4:151837/4‑140 (MQ=255)
                                                                                                                           cccgACCAGGCCACCTTGGGTGGGCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCAt                                 <  3:151837/137‑1 (MQ=255)
                                                                                                                              gACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGc                      <  3:286549/148‑1 (MQ=255)
                                                                                                                              gACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGc                      <  7:196437/148‑1 (MQ=255)
                                                                                                                                 cAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTCGa                  >  5:113706/1‑149 (MQ=255)
                                                                                                                                      cacccaAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGc                      >  2:277749/1‑140 (MQ=255)
                                                                                                                                      cacccaAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGc                      <  1:277749/140‑1 (MQ=255)
                                                                                                                                       acccaAGGTGGCATCCTCGGGTTTCACAGGTGAGAAGCGACGGCACATTTCGTAGCTGATGCCCGCCACTACCAGGCTACTGGCCCAGATAACGCCCAGTTTCTGCTCCACTGCCCGCGAGCGCGCATGGAAAATCTGCTCCATCAc              >  8:313293/1‑147 (MQ=255)
                                                                                                                                                  ccTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTCGATCACCAGGCCGATGGc  >  3:208629/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTCGATCACCAGGCCGATGGC  >  NC_002947/1622775‑1623068

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