Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,954,273 G→T 5.8% 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 5.8% 461.1 / 12.8 138E266* (GAA→TAA) ‡ilvDdihydroxyacid dehydratase
Reads supporting (aligned to +/- strand):  ref base G (47/83);  new base T (6/2);  total (53/85)
Fisher's exact test for biased strand distribution p-value = 5.43e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.49e-01

Only 100 of 140 total aligned reads displayed.

CGGGCAACGGCTCGCTGCTGGCAACCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACAGCTGTGTAAAGTTGC  >  NC_000913/3954124‑3954421
                                                                                                                                                     |                                                                                                                                                    
      aCGGCTCGCTGCTGGCAACCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTTTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACg                                                                                                                                                <  2:1066273/150‑1 (MQ=255)
      aCGGCTCGCTGCTGGCAACCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTTCAAACg                                                                                                                                                <  2:1050624/150‑1 (MQ=255)
       cGGCTCGCTGCTGGCAACCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGc                                                                                                                                               >  1:766429/1‑150 (MQ=255)
       cGGCTCGCTGCTGGCAACCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGc                                                                                                                                               <  1:580372/150‑1 (MQ=255)
        ggCTCGCTGCTGGCAACCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAATCGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGcc                                                                                                                                              <  2:38518/150‑1 (MQ=255)
        ggCTCGCTGCTGGCAACCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGcc                                                                                                                                              <  2:1449104/150‑1 (MQ=255)
        ggCTCGCTGCTGGCAACCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGcc                                                                                                                                              <  2:226365/150‑1 (MQ=255)
           tggctgcgGGCAACCCACGCCGACCGTAAGCATCTGTTCCTTAATGCTGGTAAACGCATTGTTGATTTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTTCAAACGccgcc                                                                                                                                           <  2:298678/148‑4 (MQ=255)
             gctgctGGCAACCCACGCCGACCGTAAGCAGCTGATCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGAc                                                                                                                                         <  1:1283333/150‑1 (MQ=255)
               tgctgGCAACCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGc                                                                                                                                       >  1:514875/1‑150 (MQ=255)
                 ctgGCAACCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTg                                                                                                                                     <  1:1179693/150‑1 (MQ=255)
                 ctgGCAACCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTg                                                                                                                                     >  1:1487458/1‑150 (MQ=255)
                  tgGCAACCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTgg                                                                                                                                    <  2:192287/150‑1 (MQ=255)
                   ggCAACCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGa                                                                                                                                   <  2:1528983/150‑1 (MQ=255)
                        cccACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATcg                                                                                                                              >  2:25164/1‑150 (MQ=255)
                        cccACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATcg                                                                                                                              >  2:21684/1‑150 (MQ=255)
                            cGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGAt                                                                                                                          <  2:854617/150‑1 (MQ=255)
                             gCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATg                                                                                                                         <  2:1170143/150‑1 (MQ=255)
                                gACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGt                                                                                                                      <  1:1462289/150‑1 (MQ=255)
                                 aCCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTg                                                                                                                     <  1:958987/150‑1 (MQ=255)
                                 aCCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTg                                                                                                                     <  1:1599273/150‑1 (MQ=255)
                                  ccGTAAGCAGCTGTTCCTTACTGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTgg                                                                                                                    <  1:715433/150‑1 (MQ=255)
                                  ccGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTgg                                                                                                                    <  1:344206/150‑1 (MQ=255)
                                   cGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTACTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGa                                                                                                                   <  1:448506/150‑1 (MQ=255)
                                    gTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGAt                                                                                                                  >  2:501063/1‑150 (MQ=255)
                                      aaGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGTGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCg                                                                                                                <  1:1442344/150‑1 (MQ=255)
                                      aaGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCg                                                                                                                <  1:16096/150‑1 (MQ=255)
                                      aaGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCg                                                                                                                <  2:250304/150‑1 (MQ=255)
                                        gcagcTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGAc                                                                                                              <  1:1415454/150‑1 (MQ=255)
                                        gcagcTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGAc                                                                                                              <  2:360280/150‑1 (MQ=255)
                                        gcagcTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGAc                                                                                                              >  2:418391/1‑150 (MQ=255)
                                        gcagcTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGAc                                                                                                              <  1:931347/150‑1 (MQ=255)
                                        gcagcTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACAGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGAc                                                                                                              <  2:1144508/150‑1 (MQ=255)
                                         cagcTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACt                                                                                                             <  1:682626/150‑1 (MQ=255)
                                         cagcTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACt                                                                                                             >  1:854321/1‑150 (MQ=255)
                                         cagcTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACt                                                                                                             <  1:1482108/150‑1 (MQ=255)
                                         cagcTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACt                                                                                                             <  1:988094/150‑1 (MQ=255)
                                          agcTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTa                                                                                                            <  1:799761/150‑1 (MQ=255)
                                          agcTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTa                                                                                                            <  2:1482344/150‑1 (MQ=255)
                                          agcTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTa                                                                                                            <  2:1694427/150‑1 (MQ=255)
                                          agcTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGGTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTa                                                                                                            <  1:1096430/150‑1 (MQ=255)
                                             tGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAaca                                                                                                         >  1:244315/1‑150 (MQ=255)
                                                       tGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCa                                                                                               <  2:737877/150‑1 (MQ=255)
                                                       tGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCa                                                                                               >  1:1013317/1‑150 (MQ=255)
                                                        gCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCAc                                                                                              >  1:188593/1‑150 (MQ=255)
                                                         cTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCAcc                                                                                             <  2:647644/150‑1 (MQ=255)
                                                              aaaCGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACctgctg                                                                                        <  1:1613705/150‑1 (MQ=255)
                                                              aaaCGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACctgctg                                                                                        <  2:1525332/150‑1 (MQ=255)
                                                                   cATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTggcggc                                                                                   >  1:1222739/1‑150 (MQ=255)
                                                                     ttgttgAATTGACCAAACGTTATTACGAGCAAAAAGACGAAAGTGCACTGCCGCGTACTATCGCCAGTAAGGCGGCGTTGGAAAACGCCATGACGCTGGATATCGCGAGGGGTGGATCGACTAACACCGTACTTCACCTGCTggcggcgg                                                                                 <  2:1221630/150‑1 (MQ=255)
                                                                         tgAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCa                                                                             >  1:1592938/1‑150 (MQ=255)
                                                                          gAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAg                                                                            <  1:382711/150‑1 (MQ=255)
                                                                          gAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAg                                                                            <  2:514249/150‑1 (MQ=255)
                                                                          gAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAg                                                                            >  1:1134181/1‑150 (MQ=255)
                                                                          gAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAg                                                                            <  1:1283485/150‑1 (MQ=255)
                                                                          gAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAACGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAg                                                                            <  2:501743/150‑1 (MQ=255)
                                                                           aaTTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAgg                                                                           >  2:140274/1‑150 (MQ=255)
                                                                               gACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGc                                                                       <  1:1380535/150‑1 (MQ=255)
                                                                               gACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGc                                                                       >  1:1266872/1‑150 (MQ=255)
                                                                                aCCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAcgcg                                                                      >  2:843191/1‑150 (MQ=255)
                                                                                aCCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCg                                                                      <  1:1523884/150‑1 (MQ=255)
                                                                                 ccAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCgg                                                                     <  1:608049/150‑1 (MQ=255)
                                                                                    aaCGTTATTACGAGCAAAACGACTAAAGTGCACTGCCGCGTAATATCGCCAGTAATGCGGCGTTTGAAAACGCCATGACGCTGGATAGCGCGAGGGGAGGATCGCATAACACCCTACTTCACCTCCGTGAGGCGGCGCAGGCAGCGGaaa                                                                  >  2:311027/1‑150 (MQ=255)
                                                                                      cGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATc                                                                >  2:394323/1‑150 (MQ=255)
                                                                                      cGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATc                                                                <  1:1410432/150‑1 (MQ=255)
                                                                                          attaCGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACt                                                            <  2:766429/150‑1 (MQ=255)
                                                                                           ttaCGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACtt                                                           <  1:1001342/150‑1 (MQ=255)
                                                                                                  cAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGGTGGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATg                                                    <  2:57904/150‑1 (MQ=255)
                                                                                                     aaCGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTTGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGt                                                 <  2:112354/150‑1 (MQ=255)
                                                                                                      acgacgAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTg                                                >  1:838181/1‑150 (MQ=255)
                                                                                                           gaAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATc                                           <  1:1265567/150‑1 (MQ=255)
                                                                                                           gaAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATc                                           <  2:453528/150‑1 (MQ=255)
                                                                                                            aaaGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCg                                          <  2:1502678/150‑1 (MQ=255)
                                                                                                                 gCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAg                                     >  1:569265/1‑150 (MQ=255)
                                                                                                                  cACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGCCTTCACCATGAGTGATATCGATAAGc                                    >  1:998952/1‑150 (MQ=255)
                                                                                                                     ggCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCttt                                 <  2:1268975/149‑1 (MQ=255)
                                                                                                                        cgcgTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTccc                              >  1:856281/1‑150 (MQ=255)
                                                                                                                              atatCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATTTCGCGATGGGTGGATCGACTAACACCGTACTTCCCCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATAGATTAGCTTGCCCGCAAgg                        >  1:310726/1‑150 (MQ=255)
                                                                                                                              atatCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTGAACCATGAGTGATATCGATAAGCTTTCCAGCAAgg                        >  2:1439572/1‑150 (MQ=255)
                                                                                                                                atCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGtt                      <  2:188593/150‑1 (MQ=255)
                                                                                                                                atCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGtt                      >  2:1375948/1‑150 (MQ=255)
                                                                                                                                atCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGtt                      <  2:648756/150‑1 (MQ=255)
                                                                                                                                 tCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTc                     <  1:988660/150‑1 (MQ=255)
                                                                                                                                  cGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTcc                    <  1:301869/150‑1 (MQ=255)
                                                                                                                                    ccAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCcac                  >  2:8581/1‑150 (MQ=255)
                                                                                                                                      aGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACAg                >  1:188891/1‑150 (MQ=255)
                                                                                                                                       gTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACAGc               <  2:806563/150‑1 (MQ=255)
                                                                                                                                          tcatggcgTTTTCAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACAGCtgt            >  2:529096/5‑150 (MQ=255)
                                                                                                                                          aGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACAGCtgt            >  1:1689464/1‑150 (MQ=255)
                                                                                                                                           ggcggcGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACAGCtgtg           <  1:654617/150‑1 (MQ=255)
                                                                                                                                           catggcgTTTTCAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACAGCtgtg           >  1:199338/4‑150 (MQ=255)
                                                                                                                                           catggcgTTTTCAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACAGCtgtg           >  1:1441582/4‑150 (MQ=255)
                                                                                                                                            gcggcgTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACAGCTgtgt          <  2:838181/150‑1 (MQ=255)
                                                                                                                                            atggcgTTTTCAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACAGCTgtgt          >  2:471285/3‑150 (MQ=255)
                                                                                                                                             cggcgTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACAGCTGTGTa         <  1:1144875/150‑1 (MQ=255)
                                                                                                                                              ggcgTTTTCAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACAGCTGTGTaa        >  1:1202385/1‑150 (MQ=255)
                                                                                                                                              ggcgTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACAGCTGTGTaa        >  2:1464448/1‑150 (MQ=255)
                                                                                                                                                 gTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACAGCTGTGTAAAGt     <  1:501063/150‑1 (MQ=255)
                                                                                                                                                    tGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTACCACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACAGCTGTGTAAAGTTgc  >  2:136026/1‑150 (MQ=255)
                                                                                                                                                    tGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACAGCTGTGTAAAGTTgc  >  2:631107/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGGGCAACGGCTCGCTGCTGGCAACCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACAGCTGTGTAAAGTTGC  >  NC_000913/3954124‑3954421

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