Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,793,898 Δ1 bp coding (212/552 nt) btuE ← glutathione peroxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,793,8970T.98.7% 337.2 / ‑1.8 76coding (213/552 nt)btuEglutathione peroxidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base . (36/38);  minor base A (1/0);  total (38/38)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCG  >  NC_000913/1793751‑1794041
                                                                                                                                                   |                                                                                                                                                
gcgACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCtgt                                                                                                                                                >  1:1073758/1‑147 (MQ=255)
     tGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAg                                                                                                                                          >  1:295599/1‑149 (MQ=255)
         gTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAAc                                                                                                                                      <  2:964586/149‑1 (MQ=255)
         gTCGGCGCTGCGGCAATCAATTTTTGATCAGCCGCGTTGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAAc                                                                                                                                      >  2:303292/1‑149 (MQ=255)
          tCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACt                                                                                                                                     >  1:591799/1‑149 (MQ=255)
                  gCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTCGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCAc                                                                                                                             <  2:1031859/149‑1 (MQ=255)
                    ggCAATCCAGGTTTGATACAGCGGAGGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACgg                                                                                                                           <  1:1094113/149‑1 (MQ=255)
                     gCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCAAGAAACGGGTTGCACggg                                                                                                                          >  1:528369/1‑149 (MQ=255)
                      caatccatTTTTGATACAGCGGCTGGTGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCCGACAATGGTGGCGCGGGa                                                                                                                         <  2:711207/149‑1 (MQ=255)
                        atcaatTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAAt                                                                                                                       >  1:1198085/1‑149 (MQ=255)
                         tcaatTTTTGATACAGCGAATAGCGTGCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATc                                                                                                                      <  1:253391/149‑1 (MQ=255)
                           aatTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATccc                                                                                                                    >  1:680486/1‑149 (MQ=255)
                           aatTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATccc                                                                                                                    >  2:550823/1‑149 (MQ=255)
                             tttttGATACAGCGGAGGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAg                                                                                                                  <  1:840163/149‑1 (MQ=255)
                                 tGATACTGCGGACTGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTTATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCAcc                                                                                                              <  2:594873/149‑1 (MQ=255)
                                  gATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCa                                                                                                             >  2:759008/1‑149 (MQ=255)
                                    tACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATa                                                                                                           >  1:549232/1‑149 (MQ=255)
                                     ataGCCGATCGCGTCCTGCGCCATTAACTTCAATCTTACAGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCCGTCCCAGAAACTGGTTGCACGGGAATCCCAGCACCATaa                                                                                                          <  1:1260868/147‑1 (MQ=255)
                                       aGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATaaaa                                                                                                        <  1:792853/149‑1 (MQ=255)
                                       aGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATaaaa                                                                                                        <  2:726746/149‑1 (MQ=255)
                                        gCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTCCTGTTCCAGAAACTGGTTGCACGGGAATCCCCGCACCATacacc                                                                                                       >  1:14253/1‑145 (MQ=255)
                                          cGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTATGTTCCGGAAACTGGTTGCACGGGAATCCCAGCACCAGAAAAcc                                                                                                      >  2:1144382/1‑149 (MQ=255)
                                            gaaggcggcCTTCGCCCTTAACTTCAATCTTACGGACCATCGGGAACGTCCCCCCCCATGTGGTGGTACAGGAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCtc                                                                                                    <  2:1130686/146‑1 (MQ=255)
                                            gATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCTATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAAcccc                                                                                                    >  2:905526/1‑147 (MQ=255)
                                            gATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCtc                                                                                                    <  1:1389589/149‑1 (MQ=255)
                                             aTGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCGCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCtcg                                                                                                   >  2:431046/1‑149 (MQ=255)
                                               ggCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCGGCACCATAAAACCtcgat                                                                                                 >  1:151830/1‑149 (MQ=255)
                                               ggCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCtcgat                                                                                                 <  1:761503/149‑1 (MQ=255)
                                                  ggccatCGCCATTAACGTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCtcgatcga                                                                                              <  2:753272/144‑1 (MQ=255)
                                                  gTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCtcgatcga                                                                                              >  2:1110383/1‑149 (MQ=255)
                                                      ttCGCCATTCACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCa                                                                                          <  1:592225/149‑1 (MQ=255)
                                                       tCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAg                                                                                         >  1:473110/1‑149 (MQ=255)
                                                        cGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCACGAAAGTGATTGCACGGGAATCCCAGCACCATAAAACCtcggtcgtacctgg                                                                                        >  1:38169/1‑142 (MQ=255)
                                                              tAACTCCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCt                                                                                  <  2:289628/149‑1 (MQ=255)
                                                                 cTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGaa                                                                               <  2:1291425/149‑1 (MQ=255)
                                                                  ttCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCtacat                                                                              >  2:349254/1‑145 (MQ=255)
                                                                     aaTCTTATTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAATCAGTTTTCAGCTCTTCACCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAACACCTCGATCGACCCAGGCTTTCTGAATAtt                                                                           <  2:641274/149‑1 (MQ=255)
                                                                      aTCTTATTGAACATCCGTAACGTCACCCCCCATGTGGTAGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGTACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATtc                                                                          <  2:373116/149‑1 (MQ=255)
                                                                           aCTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATGTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCaa                                                                     >  1:549131/1‑149 (MQ=255)
                                                                             tGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACt                                                                   >  1:829354/1‑149 (MQ=255)
                                                                                aCATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCt                                                                >  1:290413/1‑149 (MQ=255)
                                                                                  aTCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTTCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCa                                                              >  1:961712/1‑149 (MQ=255)
                                                                                  aTCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCAGCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCCCCAACTGCTCa                                                              >  1:962127/1‑149 (MQ=255)
                                                                                         aCGTCACCCCCCATGTGGTGGTACAGTAATTTTTAATCTCTTCATCGCTGCCCTGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTgcg                                                       <  2:1030314/149‑1 (MQ=255)
                                                                                         aCGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAATTGCTCATATTgcg                                                       <  2:967732/149‑1 (MQ=255)
                                                                                           gTCACCGCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCATGCGTTCTGAATATTCTCCAACTGCTCATATTgcggc                                                     <  2:1276022/149‑1 (MQ=255)
                                                                                            tCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTgcggcg                                                    <  1:1155767/149‑1 (MQ=255)
                                                                                              tccccccAGGTGGTGGTACAGGAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGtt                                                  <  2:279883/148‑1 (MQ=255)
                                                                                               ccccccATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTa                                                 <  2:220758/149‑1 (MQ=255)
                                                                                                cccccATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTaa                                                <  2:945783/149‑1 (MQ=255)
                                                                                                cccccATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTaa                                                >  1:244859/1‑149 (MQ=255)
                                                                                                cccccATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCGGAATATTCTCCAACTGCTCATATTGCGGCGTTaa                                                >  1:1296/1‑149 (MQ=255)
                                                                                                 ccccATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATACTGCGGCGTTAAg                                               >  1:912929/1‑149 (MQ=255)
                                                                                                  cccATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAGCTGGTTGCACGGGAATCCCAGCACCATAACACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGc                                              <  2:295599/149‑1 (MQ=255)
                                                                                                            ggtACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCAAGGGAATACAAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGa                                    >  2:1264575/1‑149 (MQ=255)
                                                                                                              tACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAgg                                  <  1:303292/149‑1 (MQ=255)
                                                                                                               aCAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAATCTGGTGGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGc                                 <  2:1332277/149‑1 (MQ=255)
                                                                                                                  gTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGAc                              >  2:426785/1‑149 (MQ=255)
                                                                                                                   tAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGACTCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGAcc                             >  2:1209077/1‑148 (MQ=255)
                                                                                                                       ttttAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGTTCGACCCAGACTTTCGGAATATTCTCCAACTGCTCCTATTGCGCCGTTAAGCCACACTTTGAGGCGACATTGa                         <  1:905526/149‑1 (MQ=255)
                                                                                                                       ttttAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGGCCCAGGCTTTCTTCTTTAACTCCAACTGCTCATATTGCAGCGATAAGCCACACTTTGAGGCGACATTGa                         >  1:611929/1‑149 (MQ=255)
                                                                                                                         ttAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACa                       >  1:143632/1‑149 (MQ=255)
                                                                                                                          tAATCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACaa                      <  2:836523/149‑1 (MQ=255)
                                                                                                                            aTCTCTTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATc                    >  2:797223/1‑149 (MQ=255)
                                                                                                                                ctTCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGGCATTGACAATCAAcc                >  1:1116317/1‑148 (MQ=255)
                                                                                                                                 ttCATCGCTGCCCGGTTCTGTTCCAGAACCTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAg               <  1:757414/149‑1 (MQ=255)
                                                                                                                                 ttCATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAg               >  2:359817/1‑149 (MQ=255)
                                                                                                                                 ttCATCGCTGCCCGGTGCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAg               <  1:822749/149‑1 (MQ=255)
                                                                                                                                 ttCAGCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAACCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAg               <  2:528369/149‑1 (MQ=255)
                                                                                                                                   cagcGCTGCCCGGGTCTGTACCATAATCTGGTTGCACGGGATGCCCAGCACCATCAACCCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACGGCTCATATGGCGGCGTTAAGCCACGCGTTGAGGCGACATTGACAATCGACGGca             <  2:1327304/146‑1 (MQ=255)
                                                                                                                                   cATCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGca             <  1:550823/149‑1 (MQ=255)
                                                                                                                                    aTCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGcac            >  1:243724/1‑149 (MQ=255)
                                                                                                                                     tCGCTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGcaca           >  1:345526/1‑149 (MQ=255)
                                                                                                                                        cTGCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACAttc        >  2:509060/1‑148 (MQ=255)
                                                                                                                                          gCCCGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTAcc      <  2:162303/149‑1 (MQ=255)
                                                                                                                                            ccGGTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCgg    <  1:759008/149‑1 (MQ=255)
                                                                                                                                              ggTTCTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCCTACCGCGCGGGTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCCCATTACCGGCg  >  1:1304294/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCG  >  NC_000913/1793751‑1794041

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