Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 720,930 G→A 8.4% P89L (CCC→CTC) ‡ 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,9300GA 8.4% 167.8 / 4.2 59P89L (CCC→CTC) ‡PP_0614bifunctional N‑carbamoyl‑beta‑alanine amidohydrolase/allantoine amidohydrolase
Reads supporting (aligned to +/- strand):  ref base G (29/25);  new base A (2/3);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

GTCCATACCACCACTTCCAGCGGGGCTTCGGTGGTCAGGTGCAGGTCGTTGAGCGTGCGCATCACTTCGAGGCCGGCCATGACGCCGTAGCAGCCATCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGCTGACCGTGCAGCCGGCCTCCTCGCACCAATGGACGAACAGGTCGCGAGCCTGGCGGTCGAGGTCGGTCAGGGCCAGGCGGCATACGCC  >  NC_002947/720783‑721078
                                                                                                                                                   |                                                                                                                                                    
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)
                                                                                                                                                   ggcggcAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGCTGACCGTGCAGCCGGCCTCCTCGCACCAATGGACGAACAGGTCGCGAGCCTGGCGGTCGAGGTCGGTCAGGGCCAGGCGGCATAcgcc  <  4:75999/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GTCCATACCACCACTTCCAGCGGGGCTTCGGTGGTCAGGTGCAGGTCGTTGAGCGTGCGCATCACTTCGAGGCCGGCCATGACGCCGTAGCAGCCATCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGCTGACCGTGCAGCCGGCCTCCTCGCACCAATGGACGAACAGGTCGCGAGCCTGGCGGTCGAGGTCGGTCAGGGCCAGGCGGCATACGCC  >  NC_002947/720783‑721078

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