Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 446,355 T→C N192S (AAC→AGC)  cyoE ← protoheme IX farnesyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913446,3550TC99.3% 482.8 / ‑4.7 150N192S (AAC→AGC) cyoEprotoheme IX farnesyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (69/80);  total (69/81)
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.32e-01

Only 100 of 155 total aligned reads displayed.

ATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGTTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATCGC  >  NC_000913/446265‑446443
                                                                                          |                                                                                        
aTCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGctt                                                                                        <  2:2503072/93‑1 (MQ=255)
 tCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATCCCGGCATGCTTg                                                                                       >  2:1510152/1‑93 (MQ=255)
 tCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTg                                                                                       >  1:3047421/1‑93 (MQ=255)
 tCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTg                                                                                       >  1:1865407/1‑93 (MQ=255)
  cAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTgc                                                                                      <  1:4275202/93‑1 (MQ=255)
  cAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTgc                                                                                      <  2:3256682/93‑1 (MQ=255)
   aGCGTGGCAACGGCAAAGGCGATTATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTgcc                                                                                     >  1:398769/1‑93 (MQ=255)
   aGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTgcc                                                                                     >  2:3609582/1‑93 (MQ=255)
    gCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTgccg                                                                                    >  1:4409446/1‑93 (MQ=255)
        ggCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTg                                                                                <  1:2140805/93‑1 (MQ=255)
        ggCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTg                                                                                <  1:4516499/93‑1 (MQ=255)
         gCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTgg                                                                               >  1:2705022/1‑93 (MQ=255)
          cAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGt                                                                              <  2:2019293/93‑1 (MQ=255)
          cAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGt                                                                              <  2:2320379/93‑1 (MQ=255)
           aaCGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTa                                                                             >  2:4119529/1‑93 (MQ=255)
            aCGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTaa                                                                            >  2:4263703/1‑93 (MQ=255)
             cGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAAt                                                                           >  2:2169976/1‑93 (MQ=255)
             cGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAAt                                                                           >  1:3013158/1‑93 (MQ=255)
             cGGCAAAGACGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAAt                                                                           <  2:2160688/93‑1 (MQ=255)
               gCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATcc                                                                         <  1:1534136/93‑1 (MQ=255)
                cAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCt                                                                        <  2:288744/93‑1 (MQ=255)
                  aaGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTa                                                                      >  2:4597592/1‑93 (MQ=255)
                    ggCGATCATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTaaa                                                                    <  1:1479329/93‑1 (MQ=255)
                       gatgatATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCg                                                                 >  1:1077918/1‑93 (MQ=255)
                       gatgatATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCg                                                                 >  1:3273373/1‑93 (MQ=255)
                       gatgatATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCg                                                                 >  2:2196410/1‑93 (MQ=255)
                       gatgatATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCg                                                                 <  1:1700530/93‑1 (MQ=255)
                         tgatATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGa                                                               <  1:1574110/93‑1 (MQ=255)
                           atatACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGaaa                                                             <  2:4595514/93‑1 (MQ=255)
                            tataCAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGaaaa                                                            <  2:2758080/93‑1 (MQ=255)
                            tataCAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGaaaa                                                            <  2:2874125/93‑1 (MQ=255)
                            tataCAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGaaaa                                                            >  2:2061/1‑93 (MQ=255)
                            tataCAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGaaaa                                                            <  2:3467123/93‑1 (MQ=255)
                             attcAGCGTGATGTTATTCTTCGCCACCGAAATTCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCCGGTAATCCTTAAAGCGGAAAAt                                                           <  1:561453/90‑1 (MQ=255)
                             ataCAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAAt                                                           <  1:3299826/93‑1 (MQ=255)
                             ataCAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAAt                                                           <  2:1329170/93‑1 (MQ=255)
                             ataCAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAAt                                                           <  2:4426271/93‑1 (MQ=255)
                             ataCAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAAt                                                           <  1:561454/93‑1 (MQ=255)
                             ataCAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAAt                                                           >  1:2021920/1‑93 (MQ=255)
                               aCAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATgg                                                         <  1:373346/93‑1 (MQ=255)
                                cAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGc                                                        <  1:3221564/93‑1 (MQ=255)
                                cAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGc                                                        <  1:1546329/93‑1 (MQ=255)
                                 aGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCg                                                       >  2:3103381/1‑93 (MQ=255)
                                  gCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAACGCGGAAACTGGCGa                                                      >  1:86471/1‑93 (MQ=255)
                                   cGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGAt                                                     <  1:4361064/93‑1 (MQ=255)
                                    gtgatgtgatTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATg                                                    <  1:3753832/93‑1 (MQ=255)
                                      gatgtgatTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGc                                                  >  1:2887554/1‑93 (MQ=255)
                                      gatgtgatTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGc                                                  >  2:2106437/1‑93 (MQ=255)
                                      gatgtgatTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATAGCGATGGc                                                  >  2:2940115/1‑93 (MQ=255)
                                         gtgatTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATa                                               <  1:2247219/93‑1 (MQ=255)
                                          tgatTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAg                                              >  2:3513331/1‑93 (MQ=255)
                                           gatTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAgg                                             <  2:3831239/93‑1 (MQ=255)
                                            atTCTTCGCTTCCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGa                                            >  2:4196684/1‑93 (MQ=255)
                                            atTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGa                                            <  1:239339/93‑1 (MQ=255)
                                                ttcGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGa                                        >  2:3700004/1‑93 (MQ=255)
                                                ttcGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGa                                        >  1:2842966/1‑93 (MQ=255)
                                                     cACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCAt                                   >  1:1495444/1‑93 (MQ=255)
                                                     cACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCAt                                   >  2:4240481/1‑93 (MQ=255)
                                                      cccGAAATGCCTTTTACCCCTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATc                                  <  1:62804/92‑1 (MQ=255)
                                                      aCCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATc                                  <  1:948264/93‑1 (MQ=255)
                                                      aCCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATc                                  <  1:4644073/93‑1 (MQ=255)
                                                       ccGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCt                                 >  2:336498/1‑93 (MQ=255)
                                                        cGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAg                                              >  1:4578901/1‑79 (MQ=255)
                                                        cGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAg                                              <  2:4578901/79‑1 (MQ=255)
                                                             tGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCcaca                           >  2:3666941/1‑93 (MQ=255)
                                                             tGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCcaca                           >  1:3506393/1‑93 (MQ=255)
                                                              gCCTTTTACTACTGGCAATACCGGAATGCTTGGCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCGTAGGAGTGAGGCATCTGCCACAg                          >  2:2603414/1‑93 (MQ=255)
                                                                cTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGc                        >  2:643372/1‑93 (MQ=255)
                                                                cTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGc                        <  1:1120217/93‑1 (MQ=255)
                                                                cTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGc                        <  1:930664/93‑1 (MQ=255)
                                                                 ttttACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCt                       <  1:1497438/93‑1 (MQ=255)
                                                                 ttttACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCt                       <  1:605759/93‑1 (MQ=255)
                                                                 ttttACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCt                       <  1:4081615/93‑1 (MQ=255)
                                                                   ttACCGCTGTCAATTCCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCAGCCACAGGCTGa                     <  2:3842045/93‑1 (MQ=255)
                                                                     aCCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAg                   >  2:4720881/1‑93 (MQ=255)
                                                                      ccACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGa                  <  2:2557949/93‑1 (MQ=255)
                                                                      ccACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGa                  <  1:1921711/93‑1 (MQ=255)
                                                                      ccACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGa                  <  2:4542508/93‑1 (MQ=255)
                                                                      ccACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGa                  >  2:4002017/1‑93 (MQ=255)
                                                                       cACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGAt                 >  2:2379562/1‑93 (MQ=255)
                                                                        aCTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATa                <  1:464670/93‑1 (MQ=255)
                                                                         cTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAg               >  1:3930954/1‑93 (MQ=255)
                                                                          tGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGc              <  2:4748810/93‑1 (MQ=255)
                                                                           ggCAATACCGGAATGCTTGCCGCCTGGCAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGcc             >  2:1510537/1‑93 (MQ=255)
                                                                             cAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCcag           <  1:1709020/93‑1 (MQ=255)
                                                                               aTACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCcagca         >  1:4087336/1‑93 (MQ=255)
                                                                               aTACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCcagca         >  1:4607357/1‑93 (MQ=255)
                                                                               aTACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCcagca         >  2:1343752/1‑93 (MQ=255)
                                                                                tACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCcagcag        >  2:639690/1‑93 (MQ=255)
                                                                                tACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCcagcag        <  2:84087/93‑1 (MQ=255)
                                                                                 aCCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAgg       <  1:2718948/93‑1 (MQ=255)
                                                                                 aCCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAgg       <  1:4548845/93‑1 (MQ=255)
                                                                                  ccGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGa      <  1:3169951/93‑1 (MQ=255)
                                                                                  ccGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGa      <  2:216456/93‑1 (MQ=255)
                                                                                   cGGAATGCTTGCCGCCTGGTAATCCTTAAATCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGAt     <  1:1391148/93‑1 (MQ=255)
                                                                                   cGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGAt     >  1:150123/1‑93 (MQ=255)
                                                                                   cGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGAt     <  2:2425762/93‑1 (MQ=255)
                                                                                    ggAATGCTTGCCTCCTGGTAATCCTTAAATCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATc    <  1:1385398/93‑1 (MQ=255)
                                                                                     gAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATCg   >  2:403570/1‑93 (MQ=255)
                                                                                      aaTGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATCGc  <  1:209541/93‑1 (MQ=255)
                                                                                          |                                                                                        
ATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGTTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATCGC  >  NC_000913/446265‑446443

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