Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 444,641 C→T intergenic (‑126/‑42) panE ← / → yajQ 2‑dehydropantoate reductase, NADPH‑specific/phage Phi6 host factor, ATP/GTP binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913444,6410CT98.8% 256.9 / ‑5.7 81intergenic (‑126/‑42)panE/yajQ2‑dehydropantoate reductase, NADPH‑specific/phage Phi6 host factor, ATP/GTP binding protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (44/36);  total (44/37)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTCATTGTTCCTCCTCACCCGCAACCACTCCGGGTGTTCAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGCGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAACGCGAGCCGCGAAGTGGAGTCCCGTTTTGACT  >  NC_000913/444509‑444776
                                                                                                                                    |                                                                                                                                       
ttttCATTGTTCCTCCTCACCCGCAACCACTCCGGGTGTTCAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGtat                                                                                                                                   >  2:1358968/1‑139 (MQ=255)
         ttCCTCCTCACCCGCAACCACTCCGGGTGTTCAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCg                                                                                                                          >  1:707949/1‑139 (MQ=255)
                      gCAACCACTCCGGGTGTTCAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGaa                                                                                                             >  1:1073179/1‑139 (MQ=255)
                      gCAACCACTCCGGGTGTTCAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGaa                                                                                                             >  2:265888/1‑139 (MQ=255)
                      gCAACCACTCCGGGTGTTCAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGaa                                                                                                             >  1:153962/1‑139 (MQ=255)
                                 gggTGTTCAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGaaa                                                                                                  <  1:265888/139‑1 (MQ=255)
                                      ttCAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATg                                                                                             >  2:1007705/1‑139 (MQ=255)
                                       tCAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGc                                                                                            <  1:327452/139‑1 (MQ=255)
                                        cAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGcc                                                                                           <  2:153962/139‑1 (MQ=255)
                                        cAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGcc                                                                                           <  2:653891/139‑1 (MQ=255)
                                        cAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGc                                                                                            <  2:1006555/138‑1 (MQ=255)
                                        cAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGc                                                                                            >  1:1006555/1‑138 (MQ=255)
                                               gCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTc                                                                                    <  1:1007705/139‑1 (MQ=255)
                                                cTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGa                                                                                                              >  2:656707/1‑112 (MQ=255)
                                                cTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGa                                                                                                              <  1:656707/112‑1 (MQ=255)
                                                             gtgCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAgg                                                                                                           >  1:149034/1‑102 (MQ=255)
                                                             gtgCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAgg                                                                                                           <  2:61567/102‑1 (MQ=255)
                                                             gtgCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAgg                                                                                                           >  1:61567/1‑102 (MQ=255)
                                                             gtgCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAgg                                                                                                           <  2:149034/102‑1 (MQ=255)
                                                             gtgCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGa                                                                      >  2:581520/1‑139 (MQ=255)
                                                             gtgCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGa                                                                      >  2:505416/1‑139 (MQ=255)
                                                              tgCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCaaa                                                                                                                   <  1:45915/93‑1 (MQ=255)
                                                              tgCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCaaa                                                                                                                   >  2:45915/1‑93 (MQ=255)
                                                              tgCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCa                                                                                                                     >  1:759981/1‑91 (MQ=255)
                                                              tgCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCa                                                                                                                     <  2:759981/91‑1 (MQ=255)
                                                              tgCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGc                                                                                            >  1:110289/1‑116 (MQ=255)
                                                              tgCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGc                                                                                            <  2:110289/116‑1 (MQ=255)
                                                              tgCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCAGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAg                                                                    >  2:832514/1‑139 (MQ=255)
                                                                cATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGATTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGt                                                                   >  2:60108/1‑139 (MQ=255)
                                                                        tGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTc                                                           >  2:1103975/1‑139 (MQ=255)
                                                                           gACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGa                                                       <  1:832514/139‑1 (MQ=255)
                                                                               aTGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCa                                                                                                                     >  1:679471/1‑74 (MQ=255)
                                                                               aTGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCa                                                                                                                     <  2:679471/74‑1 (MQ=255)
                                                                               aTGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATc                                                                                        >  1:578594/1‑103 (MQ=255)
                                                                               aTGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATc                                                                                        <  2:578594/103‑1 (MQ=255)
                                                                               aTGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATAtt                                                                              >  2:902758/1‑113 (MQ=255)
                                                                               aTGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGCGGGTATTATGCTTCGCATCTAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATAtt                                                                              <  1:902758/113‑1 (MQ=255)
                                                                                 gCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCAc                                                  >  2:105255/1‑139 (MQ=255)
                                                                                 gCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCAc                                                  >  2:939093/1‑139 (MQ=255)
                                                                                        tAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAAcgc                                           <  2:1073179/139‑1 (MQ=255)
                                                                                          aGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGa                                                                                               >  1:887171/1‑85 (MQ=255)
                                                                                          aGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGa                                                                                               <  2:887171/85‑1 (MQ=255)
                                                                                          aGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCgg                                         <  1:939093/139‑1 (MQ=255)
                                                                                            ggATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTc                                       >  2:194054/1‑139 (MQ=255)
                                                                                                 cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCa                                                                                                                     >  1:1026663/1‑56 (MQ=255)
                                                                                                 cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCa                                                                                                                     <  2:1026663/56‑1 (MQ=255)
                                                                                                 cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAgg                                                                                                           >  2:685963/1‑66 (MQ=255)
                                                                                                 cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAgg                                                                                                           <  1:685963/66‑1 (MQ=255)
                                                                                                 cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGa                                                                                                    >  1:334879/1‑73 (MQ=255)
                                                                                                 cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGa                                                                                                    <  2:334879/73‑1 (MQ=255)
                                                                                                 cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCa                                                                                          >  2:567809/1‑83 (MQ=255)
                                                                                                 cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCa                                                                                          <  1:567809/83‑1 (MQ=255)
                                                                                                 cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGa                                                                <  1:1142597/109‑1 (MQ=255)
                                                                                                 cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGa                                                                >  2:1142597/1‑109 (MQ=255)
                                                                                                 cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGc                                                    <  2:9268/121‑1 (MQ=255)
                                                                                                 cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGc                                                    >  1:9268/1‑121 (MQ=255)
                                                                                                 cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATaa                                  >  2:1355588/1‑139 (MQ=255)
                                                                                                 cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATaa                                  >  1:1138244/1‑139 (MQ=255)
                                                                                                 cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATaa                                  >  1:1353140/1‑139 (MQ=255)
                                                                                                 cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATAAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATaa                                  >  1:882267/1‑139 (MQ=255)
                                                                                                       ggTGCAGGATTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAACGCGAg                            <  1:60108/139‑1 (MQ=255)
                                                                                                       ggTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATc                                                                                        <  2:716616/79‑1 (MQ=255)
                                                                                                       ggTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATc                                                                                        >  1:716616/1‑79 (MQ=255)
                                                                                                       ggTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAACGCGAg                            <  2:707949/139‑1 (MQ=255)
                                                                                                        gTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCa                                                                                          <  2:16223/76‑1 (MQ=255)
                                                                                                        gTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCa                                                                                          >  1:16223/1‑76 (MQ=255)
                                                                                                            aGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAACGCGAGCcgcg                       <  1:1355588/139‑1 (MQ=255)
                                                                                                            aGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAACGCGAGCcgcg                       <  1:1103975/139‑1 (MQ=255)
                                                                                                             ggACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAgg                                                                                                           >  2:1199461/1‑54 (MQ=255)
                                                                                                             ggACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAgg                                                                                                           <  1:1199461/54‑1 (MQ=255)
                                                                                                                 tGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAACGCGAGCCGCGAAGTg                  <  2:1353140/139‑1 (MQ=255)
                                                                                                                   tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCa                                                                                          <  2:590764/65‑1 (MQ=255)
                                                                                                                   tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCa                                                                                          >  1:590764/1‑65 (MQ=255)
                                                                                                                   tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGt                                        >  1:349902/1‑115 (MQ=255)
                                                                                                                   tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGt                                        <  2:349902/115‑1 (MQ=255)
                                                                                                                   tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACCGGTCGATAACGCGAGCCGCGAAGTGGAg               >  2:1022909/1‑139 (MQ=255)
                                                                                                                        tATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAACGCGAGCCGCGAAGTGGAGTccc           >  1:1270385/1‑139 (MQ=255)
                                                                                                                             aaCTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGaa                                                      <  1:1367804/91‑1 (MQ=255)
                                                                                                                             aaCTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGaa                                                      >  2:1367804/1‑91 (MQ=255)
                                                                                                                              aCTCTGTGGGTATTATGCTTCGCATTAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAACGCGAGCCGCGAAGTGGAGTCCCGTTTTg     <  2:882267/139‑1 (MQ=255)
                                                                                                                                 ttGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAACGCGAGCCGCGAAGTGGAGTCCCGTTTTGACt  >  2:887825/1‑139 (MQ=255)
                                                                                                                                    |                                                                                                                                       
TTTTCATTGTTCCTCCTCACCCGCAACCACTCCGGGTGTTCAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGCGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAACGCGAGCCGCGAAGTGGAGTCCCGTTTTGACT  >  NC_000913/444509‑444776

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