Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,950,557 2 bp→GC coding (451‑452/930 nt) ilvE → branched‑chain amino‑acid aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,950,5570CG100.0% 684.7 / NA 213R151A (CGA→GCA) ilvEbranched‑chain amino‑acid aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (99/114);  total (99/114)
*NC_0009133,950,5580GC99.5% 664.3 / ‑5.9 213R151A (CGA→GCA) ilvEbranched‑chain amino‑acid aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base C (98/114);  total (99/114)
Fisher's exact test for biased strand distribution p-value = 4.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 220 total aligned reads displayed.

GATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCACGACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGT  >  NC_000913/3950466‑3950643
                                                                                           ||                                                                                     
    aTCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCa                                                                                   >  1:3630922/1‑93 (MQ=255)
    aTCGCTGCTGTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCa                                                                                   <  2:3161389/93‑1 (MQ=255)
       gctgctTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAac                                                                                >  1:3793108/1‑93 (MQ=255)
        ctgctTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAaca                                                                               >  2:766991/1‑93 (MQ=255)
         tgctTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAacac                                                                              <  1:453290/93‑1 (MQ=255)
         tgctTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAacac                                                                              >  1:4222280/1‑93 (MQ=255)
          gctTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACAcc                                                                             <  1:3769707/93‑1 (MQ=255)
          gctTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACAcc                                                                             >  2:634028/1‑93 (MQ=255)
           ctTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCa                                                                            <  2:4675144/93‑1 (MQ=255)
           ctTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCa                                                                            <  2:331697/93‑1 (MQ=255)
              tCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATcc                                                                         >  1:1853729/1‑93 (MQ=255)
              tCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATcc                                                                         <  1:4018898/93‑1 (MQ=255)
                ccGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCg                                                                       >  1:3273075/1‑93 (MQ=255)
                 cGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGa                                                                      >  1:3253108/1‑93 (MQ=255)
                  gTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGAc                                                                     >  2:1224717/1‑93 (MQ=255)
                  gTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACAATCCCGag                                                                     >  2:3438152/1‑92 (MQ=255)
                    ggggAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGAcgg                                                                   >  2:3233490/1‑93 (MQ=255)
                         gCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCa                                                              >  1:1797520/1‑93 (MQ=255)
                           gTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCaaa                                                            <  1:3347926/93‑1 (MQ=255)
                            tATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCaaaa                                                           >  1:428489/1‑93 (MQ=255)
                            tATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCaaaa                                                           <  2:3709931/93‑1 (MQ=255)
                            tATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCaaaa                                                           <  2:1719194/93‑1 (MQ=255)
                            tATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCaaaa                                                           >  1:1220855/1‑93 (MQ=255)
                             aTCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAg                                                          >  1:4469600/1‑93 (MQ=255)
                             aTCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAg                                                          <  1:854384/93‑1 (MQ=255)
                              tCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGc                                                         <  1:4566552/93‑1 (MQ=255)
                              tCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGc                                                         <  1:1138614/93‑1 (MQ=255)
                               cTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGcc                                                        <  1:926008/93‑1 (MQ=255)
                                tGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCg                                                       <  1:330135/93‑1 (MQ=255)
                                 gggCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCgg                                                      <  1:2401915/93‑1 (MQ=255)
                                  ggCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCggt                                                     <  1:4171729/93‑1 (MQ=255)
                                   gcgcAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCggtg                                                    <  2:1508311/93‑1 (MQ=255)
                                     gcAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCggtggt                                                  <  2:817427/93‑1 (MQ=255)
                                     gcAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCggtggt                                                  <  1:2588695/93‑1 (MQ=255)
                                     gcAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCggtggt                                                  <  1:3092722/93‑1 (MQ=255)
                                     gcAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCggtggt                                                  >  2:2226349/1‑93 (MQ=255)
                                       aGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTaa                                                <  2:4296306/93‑1 (MQ=255)
                                       aGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTaa                                                >  1:584703/1‑93 (MQ=255)
                                       aGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTaa                                                >  1:1422034/1‑93 (MQ=255)
                                       aGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCGGGAACCGCGCAGCACCAACCGCCATCCCGACGGCGGCAAAAGCCGGTGGTaa                                                <  2:3937953/93‑1 (MQ=255)
                                        gAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAAc                                               >  2:4068268/1‑93 (MQ=255)
                                        gAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAAc                                               <  2:3288365/93‑1 (MQ=255)
                                         aaGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACt                                              <  1:2670748/93‑1 (MQ=255)
                                           gcgcTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAATCACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTAc                                            <  2:1277327/93‑1 (MQ=255)
                                            cgcTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTAcc                                           >  1:1587594/1‑93 (MQ=255)
                                             gcTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACct                                          >  1:2887508/1‑93 (MQ=255)
                                             gcTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACct                                          <  2:1543477/93‑1 (MQ=255)
                                             gcTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGATGGCCAAAGCCGGTGGTAACTACct                                          >  1:1557930/1‑93 (MQ=255)
                                              cTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACctc                                         <  1:768096/93‑1 (MQ=255)
                                               tGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACctct                                        >  2:2465628/1‑93 (MQ=255)
                                                ggAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACctctc                                       <  1:1053495/93‑1 (MQ=255)
                                                ggAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACctctc                                       <  1:3964715/93‑1 (MQ=255)
                                                ggAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACctctc                                       <  1:4642508/93‑1 (MQ=255)
                                                gaatcAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAACCACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACctctc                                       <  2:1769184/89‑1 (MQ=255)
                                                     aGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTccc                                  >  2:2165645/1‑93 (MQ=255)
                                                     aGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTccc                                  <  2:2890317/93‑1 (MQ=255)
                                                     aGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTccc                                  <  1:3391756/93‑1 (MQ=255)
                                                      ggggATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCct                                 <  1:3061242/93‑1 (MQ=255)
                                                      ggggATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCct                                 >  2:2148768/1‑93 (MQ=255)
                                                       gggATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCctg                                >  2:1862917/1‑93 (MQ=255)
                                                         gATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCctgct                              <  1:1677842/93‑1 (MQ=255)
                                                         gATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCctgct                              <  2:1819514/93‑1 (MQ=255)
                                                         gATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGGGGTAACTACCTCTCTTCCctgct                              >  2:1871939/1‑93 (MQ=255)
                                                           tCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCTCTGCtgg                            >  1:3214348/1‑93 (MQ=255)
                                                           tCGATGCGATGGATTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCCGCtgg                            >  2:3574171/1‑93 (MQ=255)
                                                            cgatgcgatgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCtggt                           <  2:3776091/93‑1 (MQ=255)
                                                            cgatgcgatgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCtggt                           <  2:4215466/93‑1 (MQ=255)
                                                              atgcgatgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCtggtgg                         >  1:81548/1‑93 (MQ=255)
                                                              atgcgatgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCtggtgg                         >  1:3575207/1‑93 (MQ=255)
                                                                gcgatgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGt                       <  1:2473064/93‑1 (MQ=255)
                                                                gcgatgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGt                       >  1:3892885/1‑93 (MQ=255)
                                                                  gatgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAg                     <  2:1438595/93‑1 (MQ=255)
                                                                   atgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGc                    >  2:2698361/1‑93 (MQ=255)
                                                                   atgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTAACTCTCTTCCCTGCTGGTGGGTAGc                    <  2:4519138/93‑1 (MQ=255)
                                                                     ggTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGa                  <  2:1015877/93‑1 (MQ=255)
                                                                     ggTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGa                  <  2:2769543/93‑1 (MQ=255)
                                                                      gTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGaa                 <  1:3104862/93‑1 (MQ=255)
                                                                      gTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGaa                 >  1:4313523/1‑93 (MQ=255)
                                                                      gTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGaa                 <  2:1171757/93‑1 (MQ=255)
                                                                       tttCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCACAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAg                >  1:2123091/1‑93 (MQ=255)
                                                                       tttCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAg                >  2:2022485/1‑93 (MQ=255)
                                                                          cctcctGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAgcgc             >  2:3894911/1‑93 (MQ=255)
                                                                            tcctGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAgcgcgc           >  1:3500224/1‑93 (MQ=255)
                                                                            tcctGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAgcgcgc           >  1:519012/1‑93 (MQ=255)
                                                                             cctGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGcgcc          >  2:895485/1‑93 (MQ=255)
                                                                              ctGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGcgccg         >  2:132603/1‑93 (MQ=255)
                                                                              ctGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGcgccc         >  2:3598298/1‑92 (MQ=255)
                                                                               tGGAGCCGCGCAGCACCAAACCCCATCCCGACGGCGGCAAAAGCCGGTTGTAACTGCCGCTCTTCCCGGCTGGTGGGTAGCGAAGCGcgccgc        <  2:722535/93‑1 (MQ=255)
                                                                               tGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGcgccgc        >  2:4032006/1‑93 (MQ=255)
                                                                                ggAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGGGAAGCGCgccgcc       >  2:1664976/1‑93 (MQ=255)
                                                                                ggAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCgccgcc       >  2:4182072/1‑93 (MQ=255)
                                                                                 gAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCa      <  2:2376840/93‑1 (MQ=255)
                                                                                 gAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCa      >  1:1351296/1‑93 (MQ=255)
                                                                                  aaCCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCAc     <  1:4320229/93‑1 (MQ=255)
                                                                                  aaCCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCAc     >  2:4128853/1‑93 (MQ=255)
                                                                                  aaCCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCAc     <  1:3603246/93‑1 (MQ=255)
                                                                                   aCCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACg    >  2:2289274/1‑93 (MQ=255)
                                                                                    ccGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACgg   >  1:3474187/1‑93 (MQ=255)
                                                                                    ccGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACgg   <  2:2589237/93‑1 (MQ=255)
                                                                                    ccGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACgg   <  2:3152105/93‑1 (MQ=255)
                                                                                           ||                                                                                     
GATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCACGACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGT  >  NC_000913/3950466‑3950643

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