Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,610,971 G→A 5.1% R89H (CGT→CAT)  hyfI → hydrogenase 4, Fe‑S subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,610,9710GA 5.1% 255.8 / 2.4 78R89H (CGT→CAT) hyfIhydrogenase 4, Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base G (32/42);  new base A (2/2);  total (34/44)
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.91e-01

ATTGCGGCGGCTGCAACGCCTGTGAAATCGAAATTTTTGCTGCCATTACACCAGTATTCGACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTAGCGACACCATTGTCCCCATTGATGTTT  >  NC_000913/2610824‑2611114
                                                                                                                                                   |                                                                                                                                               
aTTGCGGCGGCTGCAACGCCTGTGAAATCGAAATTTTTGCTGCCATTACACCAGTATTCGACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCg                                                                                                                                                 >  2:144257/1‑148 (MQ=255)
 ttGCGGCGGCTGCAACGCCTGTGAAATCGAAATTTTTGCTGCCATTACACCAGTATTCGACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTa                                                                                                                                               >  3:210122/1‑149 (MQ=255)
  tGCGGCGGCTGCAACGCCTGTGAAATCGAAATTTTTGCTGCCATTACACCAGTATTCGACGCAGAACGTTTTTGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTa                                                                                                                                               <  2:263687/148‑1 (MQ=255)
  tGCGGCGGCTGCAACGCCTGTGAAATCGAAATTTTTGCTGCCATTACACCAGTATTCGACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTa                                                                                                                                               <  4:324775/148‑1 (MQ=255)
         gcTGCAACGCCTGTGAAATCGAAATTTTTGCTGCCATTACACCAGTATTCGACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCa                                                                                                                                       >  3:302307/1‑149 (MQ=255)
          cTGCAACGCCTGTGAAATCGAAATTTTTGCTGCCATTACACCAGTATTCGACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCAc                                                                                                                                      <  4:210122/149‑1 (MQ=255)
             cAACGCCTGTGAAATCGAAATTTTTGCTGCCATTACACCAGTATTCGACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTc                                                                                                                                   >  3:237475/1‑149 (MQ=255)
             cAACGCCTGTGAAATCGAAATTTTTGCTGCCATTACACCAGTATTCGACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTc                                                                                                                                   >  3:17493/1‑149 (MQ=255)
               aCGCCTGTGAAATCGAAATTTTTGCTGCCATTACACCAGTATTCGACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCgg                                                                                                                                 <  4:122144/149‑1 (MQ=255)
                     gtgAATTCGAAATTTTTGCTGCCATTACTCCTGTATTCGACGGAGAACGTTTTGTCATTAAGGTTGTTTCTTCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCAACCGGGCGATGCGTATGCCTTCACTTCGGGCGTAt                                                                                                                           <  8:124025/149‑1 (MQ=255)
                      tgAAATCGAAATTTTTGCTGCCATTACACCAGTATTCGACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTa                                                                                                                            <  5:29098/147‑1 (MQ=255)
                      tgAAATCGAAATTTTTGCTGCCATTACACCAGTATTCGACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTAt                                                                                                                           <  8:106194/148‑1 (MQ=255)
                        aaaTCGAAATTTTTGCTGCCATTACACCAGTATTCGACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGa                                                                                                                         <  8:17544/148‑1 (MQ=255)
                          aTCGAAATTTTTGCTGCCATTACACCAGTATTCGACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTc                                                                                                                      >  3:297546/1‑149 (MQ=255)
                           tCGAAATTTTTGCTGCCATTACACCAGTATTCGACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCt                                                                                                                     <  2:328979/149‑1 (MQ=255)
                             gAAATTTTTGCTGCCATTACACCAGTATTCGACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGc                                                                                                                   >  6:88164/1‑149 (MQ=255)
                              aaaTTTTTGCTGCCATTACACCAGTATTCGACGCAGAACGTTTTAGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGCTGCGTATGCCTGCACTTCGGGCGTATGAGTCTGcc                                                                                                                  >  5:343004/1‑149 (MQ=255)
                                          ccATTACACCAGTATTCGACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATaa                                                                                                       >  7:244929/1‑148 (MQ=255)
                                           cATTACACCAGTATTCGACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATaaaa                                                                                                     <  6:342976/149‑1 (MQ=255)
                                               acacCAGTATTCGACGGTGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTAGTGGCGCTGTCACCCGTGCGATGGGTATGCCTTCACTTCTTGCGTATGAGTCTGCCCCCGATCATAAAAttt                                                                                                  <  1:53226/148‑1 (MQ=255)
                                                cacCAGTATTCGACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTtgt                                                                                                <  1:7390/149‑1 (MQ=255)
                                                 acCAGTATTCGACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTtgtg                                                                                               >  7:68768/1‑149 (MQ=255)
                                                 acCAGTATTCGACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACAGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTtgtg                                                                                               >  6:323381/1‑149 (MQ=255)
                                                   cAGTATTCGTCGCATAACGTTTTGGCCTTTAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGtt                                                                                             <  4:17493/149‑1 (MQ=255)
                                                   cAGTATTCGACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTtgt                                                                                                >  5:225862/1‑146 (MQ=255)
                                                   cAGTATTCGACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTtgt                                                                                                <  6:225862/146‑1 (MQ=255)
                                                   cAGTATTCGACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGtt                                                                                             >  7:249197/1‑149 (MQ=255)
                                                    aGTATTCGACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGttt                                                                                            <  6:236654/149‑1 (MQ=255)
                                                           gACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACg                                                                                      <  4:297546/148‑1 (MQ=255)
                                                                 gAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTg                                                                               <  1:365909/149‑1 (MQ=255)
                                                                   aCGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCg                                                                             <  2:331924/149‑1 (MQ=255)
                                                                   aCGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGAGGCGTATGCCTGCGCTTCGGGCGTATGAGGCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCg                                                                             <  1:248238/149‑1 (MQ=255)
                                                                     gTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGgt                                                                           <  4:302307/149‑1 (MQ=255)
                                                                      ttttGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGgtg                                                                          >  4:355170/1‑149 (MQ=255)
                                                                           gCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTcggc                                                                     >  4:121696/1‑149 (MQ=255)
                                                                           gCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTcgg                                                                      <  3:155350/148‑1 (MQ=255)
                                                                              ttAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGAGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCAGACAATTTGTGTTTCCTACGGCGCGTGCGGTGTcggcggc                                                                  >  8:314445/1‑149 (MQ=255)
                                                                                 aGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGGCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGGGCGGTGTcggcggcgg                                                                <  6:97356/148‑1 (MQ=255)
                                                                                  ggTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTa                                                              >  7:344742/1‑149 (MQ=255)
                                                                                  ggTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTa                                                              >  4:168782/1‑149 (MQ=255)
                                                                                  ggTTGTTTCATCACAGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGACCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGAGTCGGCGGCGGTa                                                              >  1:39445/1‑149 (MQ=255)
                                                                                       tttCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGgcgc                                                                                  <  4:191191/124‑1 (MQ=255)
                                                                                       tttCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGgcgc                                                                                  >  3:191191/1‑124 (MQ=255)
                                                                                         tcatcaCCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTcc                                                        <  8:249197/148‑1 (MQ=255)
                                                                                              aCCGCGTCACTCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCACAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGAt                                                   <  2:275158/148‑1 (MQ=255)
                                                                                              aCCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGAt                                                   <  4:237475/148‑1 (MQ=255)
                                                                                               ccGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGAtct                                                 <  3:355170/149‑1 (MQ=255)
                                                                                               ccGCGTCACGCCGATATTTTGGTATTTTCTGGACCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGCGCGTATGAGTCTGCCCCCGAGCATAAAATTTGTGTTTCCTACGTCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGAtct                                                 <  2:39445/149‑1 (MQ=255)
                                                                                                  cgTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTa                                              <  5:88164/149‑1 (MQ=255)
                                                                                                    tCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGTCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACa                                            <  1:144303/149‑1 (MQ=255)
                                                                                                    tCACGCCGATATTTTGTCAATTACTGGCGCAGTCACCCGGGCGCTGCGTATGCCTGCACTTCGGGCGTATGAGGCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGGGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACa                                            <  6:288212/149‑1 (MQ=255)
                                                                                                       cGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCg                                         <  3:343621/149‑1 (MQ=255)
                                                                                                           gATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTtgt                                                                                                <  8:243194/90‑1 (MQ=255)
                                                                                                           gATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTtgt                                                                                                >  7:243194/1‑90 (MQ=255)
                                                                                                           gATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTg                                     <  3:121696/149‑1 (MQ=255)
                                                                                                             tatTTTGTTATTTACTGGCGCAGTCAACCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTggg                                   <  7:342910/149‑1 (MQ=255)
                                                                                                                   gTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTGGGGCgg                               >  2:234337/1‑147 (MQ=255)
                                                                                                                    ttATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTAg                            >  8:366331/1‑149 (MQ=255)
                                                                                                                    ttATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTAg                            >  4:242937/1‑149 (MQ=255)
                                                                                                                     tttttACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTAGc                           >  3:153724/3‑149 (MQ=255)
                                                                                                                       tttACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGATTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGTCGGTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTAGCGa                         >  6:320000/1‑149 (MQ=255)
                                                                                                                          aCTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTAGCGacac                      >  1:246268/1‑149 (MQ=255)
                                                                                                                            tGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCt                                      <  5:272550/131‑1 (MQ=255)
                                                                                                                            tGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCt                                      >  6:272550/1‑131 (MQ=255)
                                                                                                                            tGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTAGCGACACCa                    >  1:170595/1‑149 (MQ=255)
                                                                                                                            tGGCGCAGTCACCCGGGCGATGCATATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTc                                       <  7:210640/130‑1 (MQ=255)
                                                                                                                            tGGCGCAGTCACCCGGGCGATGCATATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTc                                       >  8:210640/1‑130 (MQ=255)
                                                                                                                              gcgcAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTAGCGACACCAtt                  <  1:144257/149‑1 (MQ=255)
                                                                                                                                       cccGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTAGCGACACCAt                   >  5:195918/1‑139 (MQ=255)
                                                                                                                                       cccGGGCGATGAGTATGGCTGCACTTCGGGCGTATGAGTGTGCCCCCGATCATAAAATTTGTGTTTCCTAGGGCGCGTGCGGTGTCGGCTGCGGTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTAGCGACACCAt                   <  6:195918/139‑1 (MQ=255)
                                                                                                                                        ccGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTAGCGACACCATTGTCCCCAtt         >  4:39332/1‑148 (MQ=255)
                                                                                                                                          gggCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTAGCGACACCATTGTCCCCATTGAt      <  1:36546/149‑1 (MQ=255)
                                                                                                                                          gggCGATGCATATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTcgg                                                                      <  6:23660/85‑1 (MQ=255)
                                                                                                                                          gggCGATGCATATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTcgg                                                                      >  5:23660/1‑85 (MQ=255)
                                                                                                                                             cGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTGTACCGCGTCGGGGGCGGTAGCGACACCATTGTCCCCATTGATGtt   >  8:291120/1‑149 (MQ=255)
                                                                                                                                              gATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTAGCGACACCATTGTCCCCATTGATGttt  <  8:68768/149‑1 (MQ=255)
                                                                                                                                              gATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTAGCGACACCATTGTCCCCATTGATGttt  <  2:147570/149‑1 (MQ=255)
                                                                                                                                              gATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTAGCGACACCATTGTCCCCATTGATGttt  <  7:366331/149‑1 (MQ=255)
                                                                                                                                              gATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTAGCGACACCATTGTCCCCATTGATGttt  >  7:262892/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
ATTGCGGCGGCTGCAACGCCTGTGAAATCGAAATTTTTGCTGCCATTACACCAGTATTCGACGCAGAACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCGTATGCCTGCACTTCGGGCGTATGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTAGCGACACCATTGTCCCCATTGATGTTT  >  NC_000913/2610824‑2611114

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