Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,622,915 A→G 6.9% L148L (TTG→CTG) ‡ PP_1422 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,622,9150AG 6.9% 211.5 / 4.7 72L148L (TTG→CTG) ‡PP_1422hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (32/35);  new base G (3/2);  total (35/37)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.88e-01

CAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTCGATCACCAGG  >  NC_002947/1622770‑1623060
                                                                                                                                                  |                                                                                                                                                  
cAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGt                                                                                                                                                 >  4:18685/1‑149 (MQ=255)
 aGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTg                                                                                                                                                <  3:167744/149‑1 (MQ=255)
   gCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAACGTGGc                                                                                                                                              <  7:286539/149‑1 (MQ=255)
     cTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCt                                                                                                                                            <  8:299979/149‑1 (MQ=255)
        gAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGt                                                                                                                                         >  2:33220/1‑149 (MQ=255)
                tGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGcc                                                                                                                                             >  5:69910/1‑137 (MQ=255)
                tGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGcc                                                                                                                                             <  6:69910/137‑1 (MQ=255)
                 gTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGGCCACCTTGGGTGGGCTGGTcaacc                                                                                                                                    <  5:169254/145‑5 (MQ=255)
                  tCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTga                                                                                                                                <  3:18685/148‑1 (MQ=255)
                  tCCCCCAGCACCGGCGGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCTGCATAGATGAGTACCGGCAGCGCGCGAATCTGTTTGACCAGCCCACGCAAGGTGGCCTGGTCGGGTTTga                                                                                                                                <  8:229294/148‑1 (MQ=255)
                    ccccAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATATATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGagag                                                                                                                             >  8:165073/1‑149 (MQ=255)
                    ccccAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGagag                                                                                                                             <  1:311532/149‑1 (MQ=255)
                             cGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGc                                                                                                                   >  6:293328/1‑149 (MQ=255)
                              ggCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCACTCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGc                                                                                                                   <  1:253273/149‑1 (MQ=255)
                               gCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGTAGCGa                                                                                                                 >  4:84708/1‑149 (MQ=255)
                                         cTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTg                                                                                                                                 >  2:297642/1‑124 (MQ=255)
                                         cTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTg                                                                                                                                 <  1:297642/124‑1 (MQ=255)
                                         cTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGAttt                                                                                                       <  5:301701/149‑1 (MQ=255)
                                            gATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCg                                                                                                     >  6:291121/1‑148 (MQ=255)
                                              tCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAg                                                                                                  <  3:210341/149‑1 (MQ=255)
                                               cacaTAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGc                                                                                                 >  3:244302/1‑149 (MQ=255)
                                                     aTGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGc                                                                                           >  5:121884/1‑149 (MQ=255)
                                                         aGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGTCAGATTTCGTAGCTGATGCCCGc                                                                                       <  6:121884/149‑1 (MQ=255)
                                                         aGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGc                                                                                       >  8:141675/1‑149 (MQ=255)
                                                         aGGCAAACCGGGTCAGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGc                                                                                       <  7:63485/149‑1 (MQ=255)
                                                           gCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCa                                                                                     >  3:83263/1‑149 (MQ=255)
                                                           gCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCa                                                                                     >  6:20878/1‑149 (MQ=255)
                                                                  gggTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAg                                                                              <  7:141675/149‑1 (MQ=255)
                                                                  gggTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAg                                                                              <  5:20878/149‑1 (MQ=255)
                                                                    gTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTCGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAgg                                                                             <  1:60596/149‑1 (MQ=255)
                                                                         acaacaGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGATGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTg                                                                       <  4:244302/149‑1 (MQ=255)
                                                                           aacaGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGc                                                                     >  5:123397/1‑149 (MQ=255)
                                                                             caGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGcc                                                                    >  2:59307/1‑148 (MQ=255)
                                                                              aGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCa                                                                  <  4:83263/149‑1 (MQ=255)
                                                                               gTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCTGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGcc                                                                    >  4:3505/1‑146 (MQ=255)
                                                                               gTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCTGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGcc                                                                    <  3:3505/146‑1 (MQ=255)
                                                                               gTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGAt                                                                                             >  1:118494/1‑121 (MQ=255)
                                                                               gTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGAt                                                                                             <  2:118494/121‑1 (MQ=255)
                                                                                  cGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGata                                                              >  1:8547/1‑149 (MQ=255)
                                                                                        gCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGc                                                                                       >  3:302477/1‑118 (MQ=255)
                                                                                        gCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGc                                                                                       <  4:302477/118‑1 (MQ=255)
                                                                                          tctTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAg                                                      >  5:28248/1‑149 (MQ=255)
                                                                                            ttCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACtt                                                                                  >  5:171623/1‑119 (MQ=255)
                                                                                            ttCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACtt                                                                                  <  6:171623/119‑1 (MQ=255)
                                                                                               ggCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTc                                                  <  1:59307/148‑1 (MQ=255)
                                                                                               ggCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTc                                                  <  5:293328/148‑1 (MQ=255)
                                                                                                 cATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGc                                               <  2:8547/149‑1 (MQ=255)
                                                                                                  aTAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAg                                                                                                  <  2:264856/97‑1 (MQ=255)
                                                                                                  aTAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAg                                                                                                  >  1:264856/1‑97 (MQ=255)
                                                                                                   tAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTc                                             <  6:123397/149‑1 (MQ=255)
                                                                                                   tAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTc                                             >  6:268139/1‑149 (MQ=255)
                                                                                                      aTGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATcgc                                                          <  8:153442/133‑1 (MQ=255)
                                                                                                      aTGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATcgc                                                          >  7:153442/1‑133 (MQ=255)
                                                                                                       tGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGcc                                                                    >  5:201037/1‑122 (MQ=255)
                                                                                                       tGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGcc                                                                    <  6:201037/122‑1 (MQ=255)
                                                                                                       ggAGGACCGGCAGCGCGCCACTCTGGGTGACCAGCCCACCCAAGGGGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATCTCGCGCAGTTTCTGCTCCACt                                         <  7:262845/148‑1 (MQ=255)
                                                                                                        gAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTg                                        >  8:139300/1‑149 (MQ=255)
                                                                                                             ccGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGc                                   >  2:238266/1‑149 (MQ=255)
                                                                                                              cGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGc                                   >  2:268945/1‑148 (MQ=255)
                                                                                                                        gcCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCAt                         <  6:28248/149‑1 (MQ=255)
                                                                                                                             tCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCa                                                                               <  3:198089/90‑1 (MQ=255)
                                                                                                                             tCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCa                                                                               >  4:198089/1‑90 (MQ=255)
                                                                                                                               aacccgACCAGGCCACCTTGGGTGGGCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAAt                  >  1:153450/6‑149 (MQ=255)
                                                                                                                                acccgACCAGGCCACCTTGGGTGGGCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGc                                   >  5:235953/5‑131 (MQ=255)
                                                                                                                                acccgACCAGGCCACCTTGGGTGGGCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGc                                   <  6:235953/127‑1 (MQ=255)
                                                                                                                                    gACCAGGCCACCTTGGGTGGGCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCt             >  1:144135/1‑149 (MQ=255)
                                                                                                                                     aCCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCTGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTc            <  7:165073/149‑1 (MQ=255)
                                                                                                                                         gCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTCGATc        <  1:33220/149‑1 (MQ=255)
                                                                                                                                         gCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACGTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTCGAt         >  2:239305/1‑148 (MQ=255)
                                                                                                                                             acccaAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTCGATc        >  3:275240/1‑145 (MQ=255)
                                                                                                                                             acccaAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTCGATc        <  4:275240/145‑1 (MQ=255)
                                                                                                                                               ccaAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTCGATCACCAgg  >  1:303403/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTCGATCACCAGG  >  NC_002947/1622770‑1623060

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: