Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,130,169 G→C 5.2% 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,1690GC 5.2% 234.4 / 4.5 77A107G (GCC→GGC) PP_5378metal‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (38/35);  new base C (2/2);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.54e-01

TGCGGCATTTGCTCAAGGAAAGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATCCTTTAGGTTCGGAACCTCGATGTCC  >  NC_002947/6130025‑6130318
                                                                                                                                                 |                                                                                                                                                     
tGCGGCATTTGCTCAAGGAAAGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGa                                                                                                                                                  >  1:229258/1‑150 (MQ=255)
         tGCTCAAGGAAAGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCGTACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGt                                                                                                                                         <  2:1182085/150‑1 (MQ=255)
               aGGAAAGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGTTGGTCGATTCAACCAATCTTGGCTAGATTGGGCGGTGCCTGa                                                                                                                                   >  2:615214/1‑150 (MQ=255)
                ggAAAGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGAt                                                                                                                                  >  2:113521/1‑150 (MQ=255)
                 gAAAGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGAAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCGGGATa                                                                                                                                 >  1:658789/1‑150 (MQ=255)
                  aaaGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATaa                                                                                                                                <  2:673319/150‑1 (MQ=255)
                     gCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAgc                                                                                                                             >  2:156433/1‑150 (MQ=255)
                               cccAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCtt                                                                                                                         <  1:166681/144‑1 (MQ=255)
                               cccAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCtt                                                                                                                         >  2:166681/1‑144 (MQ=255)
                               cccAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGa                                                                                                                   <  1:125889/150‑1 (MQ=255)
                               cccAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGa                                                                                                                   <  1:792351/150‑1 (MQ=255)
                               cccAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGa                                                                                                                   >  2:125889/1‑150 (MQ=255)
                               cccAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGa                                                                                                                   >  2:206954/1‑150 (MQ=255)
                                ccAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGc                                                                                                                                            <  1:788819/124‑1 (MQ=255)
                                ccAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGc                                                                                                                                            >  2:788819/1‑124 (MQ=255)
                                    aTGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCg                                                                                                              >  1:311549/1‑150 (MQ=255)
                                      gTACTGATCACCTATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCAAGAACGCAGCTGCGTGGTTCAT‑GATTAAACCAACCTTGGCAAGATGTGGCGCTGCCGCATCAAGCGCGTGGGTGTGCTCggg                                                                                                            >  2:228192/1‑149 (MQ=255)
                                        aCTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGttt                                                                                                          >  2:935802/1‑150 (MQ=255)
                                            aTCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCg                                                                                                              >  2:1082617/1‑142 (MQ=255)
                                            aTCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCg                                                                                                              <  1:1082617/142‑1 (MQ=255)
                                              cACCCCTGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCa                                                                                                    <  1:228192/150‑1 (MQ=255)
                                               ccccATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCAc                                                                                                   <  2:660759/149‑1 (MQ=255)
                                                 ccATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCAccc                                                                                                 >  1:109733/1‑150 (MQ=255)
                                                 ccATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCAccc                                                                                                 >  1:774783/1‑150 (MQ=255)
                                                      ctctTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCTGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTggg                                                                                            <  1:615214/150‑1 (MQ=255)
                                                      ctctTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTggg                                                                                            <  1:972416/150‑1 (MQ=255)
                                                      ctctTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTggg                                                                                            <  1:364973/150‑1 (MQ=255)
                                                       tctTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAgc                                                                                                                             <  1:1108062/116‑1 (MQ=255)
                                                       tctTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAgc                                                                                                                             >  2:1108062/1‑116 (MQ=255)
                                                        ctTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGc                                                                                          <  2:229258/150‑1 (MQ=255)
                                                              ccAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGtt                                                                                                           >  2:774789/1‑127 (MQ=255)
                                                              ccAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGtt                                                                                                           <  1:774789/127‑1 (MQ=255)
                                                              ccAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTg                                                                                    >  1:485714/1‑150 (MQ=255)
                                                                 aGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCa                                                                                 <  1:701135/150‑1 (MQ=255)
                                                                 aGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCa                                                                                 >  2:701135/1‑150 (MQ=255)
                                                                      ccATACCGCTACCTTTGATACCGTGTTTAATGTTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTCCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCaa                                                                            <  2:63729/149‑1 (MQ=255)
                                                                           ccGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGcac                                                                       <  1:206954/150‑1 (MQ=255)
                                                                               tAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCaa                                                                            >  2:241866/1‑141 (MQ=255)
                                                                               tAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCaa                                                                            <  1:241866/141‑1 (MQ=255)
                                                                                 gCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAg                                                                   <  2:115939/148‑1 (MQ=255)
                                                                                 gCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAg                                                                   >  1:115939/1‑148 (MQ=255)
                                                                                      gATACCGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTa                                                            <  1:156433/150‑1 (MQ=255)
                                                                                         aCCGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGc                                                                               <  1:692793/128‑1 (MQ=255)
                                                                                         aCCGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGc                                                                               >  2:692793/1‑128 (MQ=255)
                                                                                          ccGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTa                                                        >  1:732065/1‑150 (MQ=255)
                                                                                           cGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTac                                                       >  1:395204/1‑150 (MQ=255)
                                                                                            gtgtTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCaa                                                                            <  1:492591/128‑1 (MQ=255)
                                                                                            gtgtTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCaa                                                                            >  2:492591/1‑128 (MQ=255)
                                                                                              gtTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTaccgg                                                    <  1:994343/150‑1 (MQ=255)
                                                                                                   aTGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTaccggcaccg                                               >  2:402970/1‑150 (MQ=255)
                                                                                                      gTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATc                                                                                                      >  1:680621/1‑92 (MQ=255)
                                                                                                      gTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATc                                                                                                      <  2:680621/92‑1 (MQ=255)
                                                                                                       tCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCgg                                                                                                             >  2:183498/1‑84 (MQ=255)
                                                                                                       tCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCgg                                                                                                             <  1:183498/84‑1 (MQ=255)
                                                                                                          aGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGa                                        <  1:935802/150‑1 (MQ=255)
                                                                                                            aaCGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGt                                                             <  2:774768/127‑1 (MQ=255)
                                                                                                            aaCGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGt                                                             >  1:774768/1‑127 (MQ=255)
                                                                                                                 aGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAgggg                                 >  2:303337/1‑149 (MQ=255)
                                                                                                                    tGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTc                              <  2:658789/150‑1 (MQ=255)
                                                                                                                    tGCG‑GGGACATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTc                              >  2:175777/1‑150 (MQ=255)
                                                                                                                       g‑gggTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATc                           <  2:485714/150‑1 (MQ=255)
                                                                                                                              aTCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGa                                                                                        <  1:506314/83‑1 (MQ=255)
                                                                                                                              aTCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGa                                                                                        >  2:506314/1‑83 (MQ=255)
                                                                                                                                cGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGtt                                <  2:630538/137‑1 (MQ=255)
                                                                                                                                cGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGtt                                >  1:630538/1‑137 (MQ=255)
                                                                                                                                      gcccAATCTTGCCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATCCtt                        <  1:855698/137‑1 (MQ=255)
                                                                                                                                      gcccAATCTTGCCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATCCtt                        >  2:855698/3‑139 (MQ=255)
                                                                                                                                      gcccAATCTTGCCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATCCTTTa                      >  1:528693/3‑141 (MQ=255)
                                                                                                                                      gcccAATCTTGCCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATCCTTTa                      <  2:528693/139‑1 (MQ=255)
                                                                                                                                       aCCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATCCTTTAGGTTCGGAAc            <  2:519046/150‑1 (MQ=255)
                                                                                                                                       aCCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATCCTTTAGGTTCGGAAc            >  2:1106905/1‑150 (MQ=255)
                                                                                                                                       aCCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATACTTTAGGtttgagcc            >  2:68126/1‑144 (MQ=255)
                                                                                                                                        ccAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATCCTTTAGGTTCGGAAcc           >  2:530481/1‑150 (MQ=255)
                                                                                                                                          aaTCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAAGGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGTGTTTCATCCTTTAGGTTCGGAACCTc         <  2:311549/150‑1 (MQ=255)
                                                                                                                                           aTCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATCCTTTAGGTTCGGAACCTCg        >  1:457584/1‑150 (MQ=255)
                                                                                                                                            tCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATCCTTTAGGTTCGGAACCTCGa       >  2:869390/1‑150 (MQ=255)
                                                                                                                                             tttGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATCCTTTAGGTTCGGAACCTCGAt      <  2:732065/149‑1 (MQ=255)
                                                                                                                                                ggCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATCCTTTAGGTTCGGAACCTCGATGTc   >  2:407922/1‑150 (MQ=255)
                                                                                                                                                 gCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATCCTTTAGGTTCGGAACCTCGATGTcc  <  1:402970/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
TGCGGCATTTGCTCAAGGAAAGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCG‑GGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATCCTTTAGGTTCGGAACCTCGATGTCC  >  NC_002947/6130025‑6130318

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑