Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,954,273 G→T 6.4% E266* (GAA→TAA) ‡ ilvD → dihydroxyacid dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,954,2730GT 6.4% 323.8 / 9.9 93E266* (GAA→TAA) ‡ilvDdihydroxyacid dehydratase
Reads supporting (aligned to +/- strand):  ref base G (35/52);  new base T (3/3);  total (38/55)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

GGGCAACGGCTCGCTGCTGGCAACCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACAGCTGTGTAAAGTTGC  >  NC_000913/3954125‑3954421
                                                                                                                                                    |                                                                                                                                                    
gggCAACGGCTCGCTGCTGGCAACCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTg                                                                                                                                                      <  3:729264/149‑1 (MQ=255)
 ggCAACGGCTCGCTGCTGGCAACCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGAGGCGTTTGa                                                                                                                                                     >  4:636119/1‑149 (MQ=255)
  gCAACGGCTCGCTGCTGGCAACCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGaa                                                                                                                                                    >  1:598644/1‑149 (MQ=255)
     aCGGCTCGCTGCTGGCAACCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAAc                                                                                                                                                 <  1:656704/149‑1 (MQ=255)
         cTCGCTGCTGGCAACCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCa                                                                                                                                             <  2:722189/149‑1 (MQ=255)
          tCGCTGCTGGCAACCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTTCAAACGccgc                                                                                                                                            <  3:543071/149‑3 (MQ=255)
          tCGCTGCTGGCAACCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTTCAAACGccgc                                                                                                                                            <  1:38919/149‑3 (MQ=255)
          tCGCTGCTGGCAACCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAAGTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGCCGTTTGCAATCGcgac                                                                                                                                            <  2:11170/149‑4 (MQ=255)
              tgctgGCAACCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACg                                                                                                                                        >  2:426580/1‑149 (MQ=255)
                        ccACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGGTATcg                                                                                                                              <  1:545870/149‑1 (MQ=255)
                        ccACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATcg                                                                                                                              >  3:895/1‑149 (MQ=255)
                         cACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTTCAAACGccgc                                                                                                                                            >  3:64735/1‑132 (MQ=255)
                         cACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTTCAAACGccgc                                                                                                                                            <  4:64735/134‑3 (MQ=255)
                           cGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTGCGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCCGGATATCGCGa                                                                                                                           <  3:501230/149‑1 (MQ=255)
                            gCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGAt                                                                                                                          <  1:112983/149‑1 (MQ=255)
                            gCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGAt                                                                                                                          >  4:557600/1‑149 (MQ=255)
                              cGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATgg                                                                                                                        <  1:426580/149‑1 (MQ=255)
                                aCCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGt                                                                                                                      <  3:703278/149‑1 (MQ=255)
                                    tAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAAc                                                                                                                                                 <  3:174639/118‑1 (MQ=255)
                                    tAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAAc                                                                                                                                                 >  4:174639/1‑118 (MQ=255)
                                    tAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGAt                                                                                                                  <  3:602180/149‑1 (MQ=255)
                                    tAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGAt                                                                                                                  >  3:726157/1‑149 (MQ=255)
                                       gcagcTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGa                                                                                                               >  2:561960/1‑149 (MQ=255)
                                        cagcTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGAc                                                                                                              <  2:602460/149‑1 (MQ=255)
                                        cagcTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGAc                                                                                                              <  1:648971/149‑1 (MQ=255)
                                         agcTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACt                                                                                                             >  3:220357/1‑149 (MQ=255)
                                            tGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAac                                                                                                          <  2:428335/149‑1 (MQ=255)
                                            tGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAac                                                                                                          <  4:895/149‑1 (MQ=255)
                                             attCCTTAATGCTTGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAaca                                                                                                         >  4:490697/2‑149 (MQ=255)
                                              ttCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAacac                                                                                                        >  2:579873/1‑149 (MQ=255)
                                                 cTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGAc                                                                                                              >  1:6381/1‑140 (MQ=255)
                                                 cTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGAc                                                                                                              <  2:6381/140‑1 (MQ=255)
                                                 cTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGt                                                                                                     >  4:661506/1‑149 (MQ=255)
                                                 cTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGt                                                                                                     >  4:440227/1‑149 (MQ=255)
                                                  ttAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTa                                                                                                    <  4:329370/149‑1 (MQ=255)
                                                  ttAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTa                                                                                                    <  2:375074/149‑1 (MQ=255)
                                                       gCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCt                                                                                                                                      <  2:392375/110‑1 (MQ=255)
                                                       gCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCt                                                                                                                                      >  1:392375/1‑110 (MQ=255)
                                                        cTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCAc                                                                                              >  1:226665/1‑149 (MQ=255)
                                                             aaaCGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACctgct                                                                                         <  3:49724/149‑1 (MQ=255)
                                                                 gCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGt                                                                                                                      >  1:671040/1‑116 (MQ=255)
                                                                 gCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGt                                                                                                                      <  2:671040/116‑1 (MQ=255)
                                                                  cATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTggcgg                                                                                    <  4:726157/149‑1 (MQ=255)
                                                                       ttgAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGgcggcggcg                                                                               >  1:143691/1‑149 (MQ=255)
                                                                         gAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCa                                                                             <  2:579169/149‑1 (MQ=255)
                                                                           aTTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAgg                                                                           >  3:401588/1‑149 (MQ=255)
                                                                           aTTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAgg                                                                           <  4:390974/149‑1 (MQ=255)
                                                                               aCCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGc                                                                       >  2:301913/1‑149 (MQ=255)
                                                                                ccAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCg                                                                      <  3:738812/149‑1 (MQ=255)
                                                                                   aaCGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGaa                                                                   <  1:118348/149‑1 (MQ=255)
                                                                                   aaCGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGaa                                                                   <  1:249219/149‑1 (MQ=255)
                                                                                    aCGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGaaa                                                                  <  1:270983/149‑1 (MQ=255)
                                                                                    aCGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGaaa                                                                  <  1:40270/149‑1 (MQ=255)
                                                                                      gTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATc                                                                <  2:270249/149‑1 (MQ=255)
                                                                                       ttattaCGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCg                                                               >  1:541372/1‑149 (MQ=255)
                                                                                        tattaCGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGCAATCGa                                                              <  2:547145/149‑1 (MQ=255)
                                                                                        tattaCGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGa                                                              <  1:579873/149‑1 (MQ=255)
                                                                                        tattaCGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGa                                                              <  3:576275/149‑1 (MQ=255)
                                                                                         attaCGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGAc                                                             <  1:729795/149‑1 (MQ=255)
                                                                                             cGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCa                                                         <  2:445659/149‑1 (MQ=255)
                                                                                             cGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCa                                                         <  1:558869/149‑1 (MQ=255)
                                                                                               aGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGAc                                                             <  1:412195/143‑1 (MQ=255)
                                                                                               aGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGAc                                                             >  2:412195/1‑143 (MQ=255)
                                                                                                gCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCTACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCa                                                      >  4:235397/1‑149 (MQ=255)
                                                                                                   aaaCGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTACCACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGATTTCACCATGa                                                   <  3:48134/149‑1 (MQ=255)
                                                                                                    aaCGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAg                                                  >  1:232983/1‑149 (MQ=255)
                                                                                                     acgacgAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGt                                                 <  2:598644/149‑1 (MQ=255)
                                                                                                      cgacgaAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTg                                                >  4:335385/1‑149 (MQ=255)
                                                                                                        acgaAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGat                                              >  4:149480/1‑149 (MQ=255)
                                                                                                        acgaAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACACCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGat                                              <  3:275111/149‑1 (MQ=255)
                                                                                                          gaAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGatat                                            >  1:35838/1‑149 (MQ=255)
                                                                                                               tGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATa                                       <  3:149480/149‑1 (MQ=255)
                                                                                                                gCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATaa                                      <  2:35838/149‑1 (MQ=255)
                                                                                                                  aCTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGc                                    <  1:301913/149‑1 (MQ=255)
                                                                                                                  aCTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGc                                    <  4:452250/149‑1 (MQ=255)
                                                                                                                    tGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCtt                                  >  4:525953/1‑149 (MQ=255)
                                                                                                                      ccGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTc                                <  1:561960/149‑1 (MQ=255)
                                                                                                                        gcgTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTccc                              >  3:292186/1‑149 (MQ=255)
                                                                                                                            aaTATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGTTGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCaa                          <  3:636119/149‑1 (MQ=255)
                                                                                                                               atCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGt                       >  4:380974/1‑149 (MQ=255)
                                                                                                                                  gCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTcc                    <  3:15397/149‑1 (MQ=255)
                                                                                                                                  gCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTcc                    >  4:413767/1‑149 (MQ=255)
                                                                                                                                   ccAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCca                   <  3:413767/149‑1 (MQ=255)
                                                                                                                                     aGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGAGATCGCGATGGGTGGATCGCCTAACACCATACTCCATCTGCTGGCGGCGGCGCAGGAAGCGGATATCGACTTCACCATGAGTGATCTCGATAAGCTTTCCCGCAAGGTTCcaca                 <  2:551905/149‑1 (MQ=255)
                                                                                                                                      gTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACAg                >  3:212664/1‑149 (MQ=255)
                                                                                                                                          ggcggcGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACAGCtgt            <  3:73615/149‑1 (MQ=255)
                                                                                                                                          catggcgTTTTCAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGTAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGAAAGGTTCCACAGCtgt            >  2:66947/4‑149 (MQ=255)
                                                                                                                                           gcggcgTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACAGCtgtg           >  4:393515/1‑149 (MQ=255)
                                                                                                                                           atggcgTTTTCAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACAGCtgtg           >  2:672592/3‑149 (MQ=255)
                                                                                                                                                 tttGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTcc                    >  3:615793/1‑134 (MQ=255)
                                                                                                                                                 tttGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTcc                    <  4:615793/134‑1 (MQ=255)
                                                                                                                                                 tttGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACAGCTGTGTAAAGt     >  3:31887/1‑149 (MQ=255)
                                                                                                                                                  ttGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACAGCTGTGTAAAGtt    <  3:335385/149‑1 (MQ=255)
                                                                                                                                                  ttGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACAGCTGTGTAAAGtt    <  3:393515/149‑1 (MQ=255)
                                                                                                                                                    gAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGa                                                   <  1:188439/100‑1 (MQ=255)
                                                                                                                                                    gAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGa                                                   >  2:188439/1‑100 (MQ=255)
                                                                                                                                                    gAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACAGCTGTGTAAAGTTgc  >  1:212479/1‑149 (MQ=255)
                                                                                                                                                    gAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACAGCTGTGTAAAGTTgc  >  2:331442/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGGCAACGGCTCGCTGCTGGCAACCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACAGCTGTGTAAAGTTGC  >  NC_000913/3954125‑3954421

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