Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 720,926 G→C 9.4% I90M (ATC→ATG)  PP_0614 ← bifunctional N‑carbamoyl‑beta‑alanine amidohydrolase/allantoine amidohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947720,9260GC 9.4% 174.3 / 7.8 59I90M (ATC→ATG) PP_0614bifunctional N‑carbamoyl‑beta‑alanine amidohydrolase/allantoine amidohydrolase
Reads supporting (aligned to +/- strand):  ref base G (27/26);  new base C (2/3);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.55e-01

GTTGGTCCATACCACCACTTCCAGCGGGGCTTCGGTGGTCAGGTGCAGGTCGTTGAGCGTGCGCATCACTTCGAGGCCGGCCATGACGCCGTAGCAGCCATCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGCTGACCGTGCAGCCGGCCTCCTCGCACCAATGGACGAACAGGTCGCGAGCCTGGCGGTCGAGGTCGGTCAGGGCCAGGCGGCA  >  NC_002947/720779‑721072
                                                                                                                                                   |                                                                                                                                                  
gTTGGTCCATACCACCACTTCCAGCGGGGCTTCGGTGGTCAGGTGCAGGTCGTTGAGCGTGCGCATCACTTCGAGGCCGGCCATGACGCCGTAGCAGCCATCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGGTCatga                                                                                                                                                   >  1:126082/1‑149 (MQ=255)
    gTCCATACCACCACTTCCAGCGGGGCTTCGGTGGTCAGGTGCAGGTCGTTGAGCGTGCGCATCACTTCGAGGCCGGCCATGACGCCGTAGCAGCCATCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGGTCATGATggg                                                                                                                                               <  7:89327/149‑1 (MQ=255)
           ccaccaCTTCCAGCGGGGCTTCGGTGGTCAGGTGCAGGTCGTTGAGCGTGCGCATCACTTCGAGGCCGGCCATGACGCCGTAGCAGCCATCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGTCGCCCATCATGACCggcagc                                                                                                                                        <  2:269768/149‑2 (MQ=25)
           ccaccaCTTCCAGCGGGGCTTCGGTGGTCAGGTGCAGGTCGTTGAGCGTGCGCATCACTTCGAGGCCGGCCATGACGCCGTAGCAGCCATCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGCCCATCATGACCggcagc                                                                                                                                        <  2:261691/149‑2 (MQ=25)
            caccacTTCCAGCGGGGCTTCGGTGGTCAGGTGCAGGTCGTTGAGCGTGCGCATCACTTCGAGGCCGGCCATGACGCCGTAGCAGCCATCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGGTCATGATGGGCggcaggg                                                                                                                                       >  8:189838/1‑149 (MQ=255)
                        cGGGGCTTCGGTGGTCAGGTGCAGGTCGTTGAGCGTGCGCATCACTTCGAGGCCGGCCATGACGCCGTAGCAGCCATCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGcc                                                                                                                           >  1:212837/1‑149 (MQ=255)
                             cTTCGGTGGTCAGGTGCAGGTCGTTGAGCGTGCGCATCACTTCGAGGCCGGCCATGACGCCGTAGCAGCCATCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGCCCATCATGACCggcag                                                                                                                                         >  7:106228/1‑130 (MQ=25)
                             cTTCGGTGGTCAGGTGCAGGTCGTTGAGCGTGCGCATCACTTCGAGGCCGGCCATGACGCCGTAGCAGCCATCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGCCCATCATGACCggcag                                                                                                                                         <  8:106228/130‑1 (MQ=25)
                                 ggtggtCAGGTGCAGGTCGTTGAGCGTGCGCATCACTTCGAGGCCGGCCATGACGCCGTAGCAGCCATCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGggcggc                                                                                                                  >  8:26453/1‑149 (MQ=255)
                                     gtCAGGTGCAGGTCGTTGAGCGTGCGCATCACTTCGAGGCCGGCCATGACGCCGTAGCAGCCATCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCggcggg                                                                                                               <  7:26453/148‑1 (MQ=255)
                                         ggTGCAGGTCGTTGAGCGTGCGCATCACTTCGAGGCCGGCCAGGACGCCGTAGCAGCCATCGAACTTGCCGCCGGTGGGGTGGGTGTCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAg                                                                                                          <  8:20828/149‑1 (MQ=255)
                                          gTGCAGGTCGTTGAGCGTGCGCATCACTTCGAGGCCGGCCATGACGCCGTAGCAGCCATCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCggcgggcg                                                                                                             >  2:250398/1‑145 (MQ=255)
                                          gTGCAGGTCGTTGAGCGTGCGCATCACTTCGAGGCCGGCCATGACGCCGTAGCAGCCATCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGa                                                                                                         <  7:189838/149‑1 (MQ=255)
                                              aGGTCGTTGAGCGTGCGCATCACTTCGAGGCCGGCCATGACGCCGTAGCAGCCATCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGcggcgg                                                                                                                >  8:258827/1‑138 (MQ=255)
                                              aGGTCGTTGAGCGTGCGCATCACTTCGAGGCCGGCCATGACGCCGTAGCAGCCATCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGcggcgg                                                                                                                <  7:258827/138‑1 (MQ=255)
                                               ggTCGTTGAGCGTGCGCATCACTTCGAGGCCGGCCATGACGCCGTAGCAGCCATCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTg                                                                                                    >  8:306095/1‑149 (MQ=255)
                                                gTCGTTGAGCGTGCGCATCACTTCGAGGCCGGCCATGACGCCGTAGCAGCCATCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGc                                                                                                   <  1:250398/149‑1 (MQ=255)
                                                   gTTGAGCGTGCGCATCACTTCGAGGCCGGCCATGACGCCGTAGCAGCCATCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGa                                                                                                <  2:212837/149‑1 (MQ=255)
                                                      gAGCGTGCGCATCACTTCGAGGCCGGCCATGACGCCGTAGCAGCCATCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATg                                                                                              <  2:126082/148‑1 (MQ=255)
                                                            gcgcATCACTTCGAGGCCGGCCATGACGCCGTAGCAGCCATCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGa                                                                                       >  7:63738/1‑149 (MQ=255)
                                                               cATCACTTCGAGGCCGACCATGACGCCGTAGCAGCCATCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGc                                                                                    >  3:123521/1‑149 (MQ=255)
                                                               cATCACTTCGAGGCCGACCATGACGCCGTAGCAGCCATCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGc                                                                                    <  4:123521/149‑1 (MQ=255)
                                                                    cTTCGAGGCCGGCCATGACGCCGTCGCAGCCATCGAACGTGCCGCCGGGGGGTAGGGTGTCGAGGTGGCTGCCGGGCCTGATGGGCGGCAGGGCAGGGTGGCGCCCTGGGCGGCGGGCGAAGAGATTGCCGCGGGCGTCGCCGCTGAcc                                                                               >  1:263005/1‑149 (MQ=255)
                                                                     ttCGAGGCCGGCCATGACGCCGTAGCAGCCATCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGCTGAcag                                                                              >  1:322460/1‑147 (MQ=255)
                                                                         aGGCCGGCCATGACGCCGTAGCAGCCATCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGCTGACCGTg                                                                            >  8:67726/1‑147 (MQ=255)
                                                                          ggccggccATGACGCCGTAGCAGCCATCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGCTGACCGTGCAg                                                                         >  3:299576/1‑149 (MQ=255)
                                                                             cggccATGACGCCGTAGCAGCCATCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGCTGACCGTGCAGCCg                                                                      >  1:1549/1‑149 (MQ=255)
                                                                             cggccATGACGCCGTAGCAGCCATCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGCTGACCGTGCAGCCg                                                                      >  8:324686/1‑149 (MQ=255)
                                                                                    gACGCCGTAGCAGCCATCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGCTGAGCGTGCAGCCGGcctcct                                                               >  2:312212/1‑149 (MQ=255)
                                                                                      cGCCGTAGCAGCCATCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGCTGACCGTGCAGCCGGCCTCCTCg                                                             >  4:281699/1‑149 (MQ=255)
                                                                                       gCCGTAGCAGCCATCTAACTTGCCGCCGGTGGGTTGGGTGTCGTTGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGCTGACCGTGCAGCCGGcctcct                                                               <  6:294396/146‑1 (MQ=255)
                                                                                       gCCGTAGCAGCCATCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGCTGACCGTGCAGCCGGCCTCCTCg                                                             >  5:294396/1‑148 (MQ=255)
                                                                                         cGTAGCAGCCATCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGCTGACCGTGCAGCCGGCCTCCTCGCAc                                                          >  5:84997/1‑149 (MQ=255)
                                                                                         cGTAGCAGCCATCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGCTGACCGTGCAGCCGGCCTCCTCGCAc                                                          <  7:306095/149‑1 (MQ=255)
                                                                                             gcagcCATCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGCTGACCGTGCAGCCGGCCTCCTCGCACCaa                                                       <  3:141737/148‑1 (MQ=255)
                                                                                             gcagcCATCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGCTGACCGTGCAGCCGGCCTCCTCGCACCAAt                                                      <  6:84997/149‑1 (MQ=255)
                                                                                               agcCATCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGCTGACCGTGCAGCCGGCCTCCTCGCACCAATgg                                                    <  4:299576/149‑1 (MQ=255)
                                                                                                gcCATCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGCTGACCGTGCAGCCGGCCTCCTCGCACCAATGGa                                                   >  2:283692/1‑149 (MQ=255)
                                                                                                gcCATCGAACTTGCCGCCGGGGGGTTGGGTGTCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGCTGACCGTGCAGCCGGCCTCCTCGCACCAATGGa                                                   <  7:150708/149‑1 (MQ=255)
                                                                                                 ccttCGAACTTGCCGCCGGTGGGTTGGGTTTCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCTGATGGCGTCGACGCTGACCGTGCAGCTGGCCTCCTCGCAc                                                          <  8:52370/138‑1 (MQ=255)
                                                                                                 ccATCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGCTGACCGTGCAGCCGGCCTCCTCGCACCAATGGAc                                                  >  7:256867/1‑149 (MQ=255)
                                                                                                 ccATCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGCTGACCGTGCAGCCGGCCTCCTCGCACCAATGGAc                                                  <  7:324686/149‑1 (MQ=255)
                                                                                                 ccATCGAACTTGCCGCCGGTGGGTTGGGTGTAGATGTGGCTGCCGGTAATAATGGGCGACAGGACAGGGTTGCGCACTGGGCGGAGGGCGAAGATATTGACTATGGAGTCCACGCTGACAGTGCAGCCGGCCTCCTCGCACCatatcg                                                   >  7:52370/1‑143 (MQ=21)
                                                                                                   aTCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGCTGACCGTGCAGCCGGCCTCCTCGCACCAATGGACg                                                 >  7:284480/1‑148 (MQ=255)
                                                                                                       aaCTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGCTGACCGTg                                                                            <  2:173694/117‑1 (MQ=255)
                                                                                                       aaCTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGCTGACCGTGCAGCCGGCCTCCTCGCACCAATGGACGAACAg                                            >  5:92372/1‑149 (MQ=255)
                                                                                                                   gTGGGTTGGGTGTCTATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGCTGACCGTGCAGCCGGCCTCCTCGCACCAATGGACGAACAGGTCGCGAGCCTg                                <  8:284480/149‑1 (MQ=255)
                                                                                                                    ggggtggggtgTGGCTGTGGCTGCCGGGCAGGATGGGCGGCAGGGCAGGGGTGCGCCCTGGGCGGCGGGCGAAGATATGGCCGATGGCGTCGACGCTGCCCGTGCAGCCGGCCTCCTCGCACCAATGGACGAACAGGTCGCGAGCCTgg                               <  2:263005/148‑1 (MQ=255)
                                                                                                                     gggttgggtGTCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGCTGACCGTGCAGCCGGCCTCCTCGCACCAATGGACGAACAGCTCGCGAGCCTGGc                              >  6:137170/1‑149 (MQ=255)
                                                                                                                             tgtCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGCTGACCGTGCAGCCGGCCTCCTCGCACCAATGGACGAACAGGTCGCGAGCCTGGCGGTCGAgg                      >  3:75999/1‑149 (MQ=255)
                                                                                                                               tCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGCTGACCGTGCAGCCGGCCTCCTCGCACCAATGGACGAACAGGTCGCGAGCCTGGCGGTCGAggtc                    <  7:67726/149‑1 (MQ=255)
                                                                                                                                     tgGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGCTGACCGTGCAGCCGGCCTCGTCGCTCCAATGGACGAAGAGGTCGCGAGCCTGGCGTTCGAGGTCGGTCAg              <  5:137170/149‑1 (MQ=255)
                                                                                                                                     gcccTGCCGCCCATCATGACCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGCTGACCGTGCAGCCGGCCTCCTCGCACCAATGGACGAACAGGTCGCGAGCCTGGCGGTCGAGGTCGGTCAg              >  3:32920/4‑149 (MQ=21)
                                                                                                                                       gCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGCTGACCGTGCAGCCGGCCTCCTCGCACCAATGGACGAACAGGTCGCGAGCCTGGCGGTCGAGGTCGGTCAggg            >  1:66943/1‑149 (MQ=255)
                                                                                                                                       gCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGCTGACCGTGCAGCCGGCCTCCTCGCACCAATGGACGAACAGGTCGCGAGCCTGGCGGTCGAGGTCGGTCAgg             <  3:281699/148‑1 (MQ=255)
                                                                                                                                           ccGGTCATGATGGGCGGCATGTCAGGGTTGCGCCCTGGGCGGCGGGCTAAGATATTGCCGATGGCGTCGACGCTGACCGTGCAGCCGGCCTCCTCGCACCAATGGACGAACAGGTCGCGAGCCTGGCGGTCGAGGTCGGTCAGGGCCAg        <  2:322460/149‑1 (MQ=255)
                                                                                                                                              gTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGCTGACCGTGCAGCCGGCCTCCTCGCACCAATGGACGAACAGGTCGCGAGCCTGGCGGTCGAGGTCGGTCAGGGCCAggcg     >  2:24231/1‑149 (MQ=255)
                                                                                                                                              gTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGCTGACCGTGCAGCCGGCCTCCTCGCACCAATGGACGAACAGGTCGCGAGCCTGGCGGTCGAGGTCGGTCAGGGCCAg        <  6:204102/146‑1 (MQ=255)
                                                                                                                                              gTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGCTGACCGTGCAGCCGGCCTCCTCGCACCAATGGACGAACAGGTCGCGAGCCTGGCGGTCGAGGTCGGTCAGGGCCAg        >  5:204102/1‑146 (MQ=255)
                                                                                                                                               tCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGCTGACCGTGCAGCCGGCCTCCTCGCACCAATGGACGAACAGGTCGCGAGGCTGGCGGTCGAGGTCGGTCAGGGCCAggcgg    <  6:92372/149‑1 (MQ=255)
                                                                                                                                                  tgatgGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGCTGACCGTGCAGCCGGCCTCCTCGCACCAATGGACGAACAGGTCGCGAGCCTGGCGGTCGAGGTCGGTCAGGGCCAGGCGGCa  >  4:15592/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GTTGGTCCATACCACCACTTCCAGCGGGGCTTCGGTGGTCAGGTGCAGGTCGTTGAGCGTGCGCATCACTTCGAGGCCGGCCATGACGCCGTAGCAGCCATCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGCTGACCGTGCAGCCGGCCTCCTCGCACCAATGGACGAACAGGTCGCGAGCCTGGCGGTCGAGGTCGGTCAGGGCCAGGCGGCA  >  NC_002947/720779‑721072

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