Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,849,789 C→T 5.5% A21T (GCC→ACC)  hycB ← hydrogenase 3, Fe‑S subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,849,7890CT 5.5% 228.8 / 2.7 72A21T (GCC→ACC) hycBhydrogenase 3, Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base C (31/37);  new base T (2/2);  total (33/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-01

GGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATC  >  NC_000913/2849642‑2849931
                                                                                                                                                   |                                                                                                                                              
gggTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAggcg                                                                                                                                               >  8:148371/1‑149 (MQ=255)
  gTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAggcggc                                                                                                                                             <  1:370569/149‑1 (MQ=255)
      tGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTca                                                                                                                                         <  7:295140/149‑1 (MQ=255)
          gTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGTCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGTGGCCTCACAAg                                                                                                                                     >  3:44901/1‑149 (MQ=255)
          gTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAg                                                                                                                                     >  7:419299/1‑149 (MQ=255)
             aaCCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAgtg                                                                                                                                   <  6:337794/148‑1 (MQ=255)
              aCCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAgtgt                                                                                                                                  <  2:224989/148‑1 (MQ=255)
               ccGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCATTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTgg                                                                                                                                >  1:348639/1‑149 (MQ=255)
                        ccaccGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATa                                                                                                                       <  1:44834/149‑1 (MQ=255)
                        ccaccGCGCAGGGGGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCCGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATa                                                                                                                       <  5:229707/149‑1 (MQ=255)
                         caccGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATAc                                                                                                                      <  3:246492/149‑1 (MQ=255)
                         cacaGCGCAGGGTGAATCTTCACAGTGGAGACAGAGCTGCGGCGCAGATTCTTTTTCATTCATCATCACTCTCAGGCCCAGAATTGATTGCAGCCCGTGCAGGCGATGAGTATCTGAACAGGCGGCCTCACAAGTTTGGCAGCCGATa                                                                                                                       >  8:105123/1‑148 (MQ=255)
                           ccGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATa                                                                                                                       >  7:95502/1‑146 (MQ=255)
                                 aGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCGGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTg                                                                                                              <  3:119433/149‑1 (MQ=255)
                                    gTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACa                                                                                                                     <  2:171910/139‑1 (MQ=255)
                                    gTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACa                                                                                                                     >  1:171910/1‑139 (MQ=255)
                                     tGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGGGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGt                                                                                                          <  1:84609/149‑1 (MQ=255)
                                       cATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCa                                                                                                        <  8:122039/149‑1 (MQ=255)
                                        aTCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAg                                                                                                       <  2:296452/149‑1 (MQ=255)
                                        aTCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATCCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAg                                                                                                       >  5:191110/1‑149 (MQ=255)
                                           ttCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGAt                                                                                                                        <  4:383472/129‑1 (MQ=255)
                                           ttCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGAt                                                                                                                        >  3:383472/1‑129 (MQ=255)
                                            tCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGc                                                                                                                               <  8:3168/121‑1 (MQ=255)
                                            tCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGc                                                                                                                               >  7:3168/1‑121 (MQ=255)
                                              ccaGTGGGGACAGAGCTGCGGCGCAGTTTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGTCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTa                                                                                                 <  6:268454/148‑1 (MQ=255)
                                              acaGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAAtt                                                                                                  <  6:267569/148‑1 (MQ=255)
                                              acaGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTa                                                                                                 >  5:299734/1‑149 (MQ=255)
                                                aGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTAc                                                                                                >  2:216788/1‑148 (MQ=255)
                                                           agCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCAc                                                                                    <  3:283911/149‑1 (MQ=255)
                                                             cTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCa                                                                                  <  8:95502/149‑1 (MQ=255)
                                                                cggcgCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAgg                                                                                <  8:419299/148‑1 (MQ=255)
                                                                    gcAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGACTCACCAGGCATTc                                                                           >  7:32463/1‑149 (MQ=255)
                                                                      aGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCt                                                                         <  8:50203/149‑1 (MQ=255)
                                                                       gATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTc                                                                        <  6:101709/149‑1 (MQ=255)
                                                                         ttCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAg                                                                      <  5:8957/149‑1 (MQ=255)
                                                                         ttCTTTTTCATCCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAg                                                                      <  6:191110/149‑1 (MQ=255)
                                                                          tCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTCGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAgg                                                                     >  5:284278/1‑149 (MQ=255)
                                                                            tttttCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGTGGCCTCACAAGTGTGGCAGCCGATACaga                                                                                                                   >  6:77701/1‑101 (MQ=255)
                                                                            tttttCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGTGGCCTCACAAGTGTGGCAGCCGATACaga                                                                                                                   <  5:77701/101‑1 (MQ=255)
                                                                              tttCATTCAGCATCACTCTCTGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGTGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGa                                                                  <  4:44901/148‑1 (MQ=255)
                                                                              tttCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGtgtg                                                                                                                                 >  4:307618/1‑85 (MQ=255)
                                                                               ttcattcaGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACaa                                                                                              <  1:321801/119‑1 (MQ=255)
                                                                               ttcattcaGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACaa                                                                                              >  2:321801/1‑119 (MQ=255)
                                                                               ttcattcaGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCGGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGAt                                                                 <  4:43112/148‑1 (MQ=255)
                                                                                tcattcaGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGAtt                                                                <  6:299734/148‑1 (MQ=255)
                                                                                tcattcaGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGAtt                                                                <  7:105123/148‑1 (MQ=255)
                                                                                  attcaGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGt                                                              >  6:246605/1‑148 (MQ=255)
                                                                                  attcaGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGAGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGt                                                              >  8:220316/1‑148 (MQ=255)
                                                                                      aGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTCTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATGCACCAGGCATTCCTCAGGTGATTGTCAgtt                                                         >  5:241680/1‑149 (MQ=255)
                                                                                          tCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGa                                                     <  1:293275/149‑1 (MQ=255)
                                                                                          tCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGa                                                     >  6:300191/1‑149 (MQ=255)
                                                                                            aCTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGAc                                                    <  5:413061/148‑1 (MQ=255)
                                                                                              tctcAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAgtgt                                                                                                                                  <  1:166271/68‑1 (MQ=255)
                                                                                              tctcAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAgtgt                                                                                                                                  >  2:166271/1‑68 (MQ=255)
                                                                                                tcAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGaaaa                                               >  4:390794/1‑149 (MQ=255)
                                                                                                tcAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGaaaa                                               >  7:5150/1‑149 (MQ=255)
                                                                                                 cAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGaaaa                                               >  7:317764/1‑148 (MQ=255)
                                                                                                   ggCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACAt                                            <  1:216788/149‑1 (MQ=255)
                                                                                                     cgcgGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGc                                          <  2:348639/149‑1 (MQ=255)
                                                                                                           attgattgCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATg                                     >  6:230487/1‑148 (MQ=255)
                                                                                                            ttgattgCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAg                                                                                 <  2:209682/103‑1 (MQ=255)
                                                                                                            ttgattgCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAg                                                                                 >  1:209682/1‑103 (MQ=255)
                                                                                                              gattgCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAAt                                 >  8:377903/1‑149 (MQ=255)
                                                                                                              gattgCAGGCCGTGCTGGCGATGCGTCTATGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGGGGAGTCAGCAATTACAAAACTATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGaaa                                  <  5:363001/148‑1 (MQ=255)
                                                                                                               attgCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATg                                >  5:271971/1‑149 (MQ=255)
                                                                                                                   cAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCa                             <  5:230487/148‑1 (MQ=255)
                                                                                                                   cAGGCCGTGCTGCCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACTGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGCAAACATGCCGATGAAATGTCAt                            >  1:123856/1‑149 (MQ=255)
                                                                                                                       ccGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGACTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTc                        <  8:32463/149‑1 (MQ=255)
                                                                                                                         gTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGa                      <  8:317764/149‑1 (MQ=255)
                                                                                                                                 gATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATc              >  4:24767/1‑149 (MQ=255)
                                                                                                                                   tGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGa            <  5:300191/149‑1 (MQ=255)
                                                                                                                                       tctcTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCCTCGACACg        >  7:255859/1‑149 (MQ=255)
                                                                                                                                           tGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGccc     >  6:336587/1‑148 (MQ=255)
                                                                                                                                             aaCAGGCGGCCTCACAAGTGTGGCAGCAGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATCTCATATTCGACACTCATCGAAACGCCCATc  >  1:389650/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
GGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATC  >  NC_000913/2849642‑2849931

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