Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,730,957 A→T 10.2% I128N (ATT→AAT)  hutG ← N‑formylglutamate deformylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,730,9570AT10.2% 179.6 / 8.2 68I128N (ATT→AAT) hutGN‑formylglutamate deformylase
Reads supporting (aligned to +/- strand):  ref base A (33/27);  new base T (5/2);  total (38/30)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

AGGATCGCAACTGGCGCCGTTGAAAGTACCCAGGTTGAAGTCCGGTAACTTGCCGTCGAACAGGTGCGGGATCAGCGAACGGATCGAATGGGCATCCCACAGCAGCGCATAGCCAAACTGCTCACGCAAGCGATCCAGTTCACGACGAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAGCCCGTCTTTGAAGAGGGGTTCGCCGCCGAACAAGGTCGCCGGGTACAAGCCCGTGGTGGCGCCGACGTACAGCGGCTTGTC  >  NC_002947/5730809‑5731105
                                                                                                                                                     |                                                                                                                                                       
aGGATCGCAACTGGCGCCGTTGAAAGTACCCAGGTTGAAGTCCGGTAACTTGCCGTCGAACAGGTGCGGGATCAGCGAACGGATCGAATGGGCATCCCACAGCAGCGCATAGCCAAACTGCTCACGCAAGCGATCCAGTTCACGACGaa                                                                                                                                                         >  7:34098/1‑149 (MQ=255)
    tCGCAACTGGCGCCGTTGAAAGTACCCAGGTTGAAGTCCGGTAACTTGCCGTCGAACAGGTGCGGGATCAGCGAACGGATCGAATGGGCATCCCACAGCAGCGCATAGCCAAACTCCTCACGCAAGCGATCCAGTTCACGACGATTAGt                                                                                                                                                     >  3:16493/1‑149 (MQ=255)
      gCAACTGGCGCCGTTGAAAGTACCCAGGTTGAAGTCCGGTAACTTGCCGTCGAACAGGTGCGGGATCAGCGAACGGATCGAATGGGCATCCCACAGCAGCGCATAGCCAAACTGCTCACGCAAGCGATCCAGTTCACGACGAATAGTc                                                                                                                                                    <  6:271483/148‑1 (MQ=255)
       cAACTGGCGCCGTTGAAAGTACCCAGGTTGAAGTCCGGTAACTTGCCGTCGAACAGGTGCGGGATCAGCGAACGGATCGAATGGGCATCCCACAGCAGCGCATAGCCAAACTGCTCACGCAAGCGATCCAGTTCACGACGAATAGTCgt                                                                                                                                                  <  6:39333/149‑1 (MQ=255)
        aaCTGGCGCCGTTGAAAGTACCCAGGTTGAAGTCCGGTAACTTGCCGTCGAACAGGTGCGGGATCAGCGAACGGATCGAATGGGCATCCCACAGCAGCGCATAGCCAAACTGCTCACGCAAGCGATCCAGTTCACGACTATTCGTCgtg                                                                                                                                                 <  2:10138/149‑1 (MQ=255)
        aaCTGGCGCCGTTGAAAGTACCCAGGTTGAAGTCCGGTAACTTGCCGTCGAACAGGTGCGGGATCAGCGAACGGATCGAATGGGCATCCCACAGCAGCGCATAGCCAAACTGCTCACGCAAGCGATCCAGTTCACGACGAATAGTCgtg                                                                                                                                                 <  3:285397/149‑1 (MQ=255)
         aCTGGCGCCGTTGAAAGTACCCAGGTTGAAGTCCGGTAACTTGCCGTCGAACAGGTGCGGGATCAGCGAACGGATCGAATGGGCATCCCACAGCAGCGCATAGCCAAACTGCTCACGCAAGCGATCCAGTTCACGACGAATAGTCgtgt                                                                                                                                                <  7:290348/149‑1 (MQ=255)
             gcgcCGTTGAAAGTACCCAGGTTGAAGTCCGGTAACTTGCCGTCGAACAGGTGCGGGATCAGCGAACGGATCGAATGGGCATCCCACAGCAGCGCATAGCCAAACTGCTCACGCAAGCGATCCAGTTCACGACGAATAGTCGTGTGAt                                                                                                                                             <  7:246487/148‑1 (MQ=255)
                  gTTGAAAGTACCCAGGTTGAAGTCCGGTAACTTGCCGTCGAACAGGTGCGGGATCAGCGAACGGATCGAATGGGCATCCCACAGCAGCGCATAGCCAAACTGCTCACGCAAGCGATCCAGTTCACGACGAATAGTCGTGTGATAAGGgc                                                                                                                                       <  1:190979/149‑1 (MQ=255)
                            cccAGGTTGAAGTCCGGTAACTTGCCGTCGAACAGGTGCGGGATCAGCGAACGGATCGAATGGGCATCCCACAGCAGCGCATAGCCAAACTGCTCACGCAAGCGATCCAGTTCACGACGAATAGTCGTGTGATAAGGGCGCCAGATTTg                                                                                                                             >  1:133694/1‑149 (MQ=255)
                               aGGTTGAAGTCCGGTAACTTGCCGTCGAACAGGTGCGGGATCAGCGAACGGATCGAATGGGCATCCCACAGCAGCGCATAGCCAAACTGCTCACGCAAGCGATCCAGTTCACGACGAATAGTCGTGTGATAAGGgc                                                                                                                                       <  2:245303/136‑1 (MQ=255)
                               aGGTAGAAGTCCGGTAACTTGCCGTCGAACAGGTGCGGGATCAGCGAACGGATCGAATGGGCATCCCACAGCAGCGCATAGCCAAACTGCTCACGCAAGCGATCCAGTTCACGACGAATAGTCGTGTGATAAGGgc                                                                                                                                       >  1:245303/1‑136 (MQ=255)
                                          cGGTAACTTGCCGTCGAACAGGTGCGGGATCAGCGAACGGATCGAATGGGCATCCCACAGCAGCGCATAGCCAAACTGCTCACGCAAGCGATCCAGTTCACGACGAATAGTCgtgt                                                                                                                                                <  5:237671/116‑1 (MQ=255)
                                          cGGTAACTTGCCGTCGAACAGGTGCGGGATCAGCGAACGGATCGAATGGGCATCCCACAGCAGCGCATAGCCAAACTGCTCACGCAAGCGATCCAGTTCACGACGAATAGTCgtgt                                                                                                                                                >  6:237671/1‑116 (MQ=255)
                                            gTAACTTGCCGTCGAACAGGTGCGGGATCAGCGAACGGATCGAATGGGCATCCCACAGCACCGCATAGCCAAACTGCTCACGCAAGCGATCCAGTTCACGACGAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTg                                                                                                          >  2:307737/1‑149 (MQ=255)
                                             tAACTTGCCGGCGAACAGGTGCGGGATCAGCGAACGGATCGAATGGGCATCCCACAGCAGCGCATAGCCAAACTGCTCACGCAAGCGATCCAGTTCACGACGAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTgc                                                                                                         <  5:256828/149‑1 (MQ=255)
                                                cTTGCCGTCGAACAGCTGCGGGATCAGCGAACGGATCGAATGGGCATCCCACAGCAGCGCATAGCCAAACTGCTCACGCAAGCGATCCAGTTCACGACGAATAGTCGTGTGATAACGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGct                                                                                                      <  8:251626/149‑1 (MQ=255)
                                                    ccGTCGAACAGGTGCGGGATCAGCGAACGGATCGAATGGGCATCCCACAGCAGCGCATAGCCAAACTGCTCACGCAAGCGATCCAGTTCACGACGAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTc                                                                                                  >  2:249843/1‑149 (MQ=255)
                                                      gTCGAACAGGTGCGCGATCAGCGAACGGATCGAATGGGCATCCCACAGCAGCGCATAGCCAAACTGCTCACGCAAGCGATCCAGTTCACGACGAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGc                                                                                                <  2:133694/149‑1 (MQ=255)
                                                           aCAGGTGCGGGATCAGCGAACGGATCGAATGGGCATCCCACAGCAGCGCATAGCCAAACTGCTCACGCAAGCGATCCAGTTCACGACGAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGAc                                                                                           <  5:251889/149‑1 (MQ=255)
                                                               gTGCGGGATCAGCGAACGGATCGAATGGGCATCCCACAGCAGCGCATAGCCAAACTGCTCACGCAAGCGATCCAGTTCACGACGAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGgcg                                                                                       >  3:240184/1‑149 (MQ=255)
                                                               gTGCGGGATCAGCGAACGGATCGAATGGGCATCCCACAGCAGCGCATAGCCAAACTGCTCACGCAAGCGATCCAGTTCACGACGAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGgcg                                                                                       >  7:190326/1‑149 (MQ=255)
                                                                    ggATCAGCGAACGGATCGAATGGGCATCCCACAGCAGCGCATAGCCAAACTGCTCACGCAAGCGATCCAGTTCACGACGAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAGc                                                                                  <  3:211147/149‑1 (MQ=255)
                                                                         aGCGAACGGATCGAATGGGCATCCCACAGCAGCGCATAGCCAAACTGCTCACGCAAGCGATCCAGTTCACGACGAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAGCCCGTc                                                                             >  3:287898/1‑149 (MQ=255)
                                                                          gCGAACGGATCGAATGGGCATCCCACAGCAGCGCATTGCCAAACTGCTCACGCAAGCGATCCAGTTCACGACGAATAGTCGTGTGAAAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAGCCCGTCt                                                                            >  3:279734/1‑149 (MQ=255)
                                                                          gCGAACGGATCGAATGGGCATCCCACAGCAGCGCATAGCCAAACTGCTCACGCAAGCGATCCAGTTCACGACGAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGgcgc                                                                                      >  8:207377/1‑139 (MQ=255)
                                                                          gCGAACGGATCGAATGGGCATCCCACAGCAGCGCATAGCCAAACTGCTCACGCAAGCGATCCAGTTCACGACGAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGgcgc                                                                                      <  7:207377/139‑1 (MQ=255)
                                                                            gAACGGATCGAATGGGCATCCCACAGCAGCGCATAGCCAAACTGCTCACGCAAGCGATCCAGTTCACGACGAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAGCCCGTCttt                                                                          >  4:294830/1‑149 (MQ=255)
                                                                              aCGGATCGAATGGGCATCCCACAGCAGCGCATAGCCAAACTGCTCACGCAAGCGATCCAGTTCACGACGAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAg                                                                                   >  4:30749/1‑138 (MQ=255)
                                                                              aCGGATCGAATGGGCATCCCACAGCAGCGCATAGCCAAACTGCTCACGCAAGCGATCCAGTTCACGACGAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAg                                                                                   <  3:30749/138‑1 (MQ=255)
                                                                               cGGATCGAATGGGCATCCCACAGCAGCGCATAGCCAAACTGCTCACGCAAGCGATCCAGTTCACGACGAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAGCCCGTCTTTGaa                                                                       <  8:34098/149‑1 (MQ=255)
                                                                                ggATCGAATGGGCATCCCACAGCAGCGCATAGCCAAACTGCTCACGCAAGCGATCCAGTTCACGACGAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAGCCCGTCTTTg                                                                         >  1:23562/1‑146 (MQ=255)
                                                                                 gATCGAATGGGCATCCCACAGCAGCGCATAGCCAAACTGCTCACGCAAGCGATCCAGTTCACGACGACTAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAGCCCGTCTTTGAaga                                                                     <  4:16493/149‑1 (MQ=255)
                                                                                  aTCGAATGGGCATCCCACAGCAGCGCATAGCCAAACTGCTCACGCAAGCGATCCAGTTCACGACGAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAGCCCGTCTTTGAagag                                                                    >  3:3055/1‑149 (MQ=255)
                                                                                            cATCCCACAGCAGCGCATAGCCAAACTGCTCACGCAAGCGATCCAGTTCACGACGATTAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTgcg                                                                                                        >  6:99617/1‑103 (MQ=255)
                                                                                            cATCCCACAGCAGCGCATAGCCAAACTGCTCACGCAAGCGATCCAGTTCACGACGATTAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTgcg                                                                                                        <  5:99617/103‑1 (MQ=255)
                                                                                             aTCCCACAGCAGCGCATAGCCAAACTGCTCACGCAAGCGATCCAGTTCACGACGAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAGCCCGTCTTTGAAGAGGGGTTcgccgc                                                         >  1:38476/1‑149 (MQ=255)
                                                                                                 cacaGCAGCGGATAGCCAAACTGCTCACGCAAGCGATCCAGTTCACGACGAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAGCCCGTCTTTGAAGAGGGGTTCGCCGCCGaa                                                     >  6:258764/1‑149 (MQ=255)
                                                                                                   cagcagCGCATAGCCAAACTGCTCACGCAAGCGATCCAGTTCACGACGAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAGCCCGTCTTTGAAGAGGGGTTCGCCGCCGAACa                                                   <  2:23562/149‑1 (MQ=255)
                                                                                                       agagcATAGCCAAACTGCTCACGCAAGCGATCCAGTTCACGACGAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAGCCCGTCTTTGAAGAGGGGTTCGCCGCCGAACAAGGt                                               >  1:212699/4‑149 (MQ=255)
                                                                                                               gCCAAACTTCTCAGGCAAGCGATCCAGTTCAGGACGAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAGCCCGTCTTTGAAGAGGGGTTCTCCGCCGAACAAGGTCGCCGGGt                                       <  4:279734/149‑1 (MQ=255)
                                                                                                                ccAAACTGCTCACGCAAGCGATCCAGTTCACGACGAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGTACCACCGTTTGCGCTCTTCGCCAGACGGCGCCAGCCCGTCTTTGAAGAGGGGTTCGCCGCCGAACAAGGTCGCCgg                                         <  1:249843/149‑1 (MQ=255)
                                                                                                                ccAAACTGCTCACGCAAGCGATCCAGTTCACGACGAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAGCCCGTCTTTGAAGAGGGGTTCGCCGCCGAACAAGGTCGCCGGGTa                                      >  7:231355/1‑149 (MQ=255)
                                                                                                                 cAAACTGCTCACGCAAGCGATCCAGTTCGACGACCAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTCCGATCTTCGCCAGACGGCGCCAGCCCGTCTTTGAAGAGGGGTACGCCGCCGAACAAGGTCGCCGGGt                                       >  1:17545/1‑148 (MQ=255)
                                                                                                                   aaCTGCTCACGCAAGCGATCCAGTTCACGACGAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAGCCCGTCTTTGAAGAGGGGTTCGCCGCCGAACAAGGTCGCCGGGTACaa                                   <  3:294830/149‑1 (MQ=255)
                                                                                                                    aCTGCTCACGCAAGCGATCCAGTTCACGACGAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAGCCCGTCTTTGAAGAGGGGTTCGCCGCCGAACAAGGTCGCCGGGTACAAg                                  >  4:46501/1‑149 (MQ=255)
                                                                                                                     cTGCTCACGCAAGCGATCCAGTTCACGACGAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAGCCCGTCTTTGAAGAGGGGTTCGCCGCCGAACAAGGTCGCCGGGTACAAGc                                 >  7:179861/1‑149 (MQ=255)
                                                                                                                         tCACGCAAGCGATCCAGTTCACGACGAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAGCCCGTCTTTGAAGAGGGGTTCGCCGCCGAACAAGGTCGCCGGGTACAAGCCCgt                             >  5:222900/1‑149 (MQ=255)
                                                                                                                           aCGCAAGCGATCCAGTTCACGACGAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAGCCCGTCTTTGAAGAGGGGTTCGCCGCCGAACAAGGTCGCCGGGTACAAGccc                               >  4:74234/1‑145 (MQ=255)
                                                                                                                           aCGCAAGCGATCCAGTTCACGACGAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAGCCCGTCTTTGAAGAGGGGTTCGCCGCCGAACAAGGTCGCCGGGTACAAGccc                               <  3:74234/145‑1 (MQ=255)
                                                                                                                           aCGCAAGCGATCCAGTTCACGACGAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAGCCCGTCTTTGAAGAGGGGTTCGCCGCCGAACAAGGTCGCCGGGTACAAGCCCgtgg                           >  1:177371/1‑149 (MQ=255)
                                                                                                                            cGCAAGCGATCCAGTTCACGACGAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAGCCCGTCTTTGAAGAGGGGTTCGCCGCCGAACAAGGTCGCCGGGTACAAGCCCgtggt                          >  8:152224/1‑149 (MQ=255)
                                                                                                                               aaGCGATCCAGTTCACGACGAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAGCCCGTCTTTGAAGAGGGGTTCGCCGCCGAACAAGGTCGCCGGGTACAAGCCCGTGGTGgc                       >  6:67460/1‑149 (MQ=255)
                                                                                                                                    aTCCAGTTCACGACGAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAGCCCGTCTTTGAAGAGGGGTTCGCCGCCGAACAAGGTCGCCGGGTACAAGCCCGTGGTGGCGCCGa                  >  7:97556/1‑149 (MQ=255)
                                                                                                                                     tCCAGTTCACGACGAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGA‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAGCCCGTCTTTGAAGAGGGGTTCGCCGCCGAACAAGGTCGCCGGGTACAAGCCCGTGGTGGCGCCGAc                 >  5:36041/1‑149 (MQ=255)
                                                                                                                                       cAGTTCACGACGAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAGCCCGTCTTTGAAGAGGGGTTCGCCGCCGAACAAGGTCGCCGGGTACAAGCCCGTGGTGGCGCCGACGt               <  2:177371/149‑1 (MQ=255)
                                                                                                                                        aGTTCACGACGAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAGCCCGTCTTTGAAGAGGGGTTCGCCGCCGAACAAGGTCGCCGGGTACAAGCCCGTGGTGGCGCCGACGTa              <  2:17545/149‑1 (MQ=255)
                                                                                                                                         gTTCACGACGAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAGCCCGTCTTTGAAGAGGGGTTCGCCGCCGAACAAGGTCGCCGGGTACAAGCCCGTGGTGGCGCCGACGTa              >  5:108987/1‑148 (MQ=255)
                                                                                                                                          ttCACGACGAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAGCCCGTCTTTGAAGAGGGGTTCGCCGCCGAACAAGGTCGCCGGGTACAAGCCCGTGGTGGCGCCGACGTACa            <  2:212699/149‑1 (MQ=255)
                                                                                                                                          cacACGACTATTCGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAGCCCGTCTTTGAAGAGGGGTTCGCCGCCGAACAAGGTCGCCGGGTACAAGCCCGTGGTGGCGCCGACGTa              >  7:275504/3‑147 (MQ=255)
                                                                                                                                          cacACGACTATTCGTCGTGTGATAAGGGCGCCAGATTTGCTACAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAGCCCGTCTTTGAAGAGGGGTTCGCCGCAGAACAAGGTCGCAGCGTACAACCCCgtggt                          >  2:53855/3‑135 (MQ=255)
                                                                                                                                            cACGACTATTCGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAGCCCGTCTTTGAAGAGGGGTTCGCCGCCGAACAAGGTCGCCGGGTACAAGCCCGTGGTGGCGCCGACGTACAGc          >  2:202226/1‑149 (MQ=255)
                                                                                                                                              acgacgAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAGCCCGTCTTTGAAGAGGGGTTCGCCGCCGAACAAGGTCGCCGGGTACAAGc                                 >  2:307837/1‑125 (MQ=255)
                                                                                                                                              acgacgAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAGCCCGTCTTTGAAGAGGGGTTCGCCGCCGAACAAGGTCGCCGGGTACAAGc                                 <  1:307837/125‑1 (MQ=255)
                                                                                                                                              acgacgAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAGCCCGTCTTTGAAGAGGGGTTCGCCGCCGAACAAGGTCGCCGGGTACAAGCCCGTGGTGGCGCCGACGTACAGCg         <  4:287898/149‑1 (MQ=255)
                                                                                                                                                gacgaATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAGCCCGTCTTTGAAGAGGGGTTCGCCGCCGAACAAGGTCGCCGGGTACAAGCCCGTGGTGGCGCCGACGTACAGCGGc       >  2:155086/1‑149 (MQ=255)
                                                                                                                                                gacgaATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAGCCCGTCTTTGAAGAGGGGTTCGCCGCCGAACAAGGTCGCCGGGTACAAGCCCGTGGTGGCGCCGACGTACAGCGGc       >  4:30089/1‑149 (MQ=255)
                                                                                                                                                   gaATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAGCCCGTCTTTGAAGAGGGGTTCGCCGCCGAACAAGGTCGCCGGGTACAAGCCCGTGGTGGCGCCGACGTACAGCGGCTTg    <  6:222900/149‑1 (MQ=255)
                                                                                                                                                   gaATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAGCCCGTCTTTGAAGAGGGGTTCGCCGCCGAACAAGGTCGCCGGGTACAAGCCCGTGGTGGCGCCGACGTACAGCGGCTTg    >  1:223679/1‑149 (MQ=255)
                                                                                                                                                    aaTAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAGCCCGTCTTTGAAGAGGGGTTCGCCGCCGAACAAGGTCGCCGGGTACAAGCCCGTGGTGGCGCCGACGTACAGCGGCTTgt   >  5:142847/1‑149 (MQ=255)
                                                                                                                                                     aTAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAGCCCGTCTTTGAAGAGGGGTTCGCCGCCGAACAAGGTCGCCGGGTACAAGCCCGTGGTGGCGCCGACGTACAGCGGCTTgtc  <  4:3055/149‑1 (MQ=255)
                                                                                                                                                     aTAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAGCCCGTCTTTGAAGAGGGGTTCGCCGCCGAACAAGGTCGCCGGGTACAAGCCCGTGGTGGCGCCGACGTACAGCGGCTTgtc  <  6:108987/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
AGGATCGCAACTGGCGCCGTTGAAAGTACCCAGGTTGAAGTCCGGTAACTTGCCGTCGAACAGGTGCGGGATCAGCGAACGGATCGAATGGGCATCCCACAGCAGCGCATAGCCAAACTGCTCACGCAAGCGATCCAGTTCACGACGAATAGTCGTGTGATAAGGGCGCCAGATTTGCTCCAGGT‑‑‑AGCGTTTGCGCTCTTCGCCAGACGGCGCCAGCCCGTCTTTGAAGAGGGGTTCGCCGCCGAACAAGGTCGCCGGGTACAAGCCCGTGGTGGCGCCGACGTACAGCGGCTTGTC  >  NC_002947/5730809‑5731105

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: