Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 720,922 T→C 9.1% T92A (ACC→GCC) ‡ 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,9220TC 9.1% 157.7 / 4.9 59T92A (ACC→GCC) ‡PP_0614bifunctional N‑carbamoyl‑beta‑alanine amidohydrolase/allantoine amidohydrolase
Reads supporting (aligned to +/- strand):  ref base T (28/24);  new base C (2/3);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.85e-01

GTTGGTCCATACCACCACTTCCAGCGGGGCTTCGGTGGTCAGGTGCAGGTCGTTGAGCGTGCGCATCACTTCGAGGCCGGCCATGACGCCGTAGCAGCCATCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGCTGACCGTGCAGCCGGCCTCCTCGCACCAATGGACGAACAGGTCGCGAGCCTGGCGGTCGAGGTCGGTCAGGGCCAGGCGG  >  NC_002947/720779‑721070
                                                                                                                                               |                                                                                                                                                    
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                                                        <  7:306095/149‑1 (MQ=255)
                                                                                         cGTAGCAGCCATCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGCTGACCGTGCAGCCGGCCTCCTCGCAc                                                        >  5:84997/1‑149 (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      >  5:204102/1‑146 (MQ=255)
                                                                                                                                              gTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGCTGACCGTGCAGCCGGCCTCCTCGCACCAATGGACGAACAGGTCGCGAGCCTGGCGGTCGAGGTCGGTCAGGGCCAg      <  6:204102/146‑1 (MQ=255)
                                                                                                                                               tCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGCTGACCGTGCAGCCGGCCTCCTCGCACCAATGGACGAACAGGTCGCGAGGCTGGCGGTCGAGGTCGGTCAGGGCCAggcgg  <  6:92372/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
GTTGGTCCATACCACCACTTCCAGCGGGGCTTCGGTGGTCAGGTGCAGGTCGTTGAGCGTGCGCATCACTTCGAGGCCGGCCATGACGCCGTAGCAGCCATCGAACTTGCCGCCGGTGGGTTGGGTGTCGATGTGGCTGCCGGTCATGATGGGCGGCAGGGCAGGGTTGCGCCCTGGGCGGCGGGCGAAGATATTGCCGATGGCGTCGACGCTGACCGTGCAGCCGGCCTCCTCGCACCAATGGACGAACAGGTCGCGAGCCTGGCGGTCGAGGTCGGTCAGGGCCAGGCGG  >  NC_002947/720779‑721070

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