Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,130,113 G→C 8.5% A107G (GCC→GGC)  PP_5378 ← metal‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,130,1130GC 8.5% 166.5 / 6.2 59A107G (GCC→GGC) PP_5378metal‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (26/28);  new base C (3/2);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

ATCTGCGGCATTTGCTCAAGGAAAGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATCCTTTAGGTTCGGAACCTCGAT  >  NC_002947/6129966‑6130258
                                                                                                                                                   |                                                                                                                                                 
agatGCGGCATTTGCTCAAGGAAAGCGACAATACCCAATATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGTTAGCCTTGATACCGTGGTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTgg                                                                                                                                                   <  5:201936/145‑1 (MQ=255)
   tGCGGCATTTGCTCAAGGAAAGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGTGTCAACGCTTAACCCAATCTTGGc                                                                                                                                                  >  4:152748/1‑146 (MQ=255)
     cGGCATTTGCTCAAGGAAAGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCaa                                                                                                                                                <  6:225277/146‑1 (MQ=255)
         aTTTGGTCAAGGAAAGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTggg                                                                                                                                         >  1:50036/1‑149 (MQ=255)
         aTTTGCTCAAGGAAAGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGCCAAGATTggt                                                                                                                                         >  3:129835/1‑148 (MQ=255)
             gCTCAAGGAAAGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGt                                                                                                                                     >  1:207537/1‑149 (MQ=255)
              cTCAAGGAAAGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGCCAAGATTggtttaag                                                                                                                                    <  5:242163/149‑7 (MQ=255)
              cTCAAGGAAAGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGCCAAGATTggtttaa                                                                                                                                     <  4:129835/148‑6 (MQ=255)
               tCAAGGAAAGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTgg                                                                                                                                          <  3:264875/142‑1 (MQ=255)
               tCAAGGAAAGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTgg                                                                                                                                          >  4:264875/1‑142 (MQ=255)
                 aaGGAAAGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGcc                                                                                                                                  <  5:96688/148‑1 (MQ=255)
                 aaGGAAAGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGCGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCg                                                                                                                                 >  6:53031/1‑149 (MQ=255)
                   ggAAAGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGa                                                                                                                               >  8:128249/1‑149 (MQ=255)
                     aaaGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGAt                                                                                                                              >  3:65954/1‑148 (MQ=255)
                            cAAAACCCCAGATGTACTGATCACCCATCCTCTTCCCCCAAGCGGGAATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCt                                                                                                                      <  4:49914/149‑1 (MQ=255)
                             aaTACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACGGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCtt                                                                                                                     <  4:253153/149‑1 (MQ=255)
                             aaTACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCtt                                                                                                                     <  2:50036/149‑1 (MQ=255)
                                  cccAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCt                                                                                                                 >  4:121100/1‑148 (MQ=255)
                                  cccAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCt                                                                                                                 >  5:167632/1‑148 (MQ=255)
                                   ccAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCTGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCt                                                                                                                 >  1:30010/1‑147 (MQ=255)
                                   ccAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCt                                                                                                                 >  1:42946/1‑147 (MQ=255)
                                   ccAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCt                                                                                                                 >  1:39023/1‑147 (MQ=255)
                                   ccAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCt                                                                                                                 >  5:256503/1‑147 (MQ=255)
                                     aGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGc                                                                                                             >  2:249846/1‑149 (MQ=255)
                                           aCTGATCACCCATGCTCTTGCCTAAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGtt                                                                                                       >  7:168980/1‑149 (MQ=255)
                                            cTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGtt                                                                                                       <  8:168980/148‑1 (MQ=255)
                                                   cccATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCAc                                                                                               >  2:113710/1‑149 (MQ=255)
                                                     cATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCAccc                                                                                             >  3:214868/1‑149 (MQ=255)
                                                         ctctTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTg                                                                                          <  6:256503/148‑1 (MQ=255)
                                                         ctctTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTg                                                                                          <  3:121100/148‑1 (MQ=255)
                                                          tctTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTggg                                                                                        >  6:173395/1‑149 (MQ=255)
                                                               ccccAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAACCCACTCTTGGCAAGATTGGGCGGTGCCGGATAACGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGa                                                                                    <  3:152748/148‑1 (MQ=255)
                                                               ccccAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGa                                                                                    <  7:128249/148‑1 (MQ=255)
                                                               ccccAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGa                                                                                    <  4:214868/148‑1 (MQ=255)
                                                                 ccAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGt                                                                                 >  7:24053/1‑149 (MQ=255)
                                                                   aaGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTgg                                                                               <  1:249846/149‑1 (MQ=255)
                                                                                ggTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGcaca                                                                  <  2:30010/149‑1 (MQ=255)
                                                                                      cTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCAt                                                            >  1:75891/1‑149 (MQ=255)
                                                                                      cTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACACCTGGGGCGAGGTGGCAGCCAATGCACATAGCAt                                                            >  5:138750/1‑149 (MQ=255)
                                                                                         gATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTg                                                          <  1:113710/148‑1 (MQ=255)
                                                                                          aTACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTa                                                        <  3:69820/149‑1 (MQ=255)
                                                                                          aTACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTa                                                        <  2:39023/149‑1 (MQ=255)
                                                                                           tACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTa                                                        <  2:42946/148‑1 (MQ=255)
                                                                                                      aTGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTaccggcac                                             <  3:250978/148‑1 (MQ=255)
                                                                                                       tGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCATTCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACGCCTGGGGCGAGGTGGCAGCCAATGGACATAGCATTGTAGTTaccggcaccg                                           <  5:53031/149‑1 (MQ=255)
                                                                                                                  gCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAg                                >  4:178177/1‑149 (MQ=255)
                                                                                                                   cAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAgg                               <  3:150751/149‑1 (MQ=255)
                                                                                                                                aTCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATCCtt                    >  7:20098/1‑147 (MQ=255)
                                                                                                                                 tCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACGTAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATCCTTTAg                 <  8:24053/149‑1 (MQ=255)
                                                                                                                                 tCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATCCTTTAg                 >  3:103777/1‑149 (MQ=255)
                                                                                                                                   gCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATCCTTTAgg                <  8:20098/148‑1 (MQ=255)
                                                                                                                                    cTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATCCTTTAGGtt              <  6:138750/149‑1 (MQ=255)
                                                                                                                                     ttAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTacc                                                  <  1:103326/112‑1 (MQ=255)
                                                                                                                                     ttAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTacc                                                  >  2:103326/1‑112 (MQ=255)
                                                                                                                                       aaaCCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATCCTTTAGGTTc             <  5:237213/147‑1 (MQ=255)
                                                                                                                                        gcccAATCTTGCCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATCCTTTAGGTTCgg           >  6:164328/3‑148 (MQ=255)
                                                                                                                                        aaCCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATCCTTTAGGTTCgg           <  6:167632/148‑1 (MQ=255)
                                                                                                                                          ccAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACcggcaccggcac                                       >  3:9472/1‑118 (MQ=255)
                                                                                                                                          ccAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACcggcaccggcac                                       <  4:9472/118‑1 (MQ=255)
                                                                                                                                          ccAATCTTGCCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACACCTCGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCTGCACCGGCACGGATAAGGGTTTCATCCTTTa                  >  3:112393/1‑139 (MQ=255)
                                                                                                                                                  ggCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATCCTTTAGGTTCGGAACCTCGAt  <  3:178177/147‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
ATCTGCGGCATTTGCTCAAGGAAAGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATCCTTTAGGTTCGGAACCTCGAT  >  NC_002947/6129966‑6130258

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